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FOREWORD 



There has long been, a need in science and engineering for system- 
atic publication of research studies larger in scope than a journal 
article but less ambitious than a finished book. Much valuable work 
of this kind is now published only in a semiprivate way, perhaps as 
a laboratory report, and so may not find its proper place in the 
literature of the field. The present contribution is the eighth of 
the Technology Press Research Monographs, which we hope will 
make selected timely and important research studies readily ac- 
cessible to libraries and to the independent worker. 



J. A. Stratton 
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PREFACE 

In. 1956 J. M. Manley and H. E. Rowe of Bell Telephone Labora- 
tories analyzed [60] power flow at various frequencies in a nonlinear 
capacitor. Their conclusions were the now-famous Manley-Rowe 
frequency-power formulas, which express a constraint on the powers 
entering the capacitor at each frequency. Their timing was fortunate, 
for within a year or so, the topic of parametric amplification was re- 
vived and applied to microwaves. Scarcely a paper on parametric 
amplifiers has appeared without mentioning the Manley-Rowe formu- 
las, usually to show that the concept of parametric amplification is 
consistent with them. 

It was clear from Manley and Rowe's work that the formulas also 
held for nonlinear inductors, but no mention was made of other de- 
vices of potential interest for frequency conversion. Although early 
parametric amplifiers [103, 114] were made from gyromagnetic 
material, it was not shown until I960 that such material obeyed the 
general Manley-Rowe formulas. An interesting question is, 'What 
systems obey the Manley-Rowe formulas?' Here 1 shall report on 
an attempt to find necessary and sufficient conditions that a physical 
system obey the formulas. 

My interest in frequency-power formulas started on April 23, 1958, 
when I attended an evening lecture (one in a series on solid state elec- 
tronics, sponsored by the Boston Section of the I.R.E.) by Mr. A. G. 
Fox of Bell Telephone Laboratories, on masers and parametric ampli- 
fiers. The Manley-Rowe formulas were presented in simple terms to 
show the plausibility of parametric amplification. 

It was clear from this talk that sets of coils, coupled together, obey 
the Manley-Rowe formulas. Reasoning that the formulas should hold 
for an induction motor, I decided to build a rota ting -machine para- 
metric amplifier [79] ; later I was to learn that this device had been 
in use for many years under the name, "induction generator. " 

With an analogy shown between microwave parametric devices and 
rotating machines, it was natural to ask how far this analogy could 
be extended. How many devices can be used in this way? Specifically, 
what systems obey the Manley-Rowe frequency- power formulas? 

A check of systems known to obey the formulas revealed that they 
were all conservative, in the classical mechanics sense. That is, 
they all had energy state functions. It was a simple matter to test 
the hypothesis that systems with energy functions obey the formulas: 
not only was this hypothesis correct, but the proof of the formulas 
was simplified by using, at the outset, the energy state function. 

Meanwhile, Prof. Hermann. A. Haus had also been interested in 
extending the Manley-Rowe formulas. His motivation was micro- 
wave amplification; he had already shown that the formulas held 
for a nonlinear electromagnetic medium, [35] and (for signals much 
smaller than the pump) a gyromagnetic medium [35] and a longitu- 
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dinal electron beam [36] . Professor Haus agreed to supervise a 
study to determine necessary and sufficient conditions for the Manley- 
Rowe formulas to hold. (While Professor Haus was away from M. I. T. 
during the 1959-1960 academic year, Professor R. L. Kyhl took over 
as supervisor. ) 

The formulas of Manley and Rowe are not the only frequency -power 
formulas, however. Both Page [74] and Pantell [77] have given fre- 
quency-power formulas relating real power in nonlinear resistors. 
Mauley and Rowe had shown formulas relating real power in nonlinear 
reactances, and reactive power in nonlinear resistors. It was not 
hard to find the fourth type of frequency-power formula, relating re- 
active power in reactances. 

It was discovered that a modification of the formulas of each type 
allowed time-varying systems to be discussed. Distributed systems 
were more difficult, however, because the formulas should relate 
the power flow at the boundary of the system. This could be done for 
some systems, but apparently not for others. The key to this diffi- 
culty proved to be Hamilton's principle, a variational principle of 
classical mechanics. Distributed systems that obey Hamilton's 
principle also obey the Manley-Rowe formulas in such a way that 
the power flow can be evaluated (in a systematic way) at the bound- 
ary of the system. The study of a distributed system then reduces 
to the search for an appropriate Lagrangian. A Lagrangian for ro- 
tational fluid flow was obtained from Professor C. -C. Lin of the 
Mathematics Department in the Spring of 1959, and a suitable La- 
grangian, for gyromagnetic media was offered by Professor Alan JL. 
McWhorter early in I960. 

The importance and simplicity of the Manley-Rowe formulas has 
led other workers to extend them. Many authors have discussed the 
formulas, interpreting them in new ways, proving them by new meth- 
ods, and extending them. Of these, we mention S. Duinker [18] and 
P. A. Sturrock [99] f each of whom independently discovered that 
systems described by Hamilton 1 s canonical equations of motion obey 
the Manley-Rowe formulas, 

I am pleased to acknowledge the great help given me by many people 
at M.I. T., including Professors D. C. White, R. L,. Kyhl, P. A. Miles, 
L. J. Chu, H. H. Woodson, and H. A. Haus. The work reported here 
was done in the Energy Conversion Group at M.I. T. because of the 
liberal view of "energy conversion" held by those in charge, notably 
Professors White and Woodson. 

This book is based on a thesis submitted in partial fulfillment of 
the requirements of the degree of Doctor of Science in the Depart- 
ment of Electrical Engineering at the Massachusetts Institute of 
Technology in June, I960. 



August I960 Paul Penfield, Jr. 
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PART ONE 
THE FORMULAS 

Chapter 1 
INTRODUCTION 



This book deals with frequency -power formulas, which constrain 
weighted sums of real or reactive powers entering a device at vari- 
ous frequencies, to be either zero, or positive, or negative. Al- 
though there are four different types of formulas, the most impor- 
tant are the Manley-Rowe formulas, which set a weighted sum of 
real powers at various frequencies entering a nonlinear, time-varying 
reactance, to zero. 

A concept of fundamental importance is that of "power at each 
frequency. " This concept is not difficult to grasp when applied to 
electrical devices, but it is more subtle for acoustic systems, elec- 
tron beams, etc. Suppose we have an electrical terminal pair with 
voltage e(t) and current i(t) in the steady state. Here, as in the 
future, we mean by "steady -state operation" that the variables can 
be written as multiple Fourier series of the form 



e(t) = e + ReS^e (1.1) 

and 

i(t) = i +ReS a i Q e a (1.2) 

where a is an index. The subscripts a on the variables i^, e^ 
indicate the peak values of the variables at frequency o> a , given by 
formulas like 

i a = 2<i(t)e ^ y (1.3) 

where the brackets denote a time average. The sums in Eqs. 1. 1 
and 1.2 include each frequency only once, not each frequency and 
its negative. 

The time -aver age power through the terminal pair, P , is the 
average value of the product 

e(t)i(t) d* 4 ) 

and hence can be written in the form 
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where 



P o = Vo 



is called the M d-c power, " and where 



P a = 



(1-6) 



(1.7) 



is called the "power at frequency ^ a , n the asterisk indicating the 
complex conjugate. Notice that the instantaneous power P(t) = e(t)i(t) 
can be written in a Fourier series like Eq. 1.1, but P , as defined 
in Eq. 1.6, is not the average value of P(t), nor is P a , as defined 
in Eq. 1.7, a Fourier coefficient of P(t). 

This decomposition of power into frequency "components" is quite 
common with electrical circuits, because the expression for instan- 
taneous power, Eq. 1.4, is expressed as a well-defined product of 
two variables. Consider, however, the power-flow vector for an 
electron beam 



= IT X H + f mpV(v v) 



(1.8) 



where E is the electric field, H the magnetic field, m the elec- 
tronic mass, p the electron density, and "v the velocity. The 
second term of this formula, the mechanical part, does not appear 
as a simple product, so we are not sure how to break it apart. It 
could be broken into the two factors 2 and mv(v * V), or other ways. 
It is not obvious which of these ways (if any) has physical significance. 

Similar reasoning applies in defining the reactive power at each 
frequency, Q a . 

It is a property of linear time- invariant elements that no frequency 
conversion takes place in them. We are familiar with the classifi- 
cations of such devices in terms of the real and reactive power at 
each frequency: 



Lossless 


P Q = O 


Pure Lossy 


Q a = 


Passive 


p a 


Inductive 


Q a i 


Capacitive 


Q Q <o 



(1.9) 



These relations are conservation theorems for linear, time -invariant 
devices. 

Such conservation theorems do not apply to nonlinear or time-- 
varying devices. We cannot conclude from the fact that a nonlinear 
capacitor is lossless that all P a vanish. Such a device is useful 
as a frequency converter only because we can put power in at some 
frequencies and get it out at other frequencies. The conservation 
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theorems of Eq. 1.9 no longer hold, but might there not be milder 
conservation theorems that remain in the nonlinear case? There 
are: they are precisely the frequency-power formulas, of which 
the most famous are the Manley-Rowe equations. 

The Manley-Rowe formulas [60] state that if the excitation of a 
nonlinear capacitor is such that the current and voltage have fre- 
quency components of the form mwj + nco Q , where m and n are 
integers, then 

00 on 

nP. 



1 1 



mn 
m=^oo n=l 



and 

i i = > 



m=l n=-oo 

where P mn is the power input at frequency inoOj + no> . The result 
was originally proved for a nonlinear capacitor, but the extension 
to a nonlinear inductor was obvious. The nonlinear capacitor is 
lossless, so the sum of all P mrL must vanish. This is predicted 
by multiplying Eq. 1.10 by w , multiplying Eq. 1.11 by Wj , and 
adding. The Manley-Rowe formulas thus consist of the law of con- 
servation of energy and one additional relation. 

If the capacitor is terminated so that power flows only at fre- 
quencies to x , o> , and <>_ = w a - o> , the summations of Eqs. 1.10 
and 1.11 reduce to 

P P P 

i = =. = - L (1.12) 

<o CO- ttj 

and so the ratio of any two of the powers is the corresponding ratio 
of frequencies. Equations 1.12 are often used to show the plausi- 
bility of parametric amplification, in which power enters at fre- 
quency oij , and leaves at the other two frequencies. 

These formulas also describe basically different devices. Weiss 
[115] has pointed out the similarity to the relations for a simple 
model of a three -level solid-state maser, which are 

"P "P T> 

1 =11 = - 2. (1.13) 

V l V 2 V 3 

where each P is the power input at the corresponding frequency v t 
and where v^ + v 2 = v 3 . In addition, a differential gear described 
by the relation of the three shaft frequencies 

co, = r(o> 3 - c^) (1.14) 

satisfies the torque relations [29] 
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(1.15) 

which can be expressed in terms of the power inputs at each shaft 
as 

"P "p "p 

Furthermore, the theory of the balanced n-phase induction motor 
[116, Sec. 3.6.4; 25, Art. 9.2] predicts that the mechanical power 
output -P m is, in the absence of stator resistance, a fraction (1- s) 
of the stator power input P g . The remainder, sP g , is dissipated 
in the rotor. But since s is the slip ^, these results can be written 

P P P 

m = = - (1.17) 

w s " >r w r "s 

The similarity among Eqs. 1.12, 1.13, 1.16, and 1.17 is striking. It 
suggests that the Manley-Rowe formulas are still more general. We 
might ask how far the formulas can be extended what devices obey 
them? 

The motivation for asking this is strong. In the first place, people 
are interested in building practical frequency converters out of de- 
vices other than the simple nonlinear capacitor, for example, from 
ferrites [103, 114] , electron beams [6, 56, 1, 55] , and plasmas [47] . 
The Manley-Rowe formulas say something of importance about fre- 
quency conversion, so we naturally want to know what systems obey 
them. 

Another motivation lies in the field of energy conversion. Although 
frequency -power formulas relate power at various frequencies, they 
are useful for energy converters in which one can identify the port 
through which power flows from its frequency. As noted before, 
balanced induction motors obey a limited form of the equations. Is 
it possible that unbalanced induction motors obey the more general 
form? What energy conversion devices can be used to circumvent 
the limitations on efficiency predicted by the Manley-Rowe formu- 
las? What devices cannot be used? 

Still another motivation lies in the application of the formulas to 
problems in hydrodynamic and magnetohydro dynamic stability. We 
cite two examples: Plasma confinement by RF fields appears quite 
promising [5, 113, 119]> but as Haus [37] and Rostoker [88] have pointed 
out, there is a possibility of a parametric type of instability. And 
secondly, the excitation of cross waves by larger-amplitude water 
waves [44] is an old [24] phenomenon undoubtedly due to parametric 
action. 

Several other physical systems have displayed similar phenom- 
ena [71]. Parametric excitation has been observed in. string vibra- 
tions [65, 82, 83] and pendulums [7, Sec. 45], and quite recently para- 
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metric energy conversion has been reported between electric and 
acoustic modes of yttrium iron garnet [98]. The interpretation of 
these phenomena would be clearer if the frequency-power formulas 
were known. 

Several extensions of the Manley-Rowe formulas have been made. 
Duinker[l9] has given an extension of a limited form of the formulas 
to a lumped, finite, reactive, electrical network. Rowe [89] has dis- 
cussed an extension to a single time-varying reactance. 

Haus has proved the formulas for some distributed systems of 
practical interest the electromagnetic field in a nonlinear, dis- 
persion-free, inhomogeneous, anisotropic, stationary, reciprocal 
medium [35]; the small-signal perturbation of a heavily -pumped 
gyromagnetic medium [35]; and the small-signal perturbation of a 
parametrically-pumped, longitudinal electron beam[36] . 

The broadest extension to date is that of Page [73] , which shows 
that the formulas hold for any lossless system. Unfortunately this 
result is not true, as the counterexample of Appendix A shows. One 
must use a knowledge of the system beyond its losslessness to show 
the formulas to be valid. Some lossless devices that do not obey 
the formulas are the ideal rectifier, the switch, and the variable 
ideal transformer. 

It is interesting that losslessness is not a sufficient condition for 
the formulas to hold. One would think that it would at least be nec- 
essary, but as the time-varying example of Sec. 4.4 shows, this is 
not true either. 

We have attempted to find necessary and sufficient conditions for 
the Manley-Rowe formulas to hold for a system. To our knowledge, 
no non- trivial necessary conditions are known. However, any system 
that is described by an energy state function obeys the formulas. 
This sufficient condition is broad enough to cover cases of practi- 
cal interest, as shown in Chapter Z. 

The Manley-Rowe formulas are not the only type of frequency- 
power formulas, although they seem to be the most important. Mauley 
and Rowe [60] have given another set of formulas, relating reactive 
power at various frequencies in a nonlinear resistor. These are 



and 



where each Q mn is the reactive power into the resistor at fre- 
quency mjtoj + no> . These formulas have not been very useful in the 
past. 

More useful formulas are those that relate real power at various 
frequencies in nonlinear resistors. Such formulas have been given by 
Page [74, 75] and Pantell [77]. They are inequalities, rather than 
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equalities; for example, PantelTs formulas are 



and 



m n mn 



>0 




(1.20) 



(1.21) 



They hold for nonlinear resistors for which the incremental resist- 
ance 



R - 

R i - di 



(1.22) 



is nonnegative. We call such inequalities formulas of Type III, 
whereas the Manley-Rowe formulas are Type I, and the formulas 
relating reactive power in nonlinear resistors are Type II. 

Thus far we have seen formulas corresponding to three of the 
four boxes in the chart "below 





Real Power 


Reactive Power 


nonlinear 
reactance 


Type I 
Manley and Rowe 
Equalities 




nonlinear 
resistance 


Type HI 
Page 
Pantell 
Inequalities 


Type II 
Manley and Rowe 
Equalities 



The remaining formulas, of Type IV, relate reactive power in non- 
linear reactances; they are derived and discussed in Chapter 3. 

In Chapter 2 we show that the energy-function method of proving 
the Manley-Rowe formulas is the logical generalization of both the 
previously known methods. In Chapter 3 the frequency -power 
formulas of the four types are proved for devices that are not only 
nonlinear, but also time -vary ing. One such time-varying device 
of importance is the switch, which obeys formulas of Types II and 

m. 

The succeeding three chapters constitute Part II of the book. They 
are concerned with specific systems. Lumped systems that obey 
the first and/or fourth types of formulas are discussed in Chapter 
4, and lumped systems that obey the second and third types are 
discussed in Chapter 5. Distributed systems present a special 
challenge, because we want to evaluate the expressions for power 
on the boundary of the system at the ports. In Chapter 6 it is 
shown that this is possible for many systems of practical interest, 
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including all distributed systems that obey Hamilton's principle. 

Some applications of the frequency-power formulas are given in 
Part III. The formulas are conservation principles, and are of the 
most value when used as such as an aid to thinking and under- 
standing, as a check on analytic results, as a tool for finding fun- 
damental limits, etc. Nevertheless, there are problems in which 
the frequency-power formulas can save many steps, and sometimes 
even lead directly to the desired results. Generally speaking, this 
happens when the number of frequencies present is small. Some 
applications of this sort to rotating machines are discussed in Chap- 
ter 7, and some communications applications in Chapter 8. We can- 
not pretend to discuss these applications exhaustively. 

Frequency-power formulas of each type, under a variety of fre- 
quency constraints, are collected for reference in Appendix B, along 
with a partial list of physical devices known to obey them. In Ap- 
pendix C it is shown that for the purpose of applying formulas of 
Type I (Manley-Rowe formulas), one may lump together power at 
all harmonics of a given frequency, considering it to be power at 
the fundamental (or at any one of the harmonics). In Appendix D it 
is shown that for a rotating shaft, one can consider the d-c power 
to be power at a frequency equal to the average shaft speed. And in 
Appendix E it is shown under what conditions new information is ob- 
tained by using a different choice of independent frequencies. 

To help guide the reader, we can .say that the most important 
points to be made are: 

(1) The Manley-Rowe formulas hold for any physical sys- 
tem that has an energy- state function. 

(2) Distributed systems that obey Hamilton 1 s principle also 
obey the Manley-Rowe formulas. 

(3) There are four types of frequency-power formulas, and 
they ought to have equal theoretical standing, in spite 
of the fact that, of the four, the Manley-Rowe formulas 
(Type I) have the most practical importance. 

(4) Frequency-power formulas of each type apply (in a re- 
stricted form) to time -vary ing devices, as well as to 
time -invariant nonlinear devices. 

(5) Frequency- power formulas are conservation principles, 
and are of the most value when used as such. 



Chapter 2 
THE ENERGY -FUNCTION METHOD 

Only three essentially different ways of proving the Manley-Rowe 
formulas, Eqs. 1.10 and 1.11, have been published. One of these is 
applicable only to linear time -varying reactive systems; so there 
are two methods for a nonlinear system. Both of these methods can 
be generalized to multiport nonlinear systems, but the natural ex- 
tension of each uses conditions that are sufficient to enable one to 
define an energy state function for the system. In Chapter 3 we 
prove the formulas by using the energy state function at the outset; 
this method can be considered as a logical generalization of both 
the previously known methods. 

The state function method is preferable because it is general, yet 
simple; because it provides a relatively simple criterion (the ex- 
istence of an energy state function); and because time-varying sys- 
tems are treated very easily. 

The two previous methods of proof of the Manley-Rowe formulas 
are discussed in the first two sections of this chapter, and it is shown 
in the third section that the logical extensions of these methods lead 
to an energy state function for the system. The actual use of this 
energy function is shown in Chapter 3. 

2.1 The Manley-Rowe Formulas 

The Manley-Rowe formulas were originally proved for a non- 
linear time -invariant capacitor, but by simple analogy it was known 
the results held for other single energy- storage elements. It is 
assumed (for the full equations) that there is no loss or hysteresis, 
and that the voltage e is some function of the charge q 



(2.1) 

It is further assumed that the excitation is such that power flows 
into or out of the capacitor only at frequencies of the form 
mc^ + nw Q . We call P mn the power that flows into the capacitor 
at frequency mcjj -f nco . ' 

Because of the nonlinearity, there can be net power into the ca- 
pacitor at some frequencies, and out at other frequencies, in spite 
of the fact that the total power input S m S n P mix vanishes. There is, 
however, one constraint caused by the fact that o> and w can be 
specified independently. The easiest way to state this constraint is 
to write it together with the law of conservation of energy, in the form 

8 
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mP mn 

= (2.2) 



and 



4- nco Q 



nP 
S S _ 22_ = (2.3) 



These are the Manley-Rowe formulas. 

To prove these formulas by any method, we write the voltage in 
the form of a multiply-periodic Fourier series 



where Wj and OJ Q are independent frequencies. To say that two 
frequencies are independent, however, is equivalent to saying that 
their two phases are independent, so let us re -write Eq. 2.4 in the 
form 

where 

^ = Wjt (2.6) 

and 

4> = t (2.7) 

are independent phases. We also represent in this form the charge 
q and the current i : 



and 

^^^o) (2 . 9) 



where the current and charge Fourier coefficients I mn and 
are related by 

(2JO) 



because the current is the time -derivative of the charge. In Eqs. 2.2, 
2.3, 2.4, 2.5, 2.8, and 2.9, the summations over m and n are in- 
tended to include each frequency, but not its negative. Thus if the 
integers (m, n) are included in the sum, the set of integers (-m, -n) 
is not. Furthermore, the sum is not intended to include the d-c term, 
which is shown separately, so (m = 0, n = 0) is not included. A 
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simple way to insure that the summation is taken properly is to sum 
only over positive frequencies, or over one of the indices (m or n) 
from to oo, and the other from -oo to oo, omitting (m = 0, n = 0). 

The Fourier coefficients like E ma are the peak values of the com- 
ponents of the variables at frequency mw 1 + nw Q . They are evaluated 
by averages of the form 



*W = 2 

where the asterisk indicates the complex conjugate, and where the 
brackets indicate an average with respect to cj^ and 4> . Note that 
we consider the voltage and the charge as functions of the phases, 
through Eqs. 2.5 and 2.8. 

The Manley-Rowe formulas can now be proven easily, by the method 
originally used [60] (rather, a slight modification of it). We calculate 
from Eq. 2.8 



D^ ^ v in ^Jv" 1 ! T I1C P ; /-> io\ 

= Re Z 2/ imU-^^e * u \6.L) 

m n mn 



and so 



by Eq. 2.11. The left-hand side of Eq. 2.13 vanishes when averaged 
with respect to 4> x , for it becomes the integral of e dq , which is 
a perfect differential, evaluated at <j>j and c^ + 2-rr, or for identi- 
cal values of q . Thus the right-hand side can be equated to zero 
in the form 

mP 
S S - 2- = o (2.14) 



where 



is the power input to the capacitor at frequency mWj + nco Q . This is 
one of the Manley-Rowe formulas; the other is obtained in a similar 
way. 

"We have shown here the method originally used by Manley and 
Rowe [60] to prove their formulas. In Sec. 2.2 we discuss the 
other important method of proof, and in Sec. 2.3 we show that the 
energy-function method of Chapter 3 can. be considered the logical 
generalization of both of them. 

2, 2 Alternate Methods 

Essentially only two alternate methods of proving the Manley- 
Rowe formulas have been given. Weiss [115] has shown that the 
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equations also hold for a simplified model of a three-level, solid- 
state maser, but his argument involves the quantum nature of the 
maser process, and is not a method of proving the Manley-Rowe 
formulas. 

One of the alternate methods was used by Duinker [19] in his study 
of linear time-varying reactances for frequency conversion. He did 
not obtain both Manley-Rowe formulas, but only one. By calculating 
the actual impedance matrix of a reactive time -varying network, 
Duinker found that it satisfied one formula of the form of Eq. 2. 2. 
He postulated (correctly) that the full form of the Manley-Rowe 
formulas, both Eqs. 2.2 and 2.3, should hold for a network of non- 
linear reactive elements, but was not able to show this by his method. 
(This fact is shown in Sec. 4.10, and was also derived by Duinker 
himself [18]. Since this method is not suitable for nonlinear devices, 
we will not discuss it further. 

The other method of proof has been used by Page [73], Salzberg 
[90], Carroll [12], Clavier [16], and C alien [10] . We write the law 
of conservation of energy 



in the form 

mP rnn nP mn 

co 2 2 - + w n S S - = (2.17) 

i m n mcoj + no) Q o m n mo^ + no> 

where we have used only the fact that the nonlinear capacitor is 
lossless. It is clear that the Manley-Rowe formulas are valid pro- 
vided we set each of the two sums individually to zero. 

Page [73] reasoned that since to and w are independent fre- 
quencies, they can be varied independently, and therefore the factor 
that multiplies each in Eq. 2.17 must vanish. Thus all lossless de- 
vices should obey the Manley-Rowe formulas. This, however, over- 
looks the possibility that the suras might themselves depend on w 1 
and o> : when one of these is varied, the sums might vary also. 
Some lossless devices that do not obey the Manley-Rowe formulas 
are the ideal diode, the switch, the variable ideal transformer, and 
their mechanical counterparts, including the ideal ratchet and the 
periodically -applied clamp. 

Salzberg[90], on the other hand, noted that the ratio of P mn to 
mco + nw Q is, from Eq. 2.15, 



mco 



(2.18) 



He then stated that the Q are functions of the E mn and the 
nonlinear characteristic of the capacitor, but not of co a or co . Al- 
though the frequency of which Q mu is the Fourier component does 
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depend on c^ and w , the value of Q mn does not. Salzberg then 
was able to conclude that the only admissible solution of Eq. 2.17 for 
various values of w a and OO Q is with each sum individually vanish- 
ing. Although the expressions in his paper [90] are for quite simple 
cases, with only a few frequencies, his method (as he pointed out) 
is also used to obtain the full formulas. 

This method has been given more recently, but without signifi- 
cant improvement. Carroll [12] used essentially Salzberg's rea- 
soning, but in more mathematical form. Clavier [16] derived Eq. 
2.17 and then claimed that this implied Eqs. 2.2 and 2.3, but unfor- 
tunately did not disclose his reasoning in detail. C alien [10] also 
derived Eq. 2.17 and demonstrated with a simple example that Salz- 
berg 1 s reasoning is correct the Q mn and E mn are related by 
the capacitor, independent of co t and W Q . 

Manley and Rowe [59] have commented on this method, pointing 
out (correctly) that going from Eq. 2.17 to Eqs. 2.2 and 2.3 does 
require knowledge of the device beyond its losslessness. They con- 
sider an ensemble of systems with different coj and o> , but the 
same E mn and Q mn - 

As Manley and Rowe [59] have stated, this alternate method does 
not have any analytical advantage over the original method. The con- 
tribution, then, is one of providing insight, not rigor. 

2,3 Extension of These Methods 

Suppose we wish to extend the method of Manley and Rowe, as 
given in Sec. 2.1, to multiport devices. To demonstrate the logi- 
cal extension, let us consider two coils in a nonlinear medium, with 
coupling. We suppose there is no winding resistance, hysteresis, 
eddy currents, or other loss. The constitutive relations of the 
system then are in the form 

i t = i^Xj.V,) (2.19) 

and 



where the subscripts 1 and z refer to the coils, and where the X 
variables are the flux linkages, the time-integrals of the voltages. 
"We write each of the variables in Fourier series 

( ^ + **J (2.21) 

X 2 (t) = X 20 + Re 2 m VW (mt| ' 1 + n *o> ( 2 .22) 



i2 (t)= i20 +Re Sm Z nWn e^ + *o) <"*> 
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and 



where we show each variable as a function of the phases <j> and cj> . 
Since the voltages are the time -derivatives of the charges, each 
E imn an(i E zmn is J( mcja i + nc V times the corresponding A Imn or 



Just as in Sec. Z.I we form the partial derivatives 



and 



and take the averages 



and 



using expressions for I * and I * similar to Eq. 2.11. We 
r imn mn ^ L 



may then add Eqs. Z.29 and 2.30, obtaining 



mP mn 



where 



is the power input to the device at frequency mw 4- nw Q . 

To show that the Manley-Rowe formulas hold, we want to show 
that the left-hand side of Eq. 2.31 vanishes. At this point in Sec. 
2.1, we noted that the average with respect to 4> 2 of the corres- 
ponding quantity vanished, since the quantity edq was a perfect 
differential. Here we need to show that \^ l + i 2 dX 2 is a perfect 
differential. This is true if and only if 
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'- = * 2 *' 2 (2.33) 



8X 2 



The logical extension of the proof of Sec. 2.1 then requires that 

Eq. 2.33 hold. 

But Eq. 2. 33 is a condition that is sufficient [42, Section 6. 10; 80] 
to enable us to define an energy state function (the stored magnetic 
energy) by the line integral in the Xj -X 2 plane 

dX, (2.34) 



The proof of the Manley-Rowe formulas is simplified by working 
with the energy state function from the beginning. This "energy func- 
tion method" is given in Chapter 3; it is seen to be a logical exten- 
sion of the original proof of Manley and Rowe [60] . 

Now suppose that we wish to extend the alternate proof given in 
Sec. 2.2 to multiport devices. We shall see that this involves as- 
sumptions that are also sufficient to enable us to define an energy 
state function. 

We discuss the same device the nonlinearly- coupled coils, with 
the currents i and i functions of the flux linkages X x and X 2 . 
Then, because the system is lossless, we write 

mP mn nP mn , M 

"1 S tn 2 n moil + na> + % S m S n ^ + n% = < 2 ' 35 > 

where 

P 

1X3111 rfc 

mco -f nw imn imn zmn 2 mn' (2.36) 

Suppose the system is operating with given A imn , A 2tnn , X 10 , 
X 20 , 6>j f and o> Q . Then in particular, \ l (t) and X 2 (t) are de- 
termined, and hence so are i (t) and i 2 (t) , and hence so are I imn 
and I . Now suppose that u> is changed slightly, keeping X 10 , 
X 20> A jmn , and A zmn the same. Does this change affect the I Jinn 
or I^mn at all? It does not. By definition I jmn or I 2mn are in- 
dependent of time. But in the specification of \ l (t) and X 2 (t) , and 
therefore in the specification of i (t) and i z (t), co appears only 
multiplied by the time t . In any functional dependence coj can only 
appear multiplied by t , and therefore any dependence of Ii mil or 
Izmn must contain t if it is to contain o>i . Thus the Iimn anc * -^mn 
cannot depend on QX , but only on the Fourier components of the flux 
linkages. It is true that the frequency of which Ii mn is the Fourier 
coefficient changes, but the value of Ii m n does not. 

It is then clear from Eq. 2. 36 that the sums of Eq. 2. 35 are inde- 
pendent of Q! . If we change a>j slightly, the left-hand side can van- 
ish only if the first summation itself vanishes; Thus we have shown 
that each sum individually vanishes, or the Manley-Rowe formulas 
hold. 
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Notice that this proof, an extension of the method of Salzberg and 
others [90, 12, 16, 10, 59] requires three assumptions: 

(1) The currents depend at each instant on the flux linkages, 
as expressed in Eqs. 2.19 and 2.ZO. 

(2) The device is lossless. 

(3) The variables \ l (t] and \ 2 (t) are physically independ- 
ent, so that each can be specified as a function of time 
independent of the other. 

This last condition is necessary if we are to vary w 1 without thereby 
affecting the flux linkage Fourier coefficients. 

But these three conditions are sufficient to define an energy state 
function of the form of Eq. 2.34, for as we shall now show, they 
imply Eq. 2.33. 

Let us operate the system around the path in the \ 1 - X 2 plane 
shown in Fig. 2.1. This operation is possible since the variables 




Fig. 2.1. Path of operation in the \ l - X 2 plane. The 
position of the path is arbitrary, so the re- 
sult of Eq. 2.33 holds for any values of \ 1 
and X 

X and X are independent. Since the system is lossless, the net 
energy out over such a cycle must vanish, for otherwise we could 
get time -average power out of the device by merely repeating this 
cycle periodically. The net energy out of the system per cycle is 
the line integral around the path in the \ l - X 2 plane 

(2.37) 



(2.38) 



which is, to second order in the quantity AX , merely 



, i 

7 2 1 

(AX) 2 - 
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shown either by direct calculation or by use of Stoke 's Law. Be- 
cause this vanishes at each point in the \ l - X 2 plane, Eq. 2.33 
must hold. Therefore we may define an energy state function of 
the form of Eq. 2.34. 

Both methods [60, 90] of proving the Manley-Rowe formulas, when 
extended to multiport devices, lead logically to the existence of an 
energy state function. Thus the energy function method we use in 
Chapter 3 can be considered a logical extension of each of the methods 
of proof. 

The energy function method leads to general results with mathe- 
matical simplicity, and provides some insight into what sort of 
systems obey the Manley-Rowe formulas. The criterion of having 
an energy state function is one that is easily applied to specific 
systems. Time-varying systems are merely those that have an 
energy state function with an explicit time dependence; they are 
discussed easily by the energy function method. 



Chapte r 3 
FOUR TYPES OF FREQUENCY -POWER FORMULAS 



In this chapter we shall prove the four types of frequency-power 
formulas. The first type (Manley-Rowe formulas) will be proved 
by an energy function method; the formulas are extended in a vari- 
ety of ways: 

(1) We allow more general frequencies than those of the 

form mco + nto . 
(Z) We allow for more than two independent frequencies. 

(3) We take into account d-c power input (by a modifica- 
tion of the formulas), such as that at the shaft of a 
rotating machine. 

(4) We prove the formulas starting from a simple, gen- 
eral characteristic (the existence of an energy state 
function), and so cover many important physical sys- 
tems. 

(5) We include time-varying systems. 

Of these extensions, some are trivial. It is easy to account for 
more general frequency relations and more than two independent 
frequencies [18, 99, 12, 120]; it is only for reasons of simplicity that 
not all previous writers have done so. We use notation so general 
that it will be no bother for us to make this extension. 

Other extensions are of some importance. The extension to in- 
clude d-c power input has been useful with rotating machines, in 
which d-c mechanical power flows. The extension to other physi- 
cal systems, which we will discuss in Chapters 4 and 6, is also 
important. 

In Sec. 3.1 we lead up to the proofs of the four types of frequency- 
power formulas, by introducing the notation used. The actual 
proofs are contained in Sees. 3.2, 3.3, 3.4, and 3.5. Appendix B 
includes a summary of the four types of formulas, together with 
their forms for some special cases, and a tabulation of some of 
the devices known to obey them. 

3.1 Notation 



Generally speaking, frequency -power formulas of Types I and IV* 
hold for nonlinear time-varying reactances, and frequency-power 
formulas of Types II and III hold for nonlinear time-varying resist- 
ances. Some devices obey formulas of only one type, and some obey 
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all four types, so the preceding statement must not be taken too 
literally. But in general, reactive, energy- storage devices obey 
Types I and IV, and resistive, dissipative devices obey Types II 
and III. The choice of notation in each of the four proofs to follow 
is clear if this is kept in mind. 

Two important notational schemes are used in this chapter. The 
first is for the variables that describe the devices, either reactive 
or dissipative; the second is for the frequencies. 

The Variables. Although the most important applications of fre- 
quency-power formulas are to electrical devices, we keep the no- 
tation for variables general. We suppose there are two variables 
that, when multiplied together, give power entering a device, or 
vector power flow density. Typical variables are voltage e and 
current i; mechanical force f and velocity v^ rotational torque 
T and speed co ; electrochemical potential p. and particle current 
density I; electric and magnetic fields forming Poynting ! s vector 
E"xS; and acoustic pressure IT and velocity v. The general no- 
tation for such variables will be f^ and v^ (taken from the mech- 
anical case), where i is an index. 

Dissipative devices are characterized by relations among these 
variables: a nonlinear resistor has a voltage -cur rent curve; a 
mechanical friction device has a nonlinear force-velocity curve, 
and a dissipative medium would have a current density T (related 
to H through one of Maxwell 1 s equations) a function of Til. Thus we 
say for dissipative systems that 

fi = fi<v k ) (3. 1) 



or that each f , is a function of all the v, . 

If these n constitutive" relations, Eqs-. 3,1, are "reciprocal, " in 
the sense that 



8f (v k ) 

~ (3 - 2) 



then we may define a "dissipation function" by the line integral in 
multidimensional v^ space [42, Sec. 6.10; 80] 



G(v k ) = S f^Vi (3.3) 

where 



In other contexts this function has been called the Rayleigh dissi- 
pation function [30, Sec. 1*5], the content [68], and the co-content 
[68], If Eq. 3.2 does not hold, we cannot define this state function, 
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for the integral of Eq. 3.3 depends on the path chosen. 

Time-varying dissipative systems are those for which the con- 
stitutive relations of Eq. 3.1 are replaced by ones that depend ex- 
plicitly on time: 

f i = f i<v k , t) (3.5) 

If Eq. 3. 2 holds for the time -varying dissipative system, then the 
state function defined in Eq. 3.3 will be time-varying, that is, it 
will have an explicit time dependence. 

A typical device that has a time-varying constitutive relation is 
a variable resistor (rheostat) whose shaft is being moved as a 
function of time, so that the resistance is some function R(t) of 
time. Then the constitutive equation is 

e(i, t) = R(t)i (3.6) 

and the voltage is seen to be a function of the current and time ex- 
plicitly. 

This explicit dependence on time should not be confused with the 
actual dependence on time of the voltage, part of which is caused 
by R(t) , and part by the actual time dependence of the current. 

In Sec. 3.3 we prove frequency-power formulas of Type II by use 
of this dissipation function. In Sec. 3.4, however, in which we 
prove the formulas of Type III, we do not require the existence of 
such a state function, but instead only the functional relations of 
Eq. 3.5 and the fact that the matrix with i-k entry, 



is positive definite [43, Sec. 1.17] . 

Reactive, or energy storage, devices, on the other hand, do not 
have constitutive relations of the form of Eq. 3.5, For example, the 
instantaneous voltage of a capacitor is not a function of the current, 
but instead of the charge on the capacitor plates, which is the inte- 
gral of the current. We discuss reactive devices by defining a new 
set of variables, x^ , as the integral of the v^ variables. The in- 
tent is that the f^ variables should, for reactive devices, be a func- 
tion of the X| variables: 

fj = f (x k ) (3.8) 

Some devices with constitutive relations of this sort are (besides 
the capacitor): the inductor, with current i a function of the flux 
linkage X , whose derivative is the voltage e ; a mechanical spring, 
for which the force f is a function of the displacement x , whose 
time derivative is the velocity v ; and an electromagnetic medium 
in which "E is some function of the displacement vector U, whose 



20 f r e qu e n c y - p o we r formulas 

time -derivative is related (through one of Maxwell's equations) to 
5, and in which H is some function of the magnetic induction vec- 
tor B , whose time -derivative is related to E . 

If these constitutive relations are Ir reciprocal !! in the sense that 



(3-9) 



then the line integral in multidimensional x. space 

U(x.) = S.jXdx. (3.10) 

defines a state function [42, Sec. 6.10; 80], independent of the path of 
integration, where 

* 



Defined in this way, the function U is the energy, since its time 
derivative 



' i 8 Xi dt ' iii ' 

is the power input to the device. 

The treatment of time-varying reactive systems is similar. The 
constitutive relations of Eq. 3.8 are replaced by ones with explicit 
time dependence 



K Eq. 3.9 holds, then the energy function defined in Eq. 3.10 has an 
explicit time dependence, and Eq. 3.12 does not hold. The device 
cannot, therefore, be called lossless, since there can be net gen- 
eration or destruction of energy within the framework of the vari- 
ables considered. This is discussed further in Sec. 4.4, in con- 
nection with the time -varying capacitor. 

Other energy state functions can be derived from U by Legendre 
transformations [30, Sec. 7.1], and in Sec. 3.2 we show how the exist- 
ence of this U function, or any energy state function related to it 
by a Legendre transformation, leads to frequency power formulas 
of Type I, the Manley-Rowe formulas. Energy state functions that 
are not related to U by a Legendre transformation also lead to fre- 
quency-power formulas, but the interpretation of the expressions is 
not always clear. 

In Sec. 3.5 frequency-power formulas of Type IV are proved, for 
reactive devices that obey constitutive relations of the form of Eq. 
3.13, provided that the matrix with i-k entry 

9f 

(3.14) 
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is positive definite [43, Sec. 1.17] . We do not use the state function U. 
In summary, the important points to remember about this notation 
for variables are that products of f. and v. give power flow or powei 
input, and that the x^ variables are the time-integrals of the v- . 

The Frequencies. For generality we do not wish to use only fre- 
quencies of the form mcOj + nu> . Perhaps in practice these frequency 
constraints are sufficient; our analysis is nevertheless more gen- 
eral and simpler, because we choose not to specify what frequencies 
are present. We merely assume that there are a number of fre- 
quencies present, and that only a few of these are independent. The 
results reduce to the formulas given earlier when we specify that 
the frequencies are to be of the form mwj + nco Q with ^ and o> 
independent. They also reduce to other formulas given before [18, 
99,12,120] when the appropriate frequencies are specified. 

We assume that the device is operating in the steady state, so that 
each of the variables can be written in a multiply -periodic Fourier 
series 

jw t 
v.(t) = v io + Re S a (v.) a e a (3.15) 



f.(t) = f io+ ReZ a (f.) a e a (3.16) 



so that 



x i = x io + V + Re S a( x i>a e < 3 ' 17 > 

where the frequencies of the system are oo , a being an index. Be- 
cause each v- is the time derivative of x. , their Fourier coeffi- 
cients are related: 



Because the operation is in the steady state, the state functions (if 
defined) are bounded and can be represented in the same form: 

ju t 
G(t) = G +Re2 a G a e a (3.19) 

and 



_ 
U(t) = U o +Re a U a e a (3.20) 

The set of frequencies co is large enough to include all frequencies 
needed for Eqs. 3.15, 3.16, a 3.17, 3.19, and 3.20. In these equations, 
the summations over a are not intended to include a given frequency 



22 f r e qu ency - po we r formulas 

and its negative, and are not intended to include the d-c terms, which 
are shown separately. Thus, if some frequency GO is included in 
the summation, -co is not. One way to insure this requirement is 
to sum over only positive frequencies, although this method is not 
the only way, nor is it always the best. We will henceforth assume, 
in all summations of this sort, that we do not sum over both co and 
-co, but only over one of them (either one). 

If a negative frequency is included in the summation, then the 
Fourier coefficients associated with this negative frequency are the 
complex conjugates of the Fourier coefficients for the corresponding 
positive frequency. In addition, the real power associated with a 
negative frequency is the same as the real power associated with 
the corresponding positive frequency, whereas the reactive power 
associated with a negative frequency is minus the reactive power 
associated with the corresponding positive frequency. 

In general, the systems we discuss are time-varying. The explicit 
time dependence of the variables arises through specific functions 
of time; these functions can, in the steady state, be represented in 
terms of some frequencies coo , where (3 is an index. The set of 
frequencies (co a ), which includes all frequencies necessary for the 
expansions of Eqs. 3.15, 3.16, 3.17, 3,19, and 3.20, thus probably 
includes many of the co as well as many other frequencies. 

Not all the co frequencies are independent; there are constraints 
that enable us to define a set of independent frequencies, in terms 
of which all co a can be expressed. For example, if the frequencies 
to^ are of the form m^ + na) Q , we may define as independent fre- 
quencies tOj and CO Q , although this is not the only possible choice. 

In general, any particular frequency co a depends on some or all 
of the independent frequencies. The frequencies coo , which go toward 
specifying the explicit time dependence of the system, thus depend 
on some of the independent frequencies, possibly not all. How many 
of the independent frequencies the co g depend on is to some extent 
determined by the independent frequencies and how they are chosen. 
We have considerable freedom in choosing them; thus, with one 
choice we might find that the co frequencies depend on more of the 
independent frequencies than with another choice. Because there are 
frequency-power formulas associated with each independent frequency 
on which no co depend, it is wise to choose the independent fre- 
quencies in such a way that only one (or only a few) is necessary 
to specify all the tog . 

An example may clarify this. Suppose a capacitor is varied mech- 
anically at frequency to and its harmonics, and that a signal of 
frequency co g is put on the electrical terminals. Sidebands will be 
generated with frequencies of the form a> s + nto , where n is a 
positive or negative integer. Considering onlylhe electrical vari- 
ables, the device is represented by a time-varying capacitance, 
with power exchange at frequencies nto + co g . There are two in- 
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dependent frequencies, and we might choose for them oo s and 
oo + w . But if we do, the frequency that specifies the explicit time 
dependence of the system, GO , depends on both independent fre- 
quencies, and there are no frequency-power formulas. On the other 
hand, if we choose as independent frequencies to and oo + <*> g , 
there is a frequency -power formula associated with the independent 
frequency u> + oo g . 

Let us call oo a , where a is an index, the independent frequencies 
on which the coo , which specify the explicit time dependence of the 
system, do not depend. Then all partial derivatives of the form 



(3.21) 



are well defined, and all partial derivatives of the form 

9w (3 

(3.22) 

**a 

vanish. There is a frequency -power formula associated with each 

"a- 

To review this notation: 

(1) Frequencies ca are all frequencies present in the 
Fourier analysis of the actual time dependences of 
the variables. 

(2) Frequencies o>o are all frequencies that help to 
specify the explicit time dependence of the consti- 
tutive relations. 

(3) Frequencies o> a are all independent frequencies 
on which the cog do not depend. 

For simplicity, in the proofs of the formulas we use the phases 

^a = "a 1 * ^B = w * aEL(i ^a = ^a*' The P hases 4>p do not de " 
pend on the $ a , and any finite change a in an independent phase 

<b induces changes 
a 

8"a 

A _ _ n 

" 



in the dependent phases <(> a . 

With the variables given in the form of Eqs. 3.15 and 3.16, the 
time -average power input P can be written as 



where 



= S i P io + S a P a 



, f. (3.24) 

io io 
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is the d-c power input caused by the i-th variables, and where 

P a = iReS^f.)*^ (3.25) 

is the total power input at frequency w a . Frequency-power formu- 
las of Types I and III express relations among the P^ Q and P a . 
Similarly, the reactive power at each frequency oi> a is generally 

Q = ilmSfv (3.26) 



a 



or its negative. Frequency-power formulas of Types II and IV ex- 
press relations among the Q . 

There are four difficulties with interpreting Eq. 3.26 as reactive 
power that do not arise in the interpretation of Eq. 3. 25 as real 
power at frequency co a . First, there are in use two definitions of 
reactive power, each the negative of the other. Guillemin [33, Chap. 
7, Sec. 4] defines the 'Vector power" 

P + JQ = |E*I (3.27) 

where E and I are peak amplitudes of voltage and current re- 
spectively. According to this definition, reactive power flows out 
of inductors and into capacitors. This is the reverse of the cus- 
tomary definition, used, for example, by Manley and Rowe [60] 

P + jQ = | El* (3.28) 

According to this convention, which is the one we use in later chap- 
ters, reactive power flows out of capacitors and into inductors. 
Each scheme is self-consistent, and it is an arbitrary choice that 
must be made. 

The second difficulty with interpreting the Q a of Eq. 3.26 as 
reactive power is that we have not specified what the variables f 
and v stand for. If v stands for the current and f for the volt- 
age, then Q a is the negative of the reactive power as we use it, 
whereas if v stands for the voltage and f for the current, then 
Q a is just equal to the reactive power as we define it. 

The third difficulty is that with multiport devices, the contribu- 
tion to Q^ at one port might be the ordinary reactive power there, 
whereas the contribution at another port might be the negative of 
the reactive power. This happens, for example, in a gyrator. In 
this case, we cannot interpret the Q a of Eq. 3.26 as either the re- 
active power or minus the reactive power. 

The fourth difficulty occurs if some of the frequencies o> a are 
negative. Then the Q a associated with these negative frequencies 
are just the negatives of the Q a that would be associated with the 
corresponding positive frequencies w . 

These four points must be kept in mind in interpreting frequency- 
power formulas of Types n and IV. 

Using the notation developed in this section, we can easily prove 
the four types of frequency-power formulas. 
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3. 2 Type I. Manley-Rowe Formulas 

Here are the purely mathematical steps leading to the Mauley- 
Rowe formulas. We use the notation explained in Sec. 3. 1, and 
start with an energy state function. 

We assume we have a "reactive" device: more precisely, we 
assume that the device has a time-varying energy state function 

U = U( Xi , t) (3.29) 

with variables 

, t) 



Using the frequency notation of Sec. 3.1(The Frequencies), we write 
all f i? X, and U in the form of Eqs. 3.16, 3.17, and 3.20. In 
our applications the v- quantities, which contribute toward d~c 
power input, are in the nature of frequencies, because they appear 
multiplied by t . (For example, the position of a rotating shaft 
would be in the form of Eq. 3.17 with v^ Q the average speed). We 
assume that these v. are, for the purpose of defining independent 
frequencies, included among the o> a frequencies. Thus all quan- 
tities of the form 



8v lo 



(3.31) 



are well defined. 

Now let us make a small change 6<j>~ in one of the independent 
phases, <f> a , keeping all the others fixed. We can compute, to first 
order, the change 6U in the state function U . We find from Eq. 
3.20 



a a 
6U = Re S a jU a e 6<j> a (3.32) 

We have another way of calculating 6U . Because the explicit time 
dependence of U does not change (remember none of the 4> de- 
pend on 4> a ) 6U can be expressed in terms of its dependence on 

the x. variables as 
i 

6U = ^f.fix- (3.33) 

which becomes, from Eqs. 3.16 and 3.17 

6U = S. 



J^o* 1 
f^e a J 
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These two calculations of 6U must agree for all time. Equation. 
3.32 shows 6U as a sum of sinusoids, which vanishes, when av- 
eraged with respect to time. Applying the same test to Eq. 3.34, 
we find that its average value is 




a 

which must vanish. Thus, by use of Eqs. 3.24 and 3.25 
P. 8v. P 8co 



io 8o> a a 



= 0.36) 



for each frequency o> . These are the Manley-Rowe formulas, ex- 
tended to include d-c power input and more general frequency con- 
straints. They reduce to the ordinary form of the Manley-Rowe 
formulas, Eqs. 2.2 and 2.3, when the d-c powers vanish, when the 
device considered is the ordinary nonlinear capacitor, and when the 
frequencies of the system are in the form mcoj + nco , with o^ and 
o> chosen as independent frequencies. 

The d-c terms in Eq. 3.36 have proved useful so far only for ro- 
tating machines, so in quoting these formulas in later chapters the 
d-c terms are usually omitted. In Chapters 4 and 6 we discuss sev- 
eral physical systems that have energy state functions, and there- 
fore obey these formulas. 

3.3 Type II 

Frequency-power formulas of Type II relate reactive power in 
nonlinear dissipative devices. We prove them here with a state 
function argument similar to that used for proving the Manley-Rowe 
formulas in Sec. 3.2. 

We assume the device is "dissipative: 1 ' more precisely, we as- 
sume that it has a time -varying dissipation function 

G(v if t) (3.37) 

such that 

8G(v if t) 



8v k 



(3 ' 38) 



Using the frequency notation of Sec. 3.1 (The Frequencies), we write 
all vj, fj, and G in the form of Eqs. 3.15, 3.16, and 3.19. 

Now let us make a small change 6<j> a in one of the independent 
phases <(> a , keeping all the others fixed. We can compute, to first 
order, the change 6G in the state function G . We find from Eq. 
3.19 
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9oj a j a t 
6G = Re S a JG a e 6 * a (3.39) 

3. 

Notice that 6G consists of sinusoids, and therefore its time -average 
value vanishes. 

But we have another way of calculating 6G . Because the explicit 
time dependence of G does not change (remember that none of the 
<j> s depends on <|> ), 6G can be expressed in terms of its depend- 
ence on the v^ variables as 

6G = S i f i 6v. (3.40) 

which becomes, from Eqs. 3.15 and 3.16 

T j" tl f 8(x) a jw t -I 

SG = S. [f. o + Re S Q (f.) Q e * J JRe S a j(v i ) a e * 6* a J (3.41) 

Because the two calculations for 6G must agree for all time, the 
average value of Eq. 3.41 

800 

S^iReKv^)*^^ (3.42) 

a 

must vanish, and so by use of Eq. 3.26 

9" Q 
SQ = < 3 ' 43 > 



for each frequency u> . These are the frequency -power formulas 
of Type II; some devices that obey them are discussed in Chapter 
5. 

In interpreting these formulas, the remarks following Eq. 3. 26 
must be kept in mind. In particular, the gyrator does not obey Eq. 
3.43, even though one can define a dissipation function G for it; 
this is because the Q that result are neither the reactive powers 
nor their negatives. 

3.4 Type III 

Page [74, 75] and Pantell[77] have given frequency-power formu- 
las relating real power in positive nonlinear resistors. By "posi- 
tive, " it is meant that the incremental resistance 

R i - f < 3 -^ 

is nonnegative. 

It is interesting to note that both Page's and Pantell ! s forms of 
the formulas have advantages. PantelPs results are easily applied 
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to devices with more than one independent frequency, whereas Page's 
formulas are not written directly in terms of the independent fre- 
quencies. On the other hand, when they both apply (as they do, for 
example, to harmonic generators) Page's formulas are more power- 
ful. Here we show the logical generalization of each we retain 
all the power of Page's formulas, but express them in terms of de- 
pendent and independent frequencies to make them easier to apply. 
The method of proof used is a modification of Page's method. 

We assume the device considered is "dissipative, " in that tizne- 
varying constitutive relations relate the f^ and v., variables in the 
form of Eq. 3.5. Using the representations of Eqs. 3.15 and 3.16, 
we compute the change in the variables caused by a finite change 
6 a in one of the independent phases <|> a . This produces a finite 
change 

3oj 

9 a = -5T 9 * (3 ' 45) 

a 

in each dependent phase (j> , so the change in each v^ is 

jo; t J0 
A Vi = Re 2 a (v i ) a e a (e a - 1) (3.46) 

and 

K* J 9 a 
Af = Re S a (f.) a e a (e a - 1) (3.47) 

If we multiply Av- by Af. and sum over i, and then take the time 
average of the sum, we find 

\< 9 a> = 

= |Re S.S a (v i ) a (f i )*(l - cos9^ (3.48) 

and so the quantity h a defined in Eq. 3.48 can be written 

h a( Q a> = S a P a< 1 ' COS V (3 ' 49) 

by use of Eq. 3.25. 

So far we have not proved anything of value, for there is no esti- 
mate of the quantity k a (0 a ) - However, because the f,, are func- 
tions of the v. variables, the changes Afj can be expressed in 
terms of line integrals in v.- space: 



v k 



dv k (3.50) 

Thus S.Af.Av-, and hence its time -aver age h (0 o ), is nonnegative 
111 CL a 
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if the matrix with i-k entry 
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(3.51) 



is positive definite (or positive semidefinite). Equation 3.51 reduces 
to the previous criterion of positiveness, Eq. 3.44, when there is 
only one set of variables; it is the logical generalization fox devices 
with more variables. 

Assuming that the matrix of Eq. 3.51 is positive definite (or semi- 
definite), then frequency-power formulas of Type III are written 



800 



a 
^a 



(3.52) 



for all values of 8 & ; these are relationships of some value, because 
sums of powers are constrained to be positive. 

The requirement that the matrix of Eq. 3.51 be positive definite 
has been termed "local passivity" by Duinker [19, Ch. II]. Locally 
passive devices exhibit, at each operating point, small-signal equa- 
tions that define a passive system. An example of a device that is 
not locally passive at all points is the tunnel diode [22, 97], with a 
volt-ampere characteristic such as that of Fig. 3.1. The points of 




Fig. 3.1. Tunnel diode d-c volt-ampere curve. The 
negative -resistance region, is not a region 
of local passivity 

negative slope are not points of local passivity; the small-signal 
equivalent circuit at such points is active, not passive. 

Equations 3.52 are the basic frequency-power formulas of Type 
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HI; they are in a form similar to Page's relations [74] . Pantell 1 s 
form of the formulas are obtained in the limit as 9 & - : 

2 

S a p /^^ > (3.53) 

These are sometimes more easily applied to practical problems than 
Eqs. 3.52, but they are less powerful. 

It is interesting to note that even more general formulas can be 
written. If in the preceding proof we make a finite, but different, 
change in each independent phase cf> a , the finite change 9 a in each 
dependent phase 6 is 



a 
a 



and the result of Eq. 3.52 becomes 

F / 9co 

P 1 - cos -r-"" 9 * 
a a a. 9 a 



> (3.55) 



for all values of the parameters 9 a . Setting all 9 a equal to zero 
except one reduces this result to Eq. 3.52. In the limit as all 
9 - , but remain proportional to each other, 



2 



for all values of the parameters X a . Setting each X a in turn to 
one and all the others to zero, we obtain Eqs. 3.53. 

The formulas can be generalized even more. Following the ex- 
tension of Page [74], we let (x) De an y nonincreasing positive 
function defined from. to oo t and we define 



8w a \ P 
9^ 9 a) J 



Then 



for all values of the parameters 6 a . If (x) is the function 

fl < x < 1 

|(x) =1 (3.59) 

x > 1 

V. 

then Eq. 3.58 reduces to Eq. 3.55. The proof of Eq. 3.58 is iden- 
tical to the proof given by Page of a similar re suit [74] . 
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Frequency-power formulas of Type 131 apply to time-varying 
locally-passive dissipative systems, some examples of which are 
discussed in Chapter 5. In later chapters we usually quote the 
formulas as given by Pantell, Eqs. 3.53, with the understanding 
that the other forms, Eqs. 3.52, 3.55, 3.56, and 3.58, are also 
valid. 

3,5 Type IV 

Frequency-power formulas of Type IV relate reactive power in 
reactances. Like. the formulas of Type III, they are inequalities, 
instead of equalities. They are proved in a similar way and re- 
quire a similar condition of local passivity. 

When interpreting the formulas, one must keep in mind the 
points discussed immediately following Eq. 3.26. For inductive de- 
vices the formulas will be written (using the actual reactive powers) 
with a > sign; and for capacitive devices the formulas ought to hold 
with a <" sign. Multiport devices with both electric and magnetic 
energy storage do not obey the formulas, because of the ambiguity 
of sign, and likewise distributed systems do not obey them. Their 
main use is with single nonlinear time -varying elements, or with 
configurations of similar-type elements, with coupling (such as a 
rotating machine). 

When the formulas are quoted in future chapters, the form for in- 
ductive-type devices is given, unless otherwise stated. It is to be 
understood that if the device is actually capacitive, the sign should 
be reversed. In some places we give both forms for convenience. 

To prove the formulas, we assume the device is "reactive:" 
specifically, we assume it has constitutive relations of the form of 
Eq. 3.13, between the f. and x. variables. Using the represen- 
tations of Eqs. 3. 16 and 3. 17, we compute the finite change in the 
variables caused by a finite change 9 a in one of the independent 
phases cj> a . As in Sec. 3.4, this produces a finite change 



in each dependent phase 4> a , and so the change in each of the vari- 
ables is 

jw tj@ ,- / , * 

i . CL i CL i\ j3 Oil 

and 



Ax 



6 + Re ( - 1) 



"i - 3a)a v a 
If we multiply Af. by Ax and sum over i, and then take the time- 
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average of the sum, we obtain 

v< e a> - ^i^i^y 

= |Re S i Z a tc i ) a (f i ) a *(l - cos 9 a ) (3.63) 

or 

V(9J = S * (1 - C S 9 J (3 ' 64) 

a. a, CL (i) 

a 

where Q is given by Eq. 3.26. 

What we have so far is not of much value, because there is no 
estimate of the quantity k a f (6 a ) However, because the f^ vari- 
ables are functions of the x i variables, the changes Af^ can be 
evaluated by line -integrals in x.-space: 



(3 ' 65) 



Thus Z^Af^Ax^, and hence its time -aver age h a f (9 a ) , is nonnegative 
if the matrix with i-k entry 



8x k 



(3.66) 



is positive definite (or positive semi definite). This condition is one 
of local passivity, just like the condition for dissipative devices in 
Sec. 3.4. 

For locally -passive time-varying nonlinear reactances, then, fre- 
quency-power formulas of Type IV hold in the form 

Q a / 8co a 

' C S 



for all 8 a , provided care is used in interpreting the Q^ quantities, 
as discussed earlier. 

Just as in Sec. 3.4, we can write other forms of these formulas. 
In the limit as 6 a - , we find 



800 / - - (3 ' 68) 

u, a 

These formulas are easily applied to practical problems, although 
they are not as powerful as Eqs. 3.67. If in the preceding proof we 
make a finite change in each of the independent frequencies, we ob- 
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tain frequency -power formulas similar to Eqs. 3.55: 

Qa[ / 8 "a \1 

S 1 - cos S 9 >0 (3.69) 

a co V a 9co a / x /; 

a L a J 

for all values of the parameters 9 . Setting all 9 to zero except 
one reduces Eq. 3.69 to the simpler form of Eq. 3.67. If all 8 a ap- 
proach zero, but remain proportional to each other, we find 



< 3 ' 70 > 



for all values of the parameters X . Setting each x * n ^ urn 
one and all others to zero, we obtain Eq. 3.68. And finally, if 
g a (9 ) is defined, as before, by Eq. 3,57, then 



for all values of the parameters a . 

In later chapters we usually quote the formulas of Type IV as given 
by Eq, 3.68, with the understanding that the other forms, Eqs. 3.67, 
3.69, 3.70, and 3.71, are also valid. 

There are devices that are not locally passive, for which we would 
like to have frequency-power formulas of Type IV. One such device 
is the induction motor, in which the electrical variables for each 
fixed shaft position form a locally passive system. When the mech- 
anical variable is added, however, there can be regions of unstable 
equilibrium, and regions in which the electromechanical system is 
not locally passive. Can we find frequency -power formulas of Type 
IV for similar cases? 

In some cases we can. If a small matrix within the matrix of Eq. 
3.66 is positive definite, let us call the remaining variables (those 
that are not related by this submatrix) the ''offending" variables. 
Thus, for the induction motor, the mechanical variables would be 
the offending ones. It may happen that none of the offending vari- 
ables depends in its actual time variation on one or more of the co a 
independent variables. In this case, the corresponding h a ! is 
nonnegative, since the summation in its definition of Eq. 3.63 can 
be restricted to the non-offending variables. Thus f with this "weak 
positive-definite" system, there are still some frequency-power 
formulas applicable. 



PART II 



SPECIFIC SYSTEMS 



Chapte r 4 
LUMPED REACTIVE SYSTEMS 



In this chapter we discuss some lumped devices that obey fre- 
quency-power formulas of Types I and IV. Generally speaking, these 
are reactive systems, although in Sec. 4.11 it is observed that linear 
resistors obey formulas of Type IV. Our purpose is to establish a 
partial list of devices that obey the formulas. 

In Sec. 4.1 we discuss the energy functions that are suitable for 
proving the Manley-Rowe formulas, and show that one may use the 
internal energy (expressed in terms of the proper variables), or 
any energy function related to it by a Legendre transformation. 

The remaining sections deal with specific devices. Sections 4.2, 
4.3, and 4.5 discuss nonlinear energy storage elements capaci- 
tors, inductors, and mechanical elements. In Sec. 4.4, we show 
that a time-varying capacitor, although not lossless, obeys the 
Manley-Rowe formulas. 

Many rotating machines (those without commutators) obey the 
formulas, and in Sec. 4.6 a general model for such a machine is 
discussed. A logical generalization of this model is discussed in 
Sec. 4.7; it includes both magnetic and electric coupling to mech- 
anical circuits, and obeys the Manley-Rowe formulas. 

There are very few idealized nonlinear circuit elements. The 
ideal diode is one; others, called "traditors, " have been proposed 
by Duinker [20] , In Sec. 4.8 we show that these, unlike the ideal 
diode, obey the Manley-Rowe formulas. 

Some degenerate elements that satisfy frequency-power formulas 
of Types I or IV are discussed in Sec. 4.9. These devices gener- 
ally satisfy the formulas in a trivial way, with all P a or Q a (or 
both) equal to zero. The major purpose of discussing them at all 
is to aid the analysis of entire networks. It is shown in Sec. 4.10 
that a network of elements, each of which obeys frequency-power 
formulas of Type I, itself obeys the formulas. Thus, in particular, 
these "degenerate" elements of Sec. 4.9 can appear in such net- 
works. Included are the ideal transformer, the gyrator, the dif- 
ferential gear, and (in a formal way) the maser. 

Many distributed systems obey the formulas of Type I; these are 
discussed in Chapter 6. 

4.1 Energy Functions 

It is shown in Sec. 3.2 that any device with an energy state func- 
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tion obeys the Manley-Rowe formulas. The function used there is 
(at least in the time -invariant case) the energy of the device, ex- 
pressed in terms of the x. variables, and the expression for P a 



is easily interpreted as the power input to the device at frequency 

a The same method, however, applies to any energy function, as a 
function of any of the variables. The formulas are of the same form, 
but the expression for P is different from that of Eq. 4.1. Since 
a given physical system can be described by many different energy 
functions, it obeys many different Manley-Rowe formulas. If we 
want to interpret the P a as power inputs at frequency o> a , only a 
few energy functions are appropriate. What combinations of energy 
functions and variables are suitable? 

It is known that the internal energy U, as a function of the x i 
variables, is suitable; we show here that any function related to 
U(x- , t) by a Legendre transformation [30, Sec. 7.1] can be used. 

Because 

f. = |2 (4.2) 

the most general form for any energy function related to U(x^ , t) 
by a Legendre transformation is 

IPfr., f v , t) = U(x, , t) - S,x v f v (4.3) 

y K. i JK > -K- 

or its negative, where i is an index that runs over all the variables, 
k is an index that runs over some of the variables, and i is an in- 
dex that runs over the rest. The Legendre transformation of Eq. 
4.3 defines explicitly the variables on which the new function depends, 
and 

121 s f (4.4) 

and 

IF - -* ^ 

8f k * 

To have U 1 satisfy the condition that it be bounded, we must require 
that all v k() vanish, as is evident from Eqs. 3.17 and 4.3. 

The procedure of Sec, 3.2, when applied to U T (x. , f, , t) yields 
Manley-Rowe formulas in the form of Eq. 3.36 with 

P, a (4.6) 

ko v ' 
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and 



P ss v f 
Jfo o o (4.7) 



P a = i 



= iR 




These are exactly the forms obtained by using Ufe^ , t), so any in- 
terpretation given the P a in that case applies here also. 

It is interesting that a given physical system obeys Manley-Rowe 
formulas with many different forms for P . One can even show 
that for any (reasonable) function of the voltage w(e) there is a 
set of Manley-Rowe formulas with 



that hold for the nonlinear capacitor. By thinking up different energy 
functions and different variables for them to depend on, we come 
upon many other forms. Because our desire is to interpret the P 
as power inputs at frequency w , we will not pursue these other 
forms for P a further. We find it sufficient to use U(x^ ? t) or a 
function related to it by a Legendre transformation. 

In a similar way, a system that obeys frequency-power formulas 
of any type also obeys similar formulas with different P a or Q Q . 
Since our purpose is to interpret the P and Q as real and re- 
active power at frequency co a , we will not have occasion to discuss 
these other forms. 

4.2 Nonlinear Capacitor 

This device is probably of most practical importance; it is the 
one originally discussed by Manley and Rowe [60] . A nonlinear ca- 
pacitor without losses and without dispersion is characterized by a 
single-valued charge-voltage relation 

e = e(q) (4.10) 

The internal energy W c (q) is 

r q 

W c (q) = e{q) dq (4.11) 
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so that 



dW c (q) 
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(4.12) 



and the coenergy [14, 62] is defined by the Legendre transformation 

I) (4.13) 



W c '{e) = eq- 



These energy functions can be given a graphical interpretation, as 
shown on the voltage-charge characteristic in Fig. 4.1. 




Fig. 4,1. A typical charge -voltage characteristic of a 
nonlinear capacitor. At the operating point 
indicated, the energy and coenergy are rep- 
resented by the areas A and B respectively 

Either the energy or the coenergy can be used to obtain frequency- 
power formulas of the form (there is no d-c power) 



P a 9u 
2 

a " 9co 



with P given by 



(4.14) 



(4.15) 



or simply the power into the capacitor at frequency w . 

The nonlinear capacitor also obeys the various forms of frequency- 
power formulas of Type IV discussed in Sec. 3.5, such as 
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where 

Q a = I Im e Q i a * (4.17) 

is the reactive power into the capacitor at frequency oo . 

Note that a linear capacitor satisfies these formulas In a trivial 
way, with all P^ = and all Q a < . The frequency-power form- 
ulas lose their importance, since much stronger conservation prin- 
ciples hold. Nevertheless, the frequency-power formulas are still 
valid, so the inclusion of linear capacitors in a network does not pre- 
vent it from obeying the formulas . 

4.2 Nonlinear Inductor 

The nonlinear inductor is treated similarly. It is characterized 
by a relation between the current and the flux linkage of the form 

i * i(X) (4.18) 

and so the stored magnetic energy is 



X (4.19) 

*j 
o 

with 

i - %^ (4.20} 

The use of W^X) as an energy state function leads to frequency- 
power formulas of the form of Eq. 4.14 with P n given by Eq. 4.15. 
The magnetic coenergy W *(i) 






can also be used, with the same result. 

The constitutive relation of Eq. 4.18 leads to frequency-power 
formulas of Type IV of the form 

\ 2 Q 

-a > (4.22) 

w a 

or of any of the other forms derived in Sec, 3.5, with 

equal to the reactive power input to the inductor. Notice that in 
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the frequency-power formulas of Eq. 4.22 the sign is different from 
the corresponding formula for the nonlinear capacitor, Eq. 4.16. 

Linear inductors satisfy Eqs. 4.14 and 4.22 in a trivial way, with 
all P = and Q > . The frequency-power formulas lose their 
importance, but not their validity. 

Two or more coils coupled nonlinearly also obey the formulas. 
Suppose, for two coils, the constitutive relations are written 



and 



Then formulas of Type IV hold in the form of Eq. 4.22, with 

(4.26) 



being the reactive power at frequency u> a flowing into both coils. 
Also, if 



(4.27) 



9i 2 



as is necessary for the device to be lossless, thea the magnetic co- 
energy 

< 4 - 28 > 

is well defined [42, Sec. 6.10; 80] . This function leads to frequency- 
power formulas of Type I, in the form of Eq. 4.14, with 

(4 ' 29) 



interpreted as the power input to both coils at frequency oj a . The 
extension to more than two coils is obvious. 

4.4 Time -Varying Capacitor 

This device probably causes more confusion than any other we 
discuss. The confusion centers around the meaning of the terms, 
"time -vary ing" and "lossless. 11 

A given physical device may be accurately represented by any of 
several models, depending on circumstances. The time-varying 
capacitor is a model (not itself a device) that accurately describes 
many physical systems with some excitation previously assigned,. 
Many people consider a variable capacitor (such as a radio tuning 
coadenser) to be time-varying; our point of view is that if the mech- 
anical motion is specified in advance, then a valid model, as far as 
the electrical variables are concerned, for this device is a time- 
varying capacitor. If the capacitance is linear, the equation of 
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motion is 

q(e, t) = C(t)e (4.30) 

where C(t) is the capacitance. 

We obviously cannot predict, within the framework of the time- 
varying model, anything about the physical mechanism that provides 
the time variation. If the mechanical motion is not specified in ad- 
vance, we cannot use this model; instead we must use a nonlinear 
time -invariant model that includes mechanical as well as electrical 
variables. 

Another example might be helpful. If a nonlinear capacitor is 
pumped strongly by current at frequency co , and if there are super- 
imposed on this small currents at frequencies w g + no , where n 
is a positive or negative integer, we may wish to restrict our at- 
tention to the small- signal variables, without the pump. We may do 
so by using a linear time-varying model, with capacitance C(t) just 
equal to the incremental capacitance of the nonlinear capacitor, as 
pumped. Such a technique is universally used in descriptions of 
frequency converters, because the linear time-varying equations are 
easier to solve than the original nonlinear ones. This "lineariza- 
tion" process is described in more detail in Sec. 8.2. 

A nonlinear capacitor is not lossless. By "lossless" we mean 
that the time-average power input vanishes: 

S a P a = (4.31) 

The time-varying capacitor does not obey Eq. 4.31, so it cannot be 
said to be lossless. Net power can be generated or destroyed by 
this model. Of course, the actual physical device may be lossless, 
in which case this energy is accounted for by loss or gain of energy 
by the "prime mover" or "pump, " but this fact is not predicted by 
the time -vary ing model. 

Simply because the time-varying capacitor is not lossless, we do 
not wish to imply that it is lossy, that is, having dissipative loss. 
One connotation often given to the word "lossless" is freedom from 
dissipation, with its noise, entropy increase, etc. This is not ex- 
actly what we mean here. 

It is stated in Chapter 1 that losslessness is neither necessary nor 
sufficient for the Manley-Rowe formulas to hold. The example of 
Appendix A shows that it is not sufficient; the example here shows 
that it is not necessary: the time-varying capacitor, although not 
lossless, obeys the frequency -power formulas. 

Now we show that the nonlinear, time-varying capacitor obeys the 
formulas of Types I and IV (although, of course, only those formu- 
las corresponding to an independent frequency w , on which the 
explicit time variation does not depend). The nonlinear, time -varying 
constitutive relation is 
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e = e(q, t) (4.32) 

and the reasoning in Sec. 3.5 leads to frequency-power formulas of 
Type IV, in the form of Eq. 4.16, with the Q a given by Eq.4.17. To 
show that the device obeys the Manley-Rowe formulas, we define the 
time -varying energy function 

r q 

W c (q, t) = ^ e(q, t) dq ( 4 - 33 ) 



This energy function leads to formulas of the form of Eq. 4.14, with 
the P Q given, as expected, by Eq. 4.15. 

Similar results hold for the time-varying inductor. 

4. 5 Mechanical Elements 

Frequency-power formulas of Types I and IV hold for nonlinear, 
time-varying mechanical energy- storage elements. We do not inves- 
tigate Type IV, since reactive power is not often used with mechani- 
cal elements. 

A nonlinear, possibly time-varying spring has a constitutive re- 
lation 

f = f(x, t) (4.34) 

where f is the force exerted on the spring and x is its displace- 
ment. Although we use one -dimensional motion in the example, the 
notions are readily extended to more than one dimension, and to ro- 
tational systems. A simple example of a nonlinear (but time-invari- 
ant) constitutive relation like Eq. 4.34 is that between torque and 
angular displacement of a pendulum, 

T = mgr sin6 (4.35) 

The potential energy 

W p (x, t) = Cf(x, t)dx (4.36) 

is an energy- state function, and its use leads to frequency-power 
formulas in the form of Eq. 4.14, with P Q interpreted as the mech- 
anical power input at frequency u> a : 

P = iRe f*v (4.37) 



Q 



where v is the velocity, the time derivative of the displacement. 

In a similar way, one can prove the Manley-Rowe formulas for a 
nonlinear, time-varying point mass with a constitutive relation 

v = v(p, t) (4.38) 

where p is the momentum of the particle, or the time integral of 
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the force. The use of the kinetic-energy function 

W k (p, t) = J v(p, t)dp (4.39) 

leads to frequency-power formulas of the form of Eq. 4.14, with the 
P a given by Eq. 4.37. 

It might be argued that in reality the constitutive relation of Eq. 
4.38 for point masses is linear, in the form 

p = mv (4.40) 

where m is the mass, and so the Manley-Rowe formulas hold in a 
trivial way, with all P = . This is true for nonrelativistic parti- 
cles, but for relativistic speeds Eq. 4.40 is modified to be nonlinear: 



4.6 Rotating Machines 

Some rotating machines obey the Manley-Rowe formulas. These 
machines are those without commutators; that is, they are suitable 
for use as induction, synchronous, or reluctance machines. 

As a model for these machines, we consider a set of coils - 
some of them fixed (on the stator) and others movable (on the rotor). 
We neglect losses like those due to winding resistance, eddy cur- 
rents, mechanical friction, and hysteresis. The rotor and stator 
shape need not be specified, so salient-pole machines are included. 
Neither the number of windings nor their configuration is specified, 
so we can have any number of phases and any number of poles. The 
magnetic material can even be nonlinear (although it must be free 
from hysteresis and dispersion). This model is a bit more general 
than the "generalized machine" [116, Ch. 3], which has been proposed 
as a model for practical magnetic-field rotating devices without 
commutator s . 

The coil currents i. (k is an index) and shaft position can be 
considered as independent variables to define the magnetic coenergy 
[116, Ch. 3] by the line integral 

k di k (4 42) 

so that 

8W '(i , 40 

Xk = _^A (4.43) 

is the flux linkage in the k-th coil, and 

3W '(i , <|>) 
T = ***. (4.44) 
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is the torque applied to the shaft [116, Ch. 3] . 

The shaft speed co = "ff can have an average value, so if cj> is 
written in the form of*Eq. 3.17, v- need not vanish. There can be 



10 



d-c power input at the shaft, although not at any of the coil termin- 
als. It is primarily with this application in mind that provision is 
made for d-c power terms in the Manley-Rowe formulas. 

The magnetic coenergy is a suitable state function for deriving 
frequency-power formulas of Type I; they are of the form 

Jl 5SL- + S !L i o (4.45) 



where w is the average value of oo m , P Q is the average torque 
times CD, and each P a has both electrical and mechanical contri- 
butions: 

(4.46) 



It is evident from Eq. 4.45 that no error would be made if the d-c 
power input at the shaft were considered as being at frequency co . 
This is true in general for rotating shaft d-c power, as proved in 
Appendix D. 

These rotating machines also obey frequency-power formulas of 
Type IV, provided the frequency u> a does not appear in the actual 
time dependence of the shaft angle <|> . The overall derivative ma- 
trix of the form of Eq. 3.66 is not positive definite, because the 
diagonal term relating the mechanical variables, 

Mi., <tt 

(4 ' 4?) 



is not necessarily positive. The electric submatrix, however, with 
i-k entry 



is ordinarily positive definite, so the rotating machines satisfy the 
"weak positive-definiteness" condition discussed in Sec. 3.5, and 
therefore obey formulas of Type IV for each co that does not ap- 
pear in the actual time dependence of <|>(t) . The formulas are writ- 
ten in the form of Eq, 4.22, with Q the overall reactive power in- 
put at frequency o> , to all coils, 

i^* (4.49) 

These formulas, of Type IV, can be used [86, 85] to show that the 
induction generator is stable without some source of reactive power. 
Rotating machines with commutators do not ordinarily obey these 
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frequency-power formulas, because the commutators are switches, 
or time-varying resistors. However, in some cases where it is 
possible to treat the machine as if these switches were "outside", 
the formulas hold. This is done, for example, in Sec. 7.3. 

4.7 Generalized Conservative Energy Conversion Device 

A logical generalization of the rotating machine is the generalized 
conservative energy conversion device [116, Sec. 1.4] , which is a model 
that describes many physical systems with either magnetic field or 
electric field energy conversion, or both. This model applies to 
such physical devices as rotating machines, phonograph cartridges, 
microphones (except carbon microphones), loudspeakers, solenoids, 
meter movements, etc. The general model is described by a La- 
grangian [116, Sec. 1.4] , and this state function can be used to obtain 
the Manley-Rowe formulas. The expressions for P a must include 
mechanical as well as electrical variables. 

4.8 Traditors 

There are not many idealized nonlinear circuit elements. Duinker 
[20] has found a set of elements, "traditors, " that can be described 
by a Lagrangian, but do not store energy, so cannot be described by 
a Hamiltonian. He defines several classes of traditors, but the most 
important seem to be the "third-degree" traditors, the simplest to 
be actually nonlinear. He shows how higher-degree and lower-degree 
traditors can be synthesized from third- degree traditors. We show 
here that traditors of the third degree (and hence of any degree) 
obey frequency-power formulas of Type I, a result also obtained 
by Duinker [18]. 

Six distinct forms of third-degree traditors are discussed by 
Duinker [20], but any can be constructed from any other one and 
gyrators. Thus, it is necessary to discuss only one form here 
we take "form I. " 

A third-degree traditor of form I is a three-port device with 
equations of motion 

GI . -Aq 2 i 3 (4.50) 



where e,, i k , and q fe are respectively the voltage, current, and 
charge at the k-th port. The constant A is the only device param- 
eter. This device has an instantaneous power input of zero, just like 
the gyrator and ideal transformer. However, unlike the ideal trans- 
former and the gyrator, it is nonlinear, so the power inputs at the 
various frequencies need not vanish. 
The device has a Lagrangian 
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L = A qi q 2 i 3 (4-53) 

which is a function of ^ , q 2> and i 3 . Its partial derivatives are 

H = - e i (4-54) 

*- = -e (4.55) 

8q 2 * 

and 

r* 

dt (4.56) 



and so the use of -L as the energy state function gives frequency- 
power formulas in the form of Eq. 4.14, with 



Fa"** 



= 4Re 



[(e 1 )*(i 1 ) a + (e 2 )*(i 2 ) a + <e s )*(i,) a ] (4.57) 



simply the power input to all ports at frequency co a . Although trad- 
itors connected together will pass d-c power, this cannot happen if 
all the variables are written in the forms of Eqs. 3.15, 3.16, and 
3.17. We need not, therefore, bother with the d-c power terms of 
the Manley-Rowe formulas. 

Traditors do not, generally, obey frequency-power formulas of 
Type IV, because they are not locally passive. 

4.9 Degenerate Elements 

In this section we discuss some elements that obey the frequency- 
power formulas in a trivial way. First we discuss the ideal trans- 
former, and then the gyrator. Next we give a generalization, and 
then apply it to a gear system, such as a differential. The results 
formally apply to the maser, as we point out at the end, but they do 
not give much physical insight. 

Ideal Transformer. By the term, "ideal transformer, " we mean 
a two-port device with equations 

e 2 (t) = Ne/t) (4.58) 

and 

ij(t) = -Ni 2 {t) (4.59) 

where the constant N is called the n turns ratio. M Except for triv- 
ial cases (N = , _+ 1 , or oo), the model is approximated physically 
by a pair of coils on a common magnetic core. With perfect coupling 
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between the coils, the flux linking one also links the other, so the two 
flux linkages (and therefore the voltages) differ only by a constant 
(the turns ratio), so Eq. 4.58 holds. The constitutive relation of the 
core material gives a relation between (i 1 + Ni 2 ) and the flux link- 
age X 2 , and without hysteresis or other losses this might be approx- 
imated by a curve like Fig. 4.2. If the excitation is such that | \ 2 \ 



en 
O 
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Fig. 4.2 Idealized B-H curve for a material suitable 
for a transformer. For values of \ 2 less 
than the critical amount, the magnetic field 
is vani shingly small. To keep the flux below 
this critical value, we must avoid having a 
non-zero average value of voltage; hence the 
physical transformer (using this material) 
cannot pass d-c power 

is always less than the critical value shown in Fig. 4.2, then ^ + Ni 2 
vanishes, so Eq. 4.59 holds. If, however, \ z \ gets too large, then 
Eq. 4.59 does not hold, and the device is not an ideal transformer, be- 
cause it has saturated. The physical transformer cannot pass d-c 
power, but the ideal transformer (which is only a model) can. 

The physical transformer is merely a pair of nonlinearly-coupled 
coils, and from the results of Sec. 4.3 we know it obeys frequency- 
power formulas of Types I and IV. We now show that the ideal trans- 
former model also satisfies these formulas, but in a trivial way. 

The vector power into the ideal transformer at frequency o> a is 
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and because Eqs. 4.58 and 4.59 hold in frequency -by-frequency form, 
Eq. 4.60 must vanish. Thus all Q vanish, so any sum of them van- 
ishes, and so in particular frequency-power formulas of Type IV, with 
the inequality sign in either direction, hold. 

By a similar argument, one can easily see that the d-c power sum 
of Eq. 3.36 vanishes, and since all P a vanish, Eqs. 3.36, the Manley- 
Rowe formulas, hold. 

Gyrator. The gyrator [104, 94] is a nonreciprocal circuit model 
with two ports, and equations of motion 

e z (t) = Rij(t) (4.61) 

and 

Cl (t) = - Ri 2 (t) (4.62) 

Like the ideal transformer, the gyrator does not store or dissipate 
any energy, but unlike the transformer it can accept arbitrary amounts 
of reactive power, of either sign. It obeys the Manley-Rowe formulas, 
but does not obey frequency -power formulas of Type IV. 
The vector power at frequency w a is 

P a + J Q a= 



and therefore is pure imaginary. Thus all P^ vanish, and because 
one can show easily that the d-c part of Eq. 3.36 vanishes, the fre- 
quency-power formulas of Type I, Eq. 3.36, hold, with the ? a in- 
terpreted as the net power input to the gyrator at frequency co a . 

A Generalization. These results can be generalized somewhat. 
Suppose a device cannot store or dissipate energy so the instantan- 
eous power input 

P(t) = S.fWv-tt) = (4.64) 



ables , 

K, 



where the f. variables do not depend at all on the v. variables; 

and suppose that the v^ variables are not all independent, so that 

they can be written in terms of a smaller number of independent vari- 

(4.65) 



For simplicity let us use matrix notation; let v(t) be the column 
matrix of V|(t) ; let v f (t) be the column matrix of vHt) ; let (t) 
be the column matrix of the f^(t) ; and let C be the rectangular 
matrix with entries C., . Then, 

v = Cv 1 (4.66) 

and 
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= f . v 

>/<* 

= ^Cv* (4.67) 

where the t subscript indicates the transpose. Since the v* vari- 
ables are independent, we must have 

f.C = (4.68) 

^w /* 

If we express all variables in the form of Eqs. 3.15 and 3.16, we 
may derive from Eqs. 4.66 and 4.68 the equations relating the aver- 
age values of the variables, 



and the equations relating each set of Fourier coefficients, 

v = Cv (4.71) 



and 



la 1 " = (4 - ?2) 



where the dagger | indicates the complex conjugate transpose of 
a matrix. 

The d-c terms in the Manley-Rowe formulas of Eq. 3.36 are writ- 
ten, in our notation, 

z !k !!!i- f ^L. 

,4.73, 



a, 

= 

and the vector power at each frequency o> a is 



(4J4) 



and so the Manley-Rowe formulas, Eq. 3.36, hold in a trivial way, 
and the frequency -power formulas of Type IV hold, also in a trivial 
way, provided Q can be interpreted as reactive power. 
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This generalization clearly applies to both the ideal transformer 
and the gyrator. In the former case the Q^ can be interpreted as 
reactive powers, and in the latter case they cannot. We now show 
how this generalization applies to gear systems, and then, formally, 
to masers. 

Gear Systems. If we neglect energy storage and friction in gear 
systems, the preceeding generalization applies to them. A typical 
gear system of interest is the differential; we use it as an example. 
The shaft speeds are identified with the v. variables, and the shaft 
torques T with the f- variables. The speed constraint for the dif- 
ferential is in the form 

w, = r(co -co 2 ) (4.75) 



Letting co 2 and GO be independent variables, Eq. 4.66 becomes 

"- r r " 



1 
1 



X 



w, 



(4.76) 



and the condition on torques, Eq. 4.68, is 



-r r 
1 
1 



= [0 0] 



(4.77) 



If only d-c power is present, then the Manley-Rowe formulas re- 
duce to the form of Eq. 1.16, which is well known. 

It is interesting to note that in Eq. 1.16 all the power flow is d-c, 
yet the form of the equations is what one would expect from a-c 
power terms. This suggests that one may actually consider d-c 
power flow terms at a rotating shaft to be power at frequency equal 
to the shaft speed. This is generally possible, as we show in Ap- 
pendix D, 

Maser. The solid-state maser formally is described by the gen- 
eralization of Sec. 4.9, although this does not provide much physical 
insight. The results are obtained more easily by the method of 
Weiss[115], 

In practice, the "idler" power given off by a maser is absorbed 
by losses in the crystal, as is a majority of the pump power. If, how- 
ever, these losses are considered as being "outside" the system, then 
frequency-power formulas of Type I apply. 

Let us identify the energy level separations Ej, with the v. vari- 
ables, and the number of photons per second entering the system at 
frequency v-. = E,,/h, where h is Plank's constant, with the f. 
variables. The energy level separations are not all independent!" for 
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certain of them add together to equal certain other ones. In the sim- 
ple three-level maser, the two smaller energy level separations 
added together just equal the large one. An independent set of E^ 
may be extracted, or, what is equivalent, an independent set of v^. 
There is a constraint among the number of photons entering the sys- 
tem at each frequency; in the three-level maser the number at each 
frequency must be equal (in magnitude). The formal application of 
the generalization of Sec. 4.9 leads to frequency -power formulas 
similar to Eq. 1.13, which holds for the three-level maser. Unfor- 
tunately, this treatment does not provide much physical insight; it 
is merely by mathematical analogy. 

4.10 Network Extension 



Suppose we have a network of elements like nonlinear, time -vary- 
ing capacitors and inductors, rotating machines, ideal transformers, 
and other devices that obey the Manley-Rowe formulas, and suppose 
this network has a number of terminal pairs, through which power 
flows to the network. If the entire network is to be used for frequency 
conversion, the quantities of engineering interest are the powers at 
each frequency at the ports of the network. Even though each ele- 
ment in the network obeys the Manley-Rowe formulas, we do not as 
yet know that the network as a whole does, in such a way that the 
powers P can be evaluated at the ports. This is what is proved 
in this section. 

A similar extension of frequency-power formulas of Type IV is 
possible if the ambiguity of sign of the reactive power can be ac- 
counted for. The extension holds if the devices are either all capa- 
citive or all inductive, but it is not valid if the network contains both 
capacitive and inductive elements. 

These extensions are intuitively reasonable, and are proven here 
by general topological arguments. Let the terminal pairs of the in- 
dividual elements in the network be indexed by k . Then the overall 
vector power into all the elements at frequency w a is 



If each element obeys frequency-power formulas of Type I, then by 
adding together the formulas one arrives at formulas of the same 
form with the P of Eq. 4.78, and alternatively, if each element 
satisfies formulas of Type IV (all with the same direction of inequal- 
ity), then by adding together formulas one arrives at formulas of the 
same form with the Q a of Eq. 4.78. What we show here is that the 
P + jQ are identical to the vector powers entering the network at 
its port a s. Therefore, the network obeys the appropriate type of fre- 
quency-power formulas, with the powers (P a or QJ evaluated at 
the ports of the network. 

To show that P + jQ a is the vector power into the network at 
frequency u> , we%ise a%opological argument [33, Ch.l] . We rcpre- 
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sent multiport devices within the network by a scheme proposed by 
Guillemin[34], and for topological purposes we consider each port 
of the network to be a separate branch, with voltage e . and current 
-i. , where i is an index over the ports of the network. The branches 
of the network consist, then, of the elements and the terminal pairs. 
We call the column matrix of branch currents i, and the column 
matrix of branch voltages e . Then one of Kirchhoff ! s laws is of the 
form 



where, picking any tree [33, Ch. 1], we denote the column matrix of 
link currents by i 1 . The matrix T is the tie-set schedule [33, Ch. 1] . 
The other one of Kirchhoff ! s laws is 

Te = (4.80) 

<+* < 

Taking Fourier coefficients of Eqs. 4.79 and 4.80, we find 

it =iJtT (4.81) 

i (I f^Q, r+f 

and 

Te = (4.82) 

"""a 

We wish to show that 
This is clear, since 



= (4.84) 

Similar reasoning applies to all-mechanical and electric -mechani- 
cal networks. 

4.11 Linear Time-Invariant Resistors 

Linear time -invariant resistors are not reactive devices, but 
they are included here because they obey frequency -power formu- 
las of Type IV. No reactive power flows in at any frequency, so all 
Q a vanish. The formulas then are satisfied in a trivial way, and 
linear resistors can be included in networks that are intended to obey 
formulas of Type IV. 



Chapte r 5 
LUMPED DISSIPATIVE SYSTEMS 



In this chapter we shall discuss some systems that satisfy fre- 
quency-power formulas of Types II and III. Generally speaking, such 
systems are dissipative, although the degenerate cases discussed in 
Sees. 5.3 and 5.5 are riot. 

In Sec. 5.1 nonlinear resistors are treated, and in Sec. 5.2 time- 
varying resistors. Some degenerate devices appear in Sec. 5.3, and 
in Sec. 5.4 it is shown that a network, each element of which obeys 
frequency-power formulas of one type, itself obeys formulas of the 
same type. In Sec. 5.5 it is pointed out that linear lossless devices 
obey formulas of Type HI. 

5.1 Nonlinear Resistors 



A nonlinear resistor is characterized by a dissipation function, 
the content [68]. The voltage e is some function of the current i, 

e = e(i) (5.1) 

and so the content, defined by 



G(i) = ^ e(i) di (5.2) 

is a function of the current. In addition, one can define another 
state function, the co-content [68] , by the Legendre transforma- 
tion 

G l (e) as ie - G(i) (5.3) 

A typical volt-ampere curve for a resistor is shown in Fig. 5.1; the 

state functions correspond to the areas indicated. 

In demonstrating that the nonlinear resistor obeys frequency-power 
formulas of the second type, we may use either the content or the co- 
content. Let us use the content. Then the conditions set up in Sec. 
3,3 are satisfied, and we find 

; 3u> 

= s *r - (e *i - i *e ) 
a 4 9co a a a a' 
a 
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Fig. 5.1. A typical voltage -cur rent characteristic of a 
nonlinear resistor. At the operating point in- 
dicated, the content and co-content are rep- 
resented by the areas A and B respectively 



where 



Q a = *i 



e i * 
a a 



(5.5) 



In addition, we know from Sec. 3.4 that because of the functional 
relation of Eq. 5.1, the device satisfies frequency-power formulas 
of the third type, in any of the forms derived in Sec. 3.4, such as 
PantelPs form 



V __ ^C-> 

~ 



where 



e i * 
a a 



provided that the incremental resistance 



(5.6) 



(5.7) 



(5.8) 



is nowhere negative. 

A nonlinear resistor of great theoretical importance is the ideal 
diode, whose volt-ampere curve is shown in Fig. 5.2. This device 
is lossless (because at all times either the voltage, the current, or 
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Fig. 5.2. Volt-ampere curve of an ideal diode. Although 
the device is lossless, its operation is described 
by a relation between voltage and current, so 
it is a nonlinear resistor 

both vanish), but it is properly classed as a nonlinear resistor, be- 
cause its operation is defined in the voltage-current plane. It satis- 
fies frequency-power formulas of both Types II and III. 

As far as frequency conversion is concerned, the most important 
physical devices that satisfy these formulas are point-contact diodes, 
or crystal diodes [106, Ch. 5; 108; 66; 107], which are in wide use in 
communications systems as mixers, modulators, etc. Some other 
physical devices that are well approximated as nonlinear resistors 
are rectifiers, varistors, field-emission diodes [26] , and special- 
purpose semiconductor diodes, including very narrow base diodes 
[84] , avalanche (Zener) diodes [63] , and current limiters [ill] . All 
these obey formulas of both Types II and III. 

A few nonlinear resistors have a "negative-resistance" region, 
and therefore cannot be expected to obey formulas of Type III, al- 
though they do obey formulas of Type II. Typical are the tunnel 
diode [22, 97] (if the barrier capacitance is neglected), the p-n-p-n 
switching diode [69, 95, 96] , etc. [81] . 

Some nonreciprocal two-port nonlinear resistors obey formulas 
of Type III, but not Type II, since the dissipation fuction is not de- 
fined. Included are the Hall effect gyrator [32, 61], and other chunks 
of ohmic material, with two or more ports, immersed in a magnetic 
field. 

Nonlinear friction devices, including dashpots and ratchets, also 
obey the formulas, as one can see by analogy. 

5. 2 Time-Varying Resistors 

Time variation of a resistance can arise in many ways. Probably 
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the simplest time-variable resistor is one that is varied mechani- 
cally a potentiometer, and, as a limiting case, a switch. Since, 
however, energy is not required to vary a resistor (contrary to the 
case of a capacitor or an inductor), many physical effects can be 
used. 

Resistors that are varied mechanically include potentiometers, 
rheostats, switches, and strain gages. Carbon microphones and 
piezoresistive microphones [9] are nothing more than resistors 
whose instantaneous values depends on the sound pressure. Therm- 
istors (and in fact, all resistors) have a resistance that depends on 
temperature, and some photoconductors can be thought of as linear 
resistors whose resistance depends on the light intensity. The re- 
sistance matrix of a chunk of ohmic material immersed in a mag- 
netic field depends on the strength and orientation of the field. As 
a last example, the magnetic field dependence of superconductivity 
has been used to make time-varying resistors (cryotrons) [8], and 
these have been considered for possible frequency conversion ap- 
plications [118] . 

Just as it is not necessary in Sec. 4.4 to think of a time-varying 
reactance in terms of the physical effect that produces the varia- 
tion, so here it is not necessary to determine whether the time vari- 
ation is caused by any particular physical effect, or by heavy non- 
linear pumping (see Sec. 8.2). 

A nonlinear time-varying resistor has a voltage that is a function 
of current and also explicitly of time 

e = e(i, t) (5.9) 

and so one can define a time-varying dissipation function 

G(i, t) = fe(i, t)di (5.10) 

and the proof of the second type of frequency-power formulas, in 
Sec. 3.3, predicts that 

800 



where 



Q a = ^ Va* (5J2) 



for each independent frequency ca a that is not used to help specify 
the explicit time variation of G. 

Similarly, frequency-power formulas of the third type hold in any 
of the forms discussed in Sec. 3.4, including the form 

/ 8 \ 2 

h = S 1 - JP > (5.13) 

a- a V 8ca / a - \ i 

a. 
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where 

Pa = * Re Va* (5 ' 14) 

provided the excitation is such that at no time is the time-varying 
incremental resistance 



R . (t) , (5.15) 

negative. Equation 5.13 holds for each independent frequency o) a 
that is not used to help specify the explicit time dependence of e(i, t) . 

By analogy, time-varying dashpots and ratchets and periodically- 
applied mechanical clamps obey the formulas of Types II and III, 
just as their electrical counterparts do. 

It is interesting to note that a linear time -invariant resistor will 
accept no reactive power at any frequency. A nonlinear resistor can 
accept net reactive power, and is capable of accepting it at some 
frequencies and releasing it at others, but there is a sum of reactive 
powers that vanishes: 

2 Q = (5.16) 

a a a 

A time-varying resistor, on the other hand, can generate or accept 
net reactive power, and the sum of Eq. 5.16 is not conserved. The 
mechanism that produces the time variation of the resistance can 
be thought of as a source of reactive power. 

This situation is similar to that of reactances. A linear time- 
invariant reactance will not accept real power at any frequency; a 
nonlinear time -invariant reactance will accept real power at some 
frequencies and release it at others, although the overall real power 
input must vanish; and a time-varying reactance can generate or ac- 
cept net amounts of real power. The physical mechanism that varies 
the reactance can be thought of as a source of real power. 

5.3 Degenerate Elements 

Generally speaking, the degenerate devices discussed in Sec. 4.9 
obey frequency-power formulas of Types II and III in the same triv- 
ial way they obey formulas of Types I and IV. The ideal transformer 
obeys both Types II and III, with all P Q and Q a zero. The ideal 
gyrator obeys formulas of Type III (but not Type II) with all P a zero. 
The other degenerate elements are treated in a similar way. 

As mentioned in Sec. 4.9 (Ideal Transformer), a physical trans- 
former cannot pass d-c power. If it is placed in a circuit with a 
non-zero average voltage, it will saturate, and no longer approxi- 
mate an ideal transformer. Nonlinear resistors often have d-c 
power input, and if a transformer is placed in a circuit with them, 
some alternate d-c path must be provided if the transformer is to 
avoid saturation. 
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It is interesting to note that a time-varying transformer, such as 
an adjustable autotransformer, obeys frequency-power formulas of 
Types 11 and III. This device obeys equations of motion like Eqs. 
4.58 and 4.59 with the turns ratio N a function of time: 



and 



e 2 (t) = N(t)ej(t) 



i r (t) = - N(t)i 2 (t) 



(5.17) 



(5.18) 



It is easily verified that this device is lossless and linear, but since 
it is time-varying the power input at each frequency P a need not 
vanish. The device does not obey frequency-power formulas of 
Types I or IV, but does obey Types II and III. Thus, combinations 
of variable transformers and potentiometers [87] obey these form- 
ulas also. 

To show that this device satisfies formulas of Type III, we con- 
sider e 2 and i ]L to be functions of 6j and i 2 , given by Eqs. 5.17 
and 5.18. The matrix 



ai, 



9i 



8e, 



'i 

N(t) 

- N(t) 



(5.19) 



is positive definite, so the proof of Sec. 3.4 indicates that 

2 



(5.20) 



and any of the other forms of frequency-power formulas of Type III 
hold, where P is the power input at frequency co a , 



P a = 



(5.21) 



and where o> a is, as usual, an independent frequency on which the 
time variation of N(t) does not depend. 

To show that this device satisfies formulas of Type II, we consider 
e and -i to be functions of e 1 and i z t as in Eqs. 5.17 and 5.18. 
Nothing that 



= N(t) 



(5.22) 



we define the time-varying dissipation function 
G( ei , i t , t) = 



(5.23) 
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independent of the path of integration. The use of this state function 
yields frequency-power formulas of Type II, 



where Q a is the reactive power input at frequency w , 

Q Q = ^[(e^W + teAW] (5-25) 

In a similar way one can show that a time-varying gyrator, that 
is, a device that satisfies equations of motion of the form of Eqs. 
4.61 and 4.62 with R a function of time, obeys frequency-power 
formulas of Type III but not Type II. 

5.4 Network Extension 



If we have a network of elements, each of which individually sat- 
isfies frequency-power formulas of Type II (or III), it is of interest 
to ask whether the entire network does, in such a way that the powers 

P and Q can be evaluated at the ports of the network. The answer 

a a . 

is that it does, and the proof of this is identical to the proof given in 

Sec. 4.10, for networks of elements that obey formulas of either Type 
I or IV. 

5.5 Linear Time -Invariant Lossless Elements 

Generally speaking, the devices that obey frequency-power form- 
ulas of Types II and in are dissipative devices, or limits of dissi- 
pative devices. However, because linear time-invariant lossless 
elements cannot perform frequency conversion, the power input to 
such devices at each frequency, P Q , is zero. They therefore 
satisfy frequency-power formulas of Type III in a trivial way. 



Chapter 6 
DISTRIBUTED SYSTEMS 



In this chapter we deal with frequency-power formulas of Type I 
only, except that in Sec. 6.1 (Type III Formulas) formulas of Type 
III are derived. Formulas of Types II and IV do not apply because 
the Q a quantities appearing in the formulas cannot be interpreted 
as reactive powers. 

Distributed systems are those with important space variation. 
Energy is stored and dissipated throughout the space, but power 
enters and leaves only at the boundary, through the ports. Notice 
the similarity to a network of lumped elements, in which energy is 
stored and/or dissipated throughout the system, but power enters 
only at the ports. In Sees. 4.10 and 5.4 it is shown that, for the pur- 
pose of applying frequency-power formulas of Types I and III, the 
powers P^ can be evaluated at the ports, rather than at each ele- 
ment in the network. Here we show the same for some distributed 
systems. 

The quantities of engineering interest for frequency conversion in 
a distributed system are the powers entering and leaving the system 
at each frequency. However, the variables that describe a distrib- 
uted system are functions of time and space. The constitutive rela- 
tions and energy state functions are defined in terms of these vari- 
ables, with explicit space (and perhaps time) variation. The P a 
quantities in the frequency-power formulas are functions of space 
(but not time); their interpretation is not clear. 

Suppose, however, that each P is a divergence 






P a = - 



where F a can be expressed in terms of the variables of the system. 
Then the frequency-power formulas of Type I, 

P a 8w 
S a ^ = (6.2) 

or Type III, in any of the forms discussed in Sec. 3.4, such as 
2 



whichever are known to hold, can be integrated over the volume of 
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the system, P^ being the only function of space. The result is (by 
use of the divergence theorem) formulas of Type I 



n dS 

9cj 



= (6.4) 



or similar formulas of Type III, such as 

2 



-4S>0 (6.5) 

In Eqs. 6.4 and 6.5, the surface integral is taken over the entire 
boundary of the system, and n is a unit vector pointing inward, per- 
pendicular to the surface. It is our intent that each "F be inter- 
preted as the contribution to time-average value of the power-flow 
vector F caused by frequency w a ; thus the formulas of Eqs. 6.4 
and 6.5 relate quantities of engineering interest, the powers at var- 
ious frequencies, at the ports of the system. 

We conclude, then, that when the P a are in the form of diver- 
gences, the frequency-power formulas relate quantities that can be 
evaluated at the boundary of the system, the ports. Our problem 
then reduces to expressing the P as divergences. For some phys- 
ical systems this can be done in a simple way, directly from the 
equations of motion. Typical are the nonlinear electromagnetic 
medium discussed by Haus[35], treated here in Sec. 6.1; the irro- 
tational acoustic system discussed by Wagner [110], treated here in 
Sec. 6.2; and the irrotational electron beam (Haus [36] has discussed 
the case of most importance, the longitudinal beam with large pump 
and small signals) of Sec. 6.3. 

For systems in which it is difficult to express the P in the form 
of divergences a systematic procedure would be very helpful. Such 
a procedure is available for systems that obey Hamilton's principle 
[30,Ch.lI; 70, Part I, Ch. 3, Sec. 3.4], but it is quite complicated. 
This is not a drawback, though, because once this systematic pro- 
cedure has indicated the form of the expressions, the results can 
then be derived directly from the equations of motion. A gyromag- 
netic medium obeys the full Manley-Rowe formulas, as we show in 
Sec. 6.4. This fact was originally proved by McWhorter [64] using 
this systematic procedure (Haus [35] had earlier discussed the case 
of most practical importance, large pump and small signals), but 
now that the proper variables are known, a simpler proof can be 
used. 

Hamilton 1 s principle is discussed in Sec. 6.5, along with the sys- 
tematic procedure leading to the Manley-Rowe formulas. In Sec. 
6.6 we discuss isentropic rotational flow of a charged, compress- 
ible fluid. This system is so complicated that Hamilton's principle 
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actually simplifies the presentation. Another complicated system, 
rotational fluid flow with free currents, is discussed in Sec. 6.7; it 
is an approximation to a magnetohydrodynamic system with infinite 
conductivity. 

One difficulty with this approach is illustrated by these examples. 
Although the systematic procedure based on Hamilton 1 s principle 
guarantees that the P a are in the form of divergences, it does not 
guarantee that physically unmeasurable potentials or Lagrange mul- 
tipliers do not appear in the power expressions. Because of this, 
there is often a real problem in assigning a physical interpretation 
to the P a or T^ . 

The form of the expressions for power, when they contain such 
physically unmeasurable terms, cannot be guessed from a knowledge 
of the physical power-flow vector. Although any result obtained by 
use of Hamilton's principle can also be obtained directly from the 
equations of motion, we would not ordinarily be motivated to define 
the required physically unmeasurable terms, which are provided 
automatically by the systematic procedure. The Hamilton's prin- 
ciple method, discussed in Sec. 6.5, is complicated, but its use is 
justified for systems that cannot be treated directly. 

6.1 Electromagnetic Medium 

Here we show that the electromagnetic field in a nonlinear, dis- 
persion-free anisotropic, inhomogeneous, reciprocal, stationary 
electromagnetic medium obeys the Manley-Rowe formulas. This 
fact was first shown by Haus [35], and our results, in Sec. 6.1 (Man- 
ley-Rowe Formulas), are not significantly different. In Sec. 6.1 
(Type III Formulas) it is shown that electromagnetic material with 
nonlinear loss and linear energy storage obeys formulas of Type III. 

Equations of Motion. The variables of this system are the elec- 
tric field "E, the magnetic field H", the electric displacement vec- 
tor "D, the magnetic induction vector B", and the electric current 
density T. Each is a function of space and time. We assume that 
the medium satisfies three constitutive relations: one between "E" 
and I); one between B" and H; and the third between E and J. 

We assume that the magnetic field H is a single-valued function 
of B" at each point in space. We do not allow H" to depend on the 
time or space derivatives of B", nor on time explicitly, but only on 
the instantaneous value of B" . Thus, although we do not allow dis- 
persion or hysteresis, the functional dependence may be anisotropic 
and nonlinear, and may depend on space; for example, the material 
may be inhomogeneous, or may occupy only a portion of the system. 
We also do not allow the material to move (that is, it must be sta- 
tionary), and we require that the dependence be "reciprocal, M by 
which we mean (if Cartesian co-ordinates x, y, and z are used) 

8H 

< 6 ' 6) 
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3H 9H 



8B " 8B 
z y 



and 

8H, 



z 



9B X " 815 ' 

A. Zi 

We have, then, the constitutive relation 

H = H(B , r) (6.9) 

where we show the explicit dependence on space; because this is re- 
ciprocal we may define the "magnetic energy density" by the line 
integral [42, Sec. 6.10; 80] 



W m (B, r) = \H(B, r) - dB (6.10) 

<J 

In a similar way, we assume that the electric field E" is a single- 
valued, nonlinear, dispersion-free, anisotropic, inhomogeneous, 
reciprocal function of 15 

E = E(D, 7) (6.11) 

so that the line integral 

W e (D, 7) = \E(D, 7) dD (6.12) 



defines the electric energy density. 

In portions of the system occupied by free space, the constitutive 
relations are 

H = B" (6.13) 



and 



E ss D (6.14) 



where u is the permeability of free space and e is the permit- 
tivity of free space, and so the energy functions are 

W m = Y~ (6.15) 

and 

W e = ^- (6.16) 



Ohrn^s law is the constitutive relation between the current density 
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7 and the electric field 1?. We assume that 7 is a single-valued, 
nonlinear, dispersion-free, anisotropic, inhomogeneous function of 
"E (but one that need not be reciprocal), 

7 = J(E, 7) (6.17) 

The power dissipated per unit volume is 7 - "E; if the system is 
lossless we assume either 7 or 1C vanishes at each point. 
The equations of motion are Maxwell 1 s equations 

V - B = (6.18) 

V - D = (6.19) 

VXH = T+i (6.20) 
and 

V X U = -B" (6.21) 

These four equations are linear, but the constitutive relations for 
the system need not be. 

From these equations of motion and the definitions of the energy 
state functions one can derive Poynting's theorem (dot-multiply Eq. 
6.20 by "E and subtract it from Eq. 6.21 dot-multiplied by H) 



= - E - 7- V- (EXH) (6.22) 

This is customarily interpreted by identifying W fi + W with the 
energy density, E J with the power dissipated per unit volume, 
and EXH with the power -flow vector. Thus, the power dissipated 
and power going to energy storage must be supplied by the power- 
flow vector "E X H. 

Because Maxwell 1 s equations are linear, we may derive similar 
theorems involving only Fourier coefficients of the variables. Let 
us express each variable in the form of Eq. 3.15, where each Four- 
ier coefficient is a complex vector (a set of three complex numbers, 
or a real vector and an imaginary vector). These complex vectors 
depend on space, but not on time. Then since Eq. 6.21 holds for 
the actual vectors "E and "B", a similar equation holds for each set 
of Fourier coefficients; 

VXE a = -j<o a B a (6.23) 

Similarly, Eq. 6.20 implies equations of the form 

VxH a* =7 a*-Jaa* < 6 - 24 ) 

If we dot-multiply Eq. 6.24 by *E a and subtract it from Eq. 6.23 dot- 
multiplied by H * , we obtain 

V - (E Q X H*) = - jco a (B a . H a * - E a - D Q *) - E a . J Q * (6.25) 
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For future reference we note two things. First, if the medium is 
lossless, then either the current density vanishes and all T = , or 
the electric field vanishes and all "E a = ; in either case, the last 
term in Eq. 6.25 disappears. And secondly, if the energy storage 
constitutive relations are linear (if B" is a linear function of H" and 
D is a linear of E), then (B a - H a * - E a I5 Q *) is real, and the first 
term on the right-hand side of Eq. 6.25 therefore has no real part. 
The second statement follows from the linear constitutive relations. 
For example, the relation between B" and H must be of the form 

B = ^Hj (6 .Z6) 

where B i and H. are components of the vectors B" and H", and 
where JJL^. is the permeability tensor, and where a sum over the re- 
peated index j is implied. Because B" is a reciprocal function of 
H tKe permeability tensor must be reciprocal, so p... = JJL.. . Thus 



= < B j>a*< H j>a 



. H (6.27) 



Q 



and (B a ^ a *) * s real. In Eq. 6.27 we sum over the repeated in- 
dices i and j . Similarly, (E a "D^) is real. 

Manley-Rowe Formulas. Suppose the material is lossless; then 
T? a * "3^* vanishes. Jf we use W e + W m as an energy state function 
for the proof of the Manley-Rowe formulas, we obtain formulas like 
Eq. 6.2 with the powers P given by 



P a = 



- E Q D a *) (6.28) 



and thus by Eq. 6.25 



P a = - V ' a < 6 * 29 ) 

where 

P a = -|ReE a XH a * (6.30) 

The Manley-Rowe formulas can then be written in the form of Eq. 
6.4, and the f quantities can be interpreted as the power -flow vec- 
tor 1C xH at frequency c. We evaluate the powers at the boundary 
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of the system, as indicated in Eq. 6.4. 

Type III Formulas. The system obeys formulas of Type III pro- 
vided the energy storage is linear. Equation 6.25 becomes 

iRe (E a * J a ) = - V- F Q (6.31) 

where F" is given, by Eq. 6.30, and has the interpretation given 
above. 

The constitutive relation of Eq. 6.17 leads to frequency -power 
formulas of the form of Eq. 6.3 with 

P a = |Re (E a * . J a ) (6.32) 

Because of Eq. 6.31, these formulas can be integrated to the form 
of Eq. 6.5, and we can evaluate the powers on the boundary of the 
system. 

6.2 Ir rotational Fluid Flow 

In this section it is shown that the irrotational flow of a 
non-viscous fluid obeys the Manley-Rowe formulas. The concept of 
irrotationality is based on a well-known theorem proved in Appen- 
dix F, which states that if a system without viscosity starts flowing 
so that the curl of the momentum vanishes, then the momentum re- 
mains curl-free. Since viscosity spontaneously introduces a curl of 
the momentum, we cannot speak of irrotational flow of a viscous 
fluid. 

For fluid flow there is a choice of coordinates to use. The cus- 
tomary coordinates space and time are known as Eulerian co- 
ordinates; they are not always the most suitable. Bobroff [3] has 
discussed three other sets of coordinates, some well suited for small- 
signal descriptions of electron beams. One that is widely used in 
fluid flow analysis is the Lagrangian set, which consists of time 
and the initial position of the particles, whereas the one that is best 
suited for a description of fields is the Eulerian set. Various prob- 
lems call for various sets of coordinates, and it is of interest to ex- 
press the Manley-Rowe formulas in terms of these various sets. 

In Sec. 6.2 (Eulerian Formulation), following the work of Wagner 
[110], we prove the formulas using Eulerian coordinates. The ex- 
pressions for power flow are evaluated on a boundary that is fixed in 
space. Material may pass through this boundary, and in applying the 
formulas this must be taken into account. The resulting formulas 
are suited for study of stability of hydrodynamic systems, but they 
are not as appropriate for problems in nonlinear acoustics, because 
in coupling to acoustic transducers the powers should be evaluated 
over a boundary that moves in space with the particles. The La- 
grangian coordinates are thus the natural ones to use; in Sec. 6. 2 
(One -Dimensional Lagrangian Formulation) we prove the formulas 
for a one -dimensional acoustic system using Lagrangian coordinates. 
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The expressions for power flow are different in the two cases. 

Undoubtedly the Manley-Rowe formulas can be expressed in other 
sets of coordinates also. 

Eulerian Formulation. The variables of the fluid-flow system 
are the particle density p , the particle velocity v", the momentum 
per particle p, and the pressure TT . It is assumed that each of 
these can be defined as a continuous function of space and time, mac- 
ro s c opically . 

These are related by the constitutive relations of the fluid. The 
velocity and momentum of each particle are related by an expression 
of the form 

v = v(p) (6.33) 

and therefore we can define a new variable, the kinetic energy per 
particle T , by the line integral 

T = ^vfp) - dp (6.34) 

Thus, any variation of T is related to the variation of "p 

6T = v 6p" (6.35) 

and in particular 

VT = (v V)p+vX(VXp) . (6.36) 

Normally the relation of Eq. 6.33 is 

"p = mv (6.37) 

where m is the mass of the particle, but for relativistic flow this 
is modified to the form of Eq. 4.41. There is no need to specify 
what the relation of Eq. 6.33 is, so we do not. 

Another constitutive relation is that between the pressure and the 
particle density 

tr = -rr(p) (6.38) 

from which we can calculate the intrinsic energy per particle by the 
integral 

T = (- dp (6.39) 

so that 

V (T + - ) = (6.40) 

P P 

where (T + ) is the enthalpy per particle. 

One of the equations of motion is the force equation, 



= p- + v 

= Force per particle 

= -3U*L (6.41) 
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where D/Dt is the substantive derivative (the time derivative taken 
in a reference frame moving along with the particles), and where F V 
is the viscous force. From Eqs. 6.36, 6.40, and 6.41 one obtains 

V(T +T + ^) = -t + vX(VXp) + F V (6.42) 

which is an alternate form of the force equation. 

For systems with no viscosity and irrotational motion, we set 1^. 
and V X p to zero, obtaining 

V(T + T + ~) = -f (6.43) 

The other equation of motion is the continuity equation, or law of 
conservation of particles, 

V- (pv) = -p (6.44) 

From these equations of motion and the definitions of the variables 
T and T one can prove the power-flow theorem 



= -V- (pv) (T + T + ) (6.45) 

This is interpreted by identifying pT + pT with the energy density, 
and (pv) (T + T + j) with the power-flow vector F" . The system is 
lossless; the energy density can increase only insofar as the power- 
flow vector has a nonzero divergence. The energy density is the sum 
of kinetic energy and intrinsic energy per unit volume, and the power- 
flow vector has three terms: the term pvT is transport of kinetic 
energy by flow of particles; the term pvT is transport of intrinsic 
energy by flow of particles; and the final term, -rfv, represents a 
flow of power caused by movement of particles overcoming a pres- 
sure force. 

Observe that the power -flow vector "F = pv"T + pvT + TTV is not 
written, except for the last term, as a product of two variables. The 
electrical power -flow vector "EXH is written as a product, and so 
the concept of "power at frequency o> " is very useful, since it is 
clear how to split apart the power-flow vector into two factors. The 
fluid-flow terms can be split apart more than one way; one of the 
results of our analysis is to give physical significance to one of these 
ways. 

We may derive an expression similar to Eq. 6.45 involving the 
Fourier coefficients of the variables. We write the variables T, 
T, , "p, p, and (pv) in multiple Fourier series as in Eq. 3.15. 
Then Eq. 6.43 implies a similar relationship among the Fourier co- 
efficients, 
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(6 . 46) 



and Eq. 6.44 implies 

V- ( P v) a = -Jfc> Q Pa (6.47) 



If we multiply Eq. 6.47 by [T a * + X a * -f (?) a *J * and dot-multiply 
Eq. 6.46 by (pv) , and add, we find 



(6.48) 

The Manley-Rowe formulas now follow easily. If we consider 
pT + pT to be an energy state function depending on p and "p, the 
Manley-Rowe formulas are written in the form of Eq. 6.2 with the 
powers P given by 

*>a - i** J a [\* + X a * + (i) a *] P Q + (pv) a * - 



{[ 



T a* + V + ?* " ' <^ ' 



or, by Eq. 6.48, 

P Q =-V-T r a (6.50) 

where the vector IF is the power-flow vector at frequency ji> a , 






F a , iRe ( P v) Q T Q * + T Q * + a * (6.51) 

Thus the formulas can be written in the form of Eq. 6.4, and the 
powers can be evaluated on the boundary of the system. 

One -Dimensional Lagrangian Formulation. In a one -dim en si onal 
system the variables depend on only one coordinate, x. At some 
initial time, t , we label each particle with its position X Q , and as 
the fluid moves, the particles retain a knowledge of their initial po- 
sition, X Q , even though they move away from it. The variables of 
the system, p, v, pv, T, T, etc., can be expressed either in 
terms of actual position x and time t, (the Eulerian set of coor- 
dinates) or in terms of the initial position X Q of the particle and 
time t (the Lagrangian set). For notation, we denote partial de- 
rivatives with respect to the Eulerian coordinates by the 8 symbol, 
and partial derivatives with respect to the Lagrangian coordinates 
by the D symbol. The time -differentiation D/Dt is the customary 
substantive derivative, 

= JL +V .V (6.52) 
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We must express the equations of motion in terms of the Lagran- 
gian coordinates. Eq. 6.41 becomes (neglecting viscosity) 

I>P -. 1 iZL - _ J- 2L (6.53) 

Dt P 3x P Dx Q 

and Eq. 6.44 implies 

JL(! 
Dt\P 



JL Dv 

p o ^o 



(6.54) 



where p Q , the initial particle density, is independent of time in the 
Lagrangian set of coordinates. 

The kinetic energy per particle, T, is a function of p given by 
Eq. 6.34, and we consider the intrinsic energy per particle to be 
given by 



which is equivalent to Eq. 6.39, except that it is now a function of 
(1/p) , instead of p , 

With the aid of Eqs. 6.53 and 6.54 we can compute the power-flow 
theorem 



D/l 



DIT Dv 

V - IT- 



Dx Dx 
= -~-~-(TTV) (6.56) 

which is interpreted by identifying p Q (T+T) with the energy density, 
and TTV with the power-flow, using Lagrangian coordinates. 

We now express each of the variables in a multiple Fourier series, 
where the Fourier coefficients depend on X Q , and rewrite Eqs. 6.53 
and 6.54 as 

j v> P Q * = \* ( 6 - 5? ) 



and 

JVo (?) a = ~ v a (6>58) 



If the energy function p (T + T) , as a function of p and -p, is 
used to generate Manley-ilowe formulas, they are written in the 
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form of Eq, 6.2 with the expression for P given by 



(6.59) 

which is, in view of Eqs. 6.57 and 6.58, 



where F a is the component of power-flow TTV caused by frequency 



co 

a 



F = iRe ir*v (6.61) 



a 



The Manley-Rowe formulas of Eqs. 6.2 can now be integrated over 
the length of the system, and in view of Eq. 6.60 the power can be 
evaluated at the ends of the system, as the power input at frequency 
co . Notice that F is evaluated by following along the particle at 
the end of the system, and noting its pressure and velocity as a 
function of time. Alternatively, if there is an acoustic transducer 
at the end of the system, the power may be evaluated as the force 
times the velocity of the diaphragm of the transducer. This Lagran- 
gian formulation is therefore well suited for nonlinear one-dimen- 
sional acoustical systems. 

We have attempted, without success, to find Manley-Rowe form- 
ulas for a three-dimensional acoustic system, using Lagrangian 
coordinates. It would be of interest to know under what conditions 
such formulas exist. 

6.3 Irrotational Electron Beam 

In this section, we show that an irrotational electron cloud in 
arbitrary electromagnetic fields obeys the Manley-Rowe formulas. 
One device this treatment applies to is the longitudinal electron 
beam with heavy pumping and small signal levels; this was discussed 
byHaus[36], who derived the Manley-Rowe formulas. His result 
is generalized here to irrotational beams and arbitrary signal levels. 

The concept of "irrotational flow" is based on a well-known theorem 
proved in Appendix F, which states that if a beam starts out with 
zero curl of the momentum (with electromagnetic fields, this is the 
"canonical" momentum [30, Sec. 2, 6]), it continues to move with 
zero curl [28], Some practical conditions under which this holds 
[28] are discussed in Appendix F. 

This system is, in a sense, a combination of the previous two 
systems, in which the particles have a charge q. For simplicity 
we neglect pressure effects, although they can be included very 
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easily [110] . Although, we speak as though the particles were elec- 
trons, the analysis clearly applies to any charged particles. In ad- 
dition, it can be extended very easily to a system with two or more 
species of particles, each with its own mass and charge, provided 
each individually satisfies its own continuity equation, and provided 
they interact only through long range electromagnetic forces. Thus, 
in particular, macroscopic plasma models without viscosity and col- 
lisions might be included. 

Our result is an extension of Chu's Kinetic Power Theorem [15], 
which is known to hold for longitudinal [41] and irrotational [48] elec- 
tron beams. The gist of this theorem, in our present notation, is 
that for small signals each P individually vanishes.lt is a small- 
signal power theorem, appropriate for small variations of a time- 
invariant device, with no frequency conversion taking place. It 
breaks down, however, when the disturbances are so large that non- 
linear effects occur. The Manley-Rowe formulas then express the 
constraints that remain in the nonlinear case. 

We use Eulerian coordinates to describe the system, even though 
many electron beam analyses are done using other sets [3] . The 
fact that small -signal power theorems are available in other sets of 
coordinates [39, 100, 4] leads one to suspect that the Manley-Rowe 
formulas can also be expressed in these other coordinates, although 
we do not discuss them here. 

The variables of this system are the electromagnetic field quan- 
tities "E, II, "D, and B", and the particle variables p, v, and 
p. The field variables "E and B" are related by two of Maxwell's 
equations, Eqs. 6.18 and 6.21; therefore, we may define, in the usual 
way [30, Sec. 1.5], the vector and scalar potentials A and $ so that 

B" = VXX (6.62) 

and 

E = -VS-X (6.63) 

The momentum per particle, p, is the canonical momentum [30, 
Sec. 2. 6], which is, nonrelativistically, 

"p s mv -f qA (6.64) 

where m is the mass of each particle. The relativistic form for 



p is 



P = 



mv 4- qA (6.65) 

Wl-( v /c) 2 

where m is the rest mass of the particles. In each case, the veloc- 
ity is a function of (p - qA), and the kinetic energy per particle is 



T = v(5 - qA) d(p - qA) (6.66) 
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so that any variation 6T of the kinetic energy is related to varia- 
tions of "p and JL 

6T = v dp - qv 5A (6.67) 

and in particular 

VT = (v V)p - q(v V)A + v X (V X p) - qv X (V X A) (6.68) 

The customary constitutive relations for the electromagnetic var- 
iables are those for free space, Eqs. 6 13 and 6.14. However, to 
allow for the possible presence of nonlinear, dispersion-free, aniso- 
tropic, inhomogeneous, lossless, reciprocal, stationary media in 
some parts of the system, we merely require relations of the form 
discussed in Sec. 6.1, so that the electric and magnetic energy den- 
sities W and W can be defined as in Eqs. 6.12 and 6.10. 
e m 

The equations of motion include Maxwell's equations 

V - 1" = (6.69) 

V D = qp (6.70) 

V X H = qpv +T> (6.71) 
and 

VX U = - 1" (6.72) 

where we have identified the currents and charges as those of the 
electrons* 

The mechanical equations of motion include the force equation, 

JQ JL ~ 

(p - qA) = p - qA + (v V)p - q(v V)A 
= Force per particle 

= qE + qvXB" (6.73) 

which, with the aid of Eqs. 6.62 and 6. 68, becomes 

VT = -"5 + qA + qE + v X (V X p) (6.74) 

If the motion is irrotational, we set V X "p = , to obtain 

* 
VT = - p + qA + qE (6,75) 

The other mechanical equation of motion is the continuity equation, 
V (pv) = - p (6*76) 

From these equations of motion and the definitions of the energy 
state functions one can derive a power-flow theorem 

"1 < W e + W m + pT) = T P + (P^ ' ^ - fapv) * 5 + "E ' D + H 5" 
=-V ("E X H + pvT) (6.77) 
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We interpret tMs by identifying W e + W m + pT with the stored en- 
ergy density, and "E X H + pvT with the power-flow vector F , 

In the same way that Eqs. 6.25 and 6.48 were derived, we can 
find a formula similar to Eq. 6.77, involving Fourier coefficients 
of the variables. It is 

- A a - E Q . D a * + H* B Q ] 
= -V EXH* + (pv)*T (6.78) 



The Manley-Rowe formulas are now proved using W" e + W^ + pT 
as the required energy state function, as a function of the variables 
"D, B", p, and (p - qA) , The expression for P a is simply a com- 
bination of the expressions obtained in Sees. 6,1 (Manley-Rowe Form- 
ulas) and 6,2 (Eulerian Formulation): 

P = - V (6.79) 

a a 

where 

(6.80) 



is the contribution to time -average power-flow vector F = E X H + pvT 
caused by frequency u> . Because the Manley-Rowe formulas then 
are written in the form of Eq. 6 A, the powers can be evaluated at 
the boundary of the system. 

6.4 Gyromagnetic Medium 

We shall now show that a gyromagnetic medium obeys the Manley- 
Rowe formulas. The case of most practical interest is that of large 
pump and small signals, for which Haus [35] has proved the formu- 
las. The more general case was first done by McWhorter [64], who 
used the Hamilton 1 s principle approach to determine the correct var- 
iables to use. More precisely, McWhorter used the Hamiltonian 
state function, instead of the Lagrangian, with Hamilton 1 s canoni- 
cal equations as the equations of motion. Because the Hamiltonian 
density and the Lagrangian density are related by a Legendre trans- 
formation, the discussion of Sec. 4.1 shows that the same P is ob- 
tained with each. The important contribution was recognition of the 
appropriate variables. When these are known, the formulas can be 
proved directly from the equations of motion, as we show here. 

One often uses the equations of motion of a gyromagnetic medium 
to describe ferrites [2, 102], if losses are neglected. Although the first 
microwave parametric amplifiers were made from ferrites [103, 114], 
not until I960 was it known that the gyromagnetic model obeys the full 
Manley-Rowe formulas. 

The variables of the gyromagnetic medium are the electromagnetic 
field quantities E, H, D, and B, and the magnetization vector *M. 
The constitutive relation between D and E is assumed to be 
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D = eE (6.81) 

where e is the dielectric constant of the material. We can equally 
well assume that the relation between D and E is single -valued, 
nonlinear, dispersion-free, anisotropic, inhomogeneous, and recip- 
rocal, as in Sec. 6. 1. The important point is that the electric 
energy density 

W e (D, r) = CE(D, Y) dD (6.82) 

be defined. 

The constitutive relation between the magnetic vectors, on the other 
hand, is (using, as usual, M.K.S. units) 

B" = ii (H + M) (6.83) 

where [i Q is the permeability of free space. 

The equations of motion include the four Maxwell's equations, 

V B" = (6.84) 

V " D" - (6.85) 

VXH = 1) (6.86) 

and 

VXE = -B" (6.87) 

and the magnetization equation, 

M = -Y(MXB) (6.88) 

where Y is a constant. Usually this relation is written with B re- 
placed by S, but by Eq. 6.83 we see that to write it using B we 
need merely redefine the constant Y . 

It is clear from Eq. 6.88 that the magnitude of M is constant, for 

-|r (M 2 ) = 2M M = (6.89) 

9t 

Thus only two variables are necessary to determine M, not three. 
The recognition of the proper variables to use [64] is the key to 
proving the Manley-Rowe formulas. We use the component of M 
in an arbitrary direction (call it the z-direction), M Z , and the angle 
<j> made by the projection of 3J in the x-y plane. Thus 

M x = N/M 2 - M z 2 cos 4> (6.90) 



M = N/M 2 - M z 2 sincj> (6.91) 

and 

M z = M z (6.92) 
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We now express Eq. 6.88 In terms of M Z and <f> . First, we ob- 
serve from Eqs. 6.90 and 6,91 that 

3M 8M x 8M y 

B * = B + B 



= (MXB) z ( 6 -93) 

and therefore, from Eq. 6.88, 

M z = _ Y B - -^ (6.94) 

Also, from Eq. 6.88, 

and therefore, by Eq. 6.94, 

5> - VB" 8 "fe (6.96) 

Equations 6.94 and 6.96 are the equations of motion expressed in 
terms of the variables <j> anc * ^ z * 

In proving the Manley-Rowe formulas, we use the state function 

R 2 

U = W e (D , r) + y_ - E* M (6.97) 

which is a function of T5 , B" , cj> , and M Z . It is clear from Eqs. 
6.94 and 6.96 that 

M 
12 - 5. (6.98) 

a<t> \ 

and 

* 

The Manley-Rowe formulas generated by U are written in the form 
of Eq. 6.2 with 

P = -=rRe jcol "E * "D + ^- B" a - "M a * B" 
a ^o 
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where 

F a = fRe (E a XH a *) (6.101) 

where we have used Eqs. 6.86, 6.87, 6.98, and 6.99, written in 
terms of the Fourier coefficients alone. The Manley-Rowe form- 
ulas can now be integrated over the system and put in the form of 
Eq. 6.4. The power quantities are the ordinary electromagnetic 
power flows at each frequency. 

6,5 Hamilton's Principle 

The distributed systems discussed in Sees. 6.1, 6.2, and 6.3 are 
all so simple that the Manley-Rowe formulas can be proved directly 
from the equations of motion. This cannot be done easily for all dis- 
tributed systems, however, and a systematic procedure for obtaining 
the formulas is helpful. Such a method is available for systems that 
obey Hamilton's principle; in this section we demonstrate that Ham- 
ilton^ principle leads to Manley-Rowe formulas with the powers P a 
automatically in the form of divergences. This technique is applied 
to physical systems in Sees. 6.6 and 6.7. 

Hamilton 1 s principle for distributed systems [30, Ch. 11; 70, Part I, 
Ch. 3, Sec. 3.4] is not as well known as the corresponding principle for 
lumped systems. It is a cumbersome way of obtaining the equations of 
motion, and invariably is stated only after the equations of motion are 
known. Its major uses in the past have been as a prelude to extending 
known equations of motion, either quantum- mechanically or relativist!- 
cally. We use it here to generate Manley-Rowe formulas; a similar use 
in the past has been to generate, in roughly the same way, small- signal 
power theorems [100, 57], 

For most systems the kinematic constraints are nonholonomic[30, 
Sec. 1.3], and either potentials must be defined, or Lagrange multi- 
pliers used. There is no assurance that these physically unmeasur- 
able quantities do not appear in the final expressions for the P a . 

The full Hamilton's principle method, with several variables in 
three dimensions, is complicated, so we give first, in Sec. 6.5 (A 
Simple Case), a simple derivation with only one variable and only 
one spacial dimension. The concepts are all illustrated by this sim- 
ple case; the extension to any number of variables, and explicit 
space and time variation of the Lagrangian, in three dimensions, is 
straightforward, and is done in Sec. 6.5 (Extension). 

Although we use Eulerian coordinates exclusively, there is no 
reason why Hamilton's principle cannot be formulated[lOO] in other 
coordinates [3] . 

A Simple Case. In distributed systems the energy state functions 
are densities; the Lagrangian density is used in Hamilton's princi- 
ple. One of the variables of the system, 7) , together with its space 
and time derivatives, defines the Lagrangian 
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L = Lfo, ri, r? x ) (6.102) 

where we denote, temporarily, 8r?/8x by TJ X . Hamilton's principle 
states that the integral of L over both time and space is stationary 
with respect to a variation of rj that vanishes on the boundary of the 
system and at the initial and final instants, or 



tt 



L dt dx = (6.103) 



We introduce here some convenient notation. We need symbols 
for four types of differentiation: 

(1) Regular partial differentiation of a variable or function 
with respect to time or space coordinates, as 8/8x, 
8/8y, 8/8z, and 8/8t. The dot notation is reserved 
for this partial derivative with respect to time, and the 
vector operator V is used for these space partial de- 
rivatives where convenient. 

(2) In lumped systems, where t is the only independent 
coordinate, we have the total derivative d/dt. 

(3) In fluid work, we use the substantive derivative D/Dt 
= 8/8t + (v V), where v~ is the particle velocity, to 
indicate time differentiation in a reference frame mov- 
ing along with the particle. 

(4) We need notation to indicate partial differentiation of 
an energy state function with respect to one of its vari- 
ables, keeping the other variables fixed. Generally, 
such functions depend not only on the variables of the 
system, but on the independent coordinates (such as 
time) as well. Thus the usual 8 notation is ambigu- 
ous. We denote such differentiation by A; for ex- 
ample, A/Arj means the partial derivative with re- 
spect to 77, keeping the other variables on which the 
function depends constant. Thus the two quantities 
AL/At and L = 8L/8t are quite different. 

Using this notation, any variation in L is, to first order, 



One more convenient notation is the symbol * which will mean 
that the difference between the quantities it separates integrates out 
to zero. Thus, for example, for any scalar , 

6r? = - 6r) (6.105) 

since 



J 



|6r?)dt = -(|6T7) dt (6.106) 
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vanishes, inasmuch as 677 vanishes at times tj and t 2 . Similarly, 
V (677) = (6.107) 

for any vector ![ . 

Using this notation, Hamilton's principle states that 

6L = (6.108) 

or, 

AL AL . AL 

= 6L = - 677 + -7-7 6r? + - - 677 
Ar? Ar} ' Ar? 'x 



The quantity on the last line of Eq. 6.109 can integrate to zero only 
if, for all time and for all space, the "Euler-Lagrange" equation 

AL _ 3 AL 8 AL 



remains valid. 

The Hamiltonian state function I is defined by the Legendre 
transformation 

I = ^ 77 - L (6.111) 

and so, by using Eq. 6.110, we obtain 
* - JLL - JL ( AL *\ 9L 

1 = at ~ at \ AH n } ' Bt 

* JL Ak AL >, AL . AL ., AL . 
9t Ar] Ar5 Ar7 A 77 Arj x 



If the Hamiltonian is to be identified with the energy density (as it 
often can be), Eq. 6.1 1Z is an expression of the law of conservation 
of energy, with 

AL . (6.113) 



playing the role of the power-flow vector (x- component). 

To derive frequency-power formulas of Type I, we use either L 
or H as the required state function. The formulas are in the form 
of Eq. 6.2, with the P a given by 
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But because Eq,, 6,110 predicts 

L\* 8 




I/LI, C X 

Q a (6J15) 

Eq. 6.114 becomes 

Thus each P a is put into the form of a divergence, and the form is 
expected from a knowledge of the actual power -flow vector of Eq. 
6.113. 

We have seen in this simple case that Hamilton 1 s principle leads 
to P a in the form of a divergence. This is extended in the next 
section. 

Extension. Here we extend the treatment of the previous section 
to describe systems with several variables, and with explicit space 
and time variations, in three dimensions. It is assumed there are 
a number of variables defined as a function of space "r" and time t . 
In case some of the variables are vectors, we can treat each com- 
ponent separately, or we can use vector notation. The Lagrangian 
depends on variables 77^ and on their derivatives 8r]j/8t and 
97?i/0xj as well as on space and time explicitly. 

Bri i 
L = Lfo., V^ r, t) (6.117) 

where (as in the future) i is an index running over the variables, 
and j is an index running over the three directions in space. 

Hamilton's principle states that the integral of L is stationary 
with respect to variations of the r). variables that are consistent 
with a set of "kinematic constraints. rl These constraints are often 
obvious from geometrical considerations for lumped systems, but 
in distributed systems they are far from obvious; it is necessary to 
introduce them purposefully. Virtually all kinematic constraints for 
distributed systems are nonholonomic [30, Sec. 1.3]; they are taken 
account of by two methods. 

One is the definition of potential functions. For example, sup- 
pose one constraint is V X A = for some vector variable A. We 
can account for this constraint by defining a scalar function \\t such 
that A = Vijj. We then express the Lagrangian in terms of \\i and 
its four derivatives, eliminating A as a variable. 

The other method is that of Lagrange multipliers. If the nonholo- 
nomic constraint is 
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we define a Lagrange multiplier \ and a new Lagrangian 

L 1 = L + \f (6.119) 

The two Lagrangians are, of course, equal numerically, but are 
quite different functional forms; thus, in general, 

AL , AL 1 ,, 1om 

-= (6.120) 

Ar}. Arj. 

etc. 

If there are several nonholonomic constraints, we may use the 
potential method on some and the multiplier method on the rest, if 
that is convenient. Usually neither the potentials nor the Lagrange 
multipliers are physically measurable quantities, although some- 
times a physical interpretation can be given to them. 

After the kinematic constraints (both holonomic and nonholonomic) 
have been put into the Lagrangian, either by direct elimination of 
variables, introduction of potentials, or addition of Lagrange mul- 
tiplier terms, we remain with a Lagrangian that is a function of some 
of the original 77. variables and their four first- order derivatives, 
some potentials and their four first-order derivatives, some La- 
grange multipliers, and space and time coordinates. We consider 
the potentials and Lagrange multipliers to be bona fide variables 
from now on, so that L is a function of some rj. , their first-order 
derivatives, and space and time coordinates, 



Using this Lagrangian, Hamilton's principle states that 

6L = (6.122) 



or, 



o , 6L = 2 6ii . 

i Ar?. i Arj. i j 

1 1 



ax. 
. .. (6-123) 

Ar?. 8t Ai7. j 9x. 



. . . 



8x. 
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The expression in Eq. 6.123 vanishes, when integrated over space 
and time, if and only if, for all i , 

L - A Ak + s J_ AL (6.124) 

Ar?. Bt AT}, j 3x. 877. 

1 x J A ^. 

8x. 

fox all time and at each point in space. Equations 6.124 are the 
Euler-Lagrange equations. The equations corresponding to the 
Lagrange multipliers are merely the kinematic constraints that the 
Lagrange multipliers were defined for. 

If Hamilton's principle is correct, the Euler-Lagrange equations 
should predict the ordinary equations of motion, and vice versa. It 
is not a trivial task to show that this is true. Equations 6.124 con- 
tain potentials and Lagrange multipliers; these must be eliminated 
to derive the ordinary equations of motion, and must be defined in 
any attempt to show the generality of the Euler-Lagrange equations* 

The Hamiltonian H is defined by the Legendre transformation 



and so by use of Eq. 6.124, the time -derivative of H is 



I _SL AL ALj - AL . _ AL __ 2 AL 7 i _ AL 

"" i i 3t An. An. i Arj. i Af}. i j 877. 9x. At 

- A ~^ JJ 

8x. 

S I n S - u *,-" y; ^-^ ^ i AL 

-* r - J 8x. 8n. " j 8r?. 8x. I " At 

A "8xT A "ixT 

J A J 

AL , , ^ (6J26) 



Tj dx. I 8r?. "i) At 
' ^ A 8^ 



so we can "write 

AL 



I + v . 



At (6.127) 

where the vector F" has components 

F ; = S ; -^- fl : (6.128) 
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If fi is the internal energy density, then we interpret Eq. 6.127 
as the law of conservation of energy, where IF is the power-flow 
vector. More often, however, H is not the energy density, and 
one must use care in interpreting 7 . 

We can obtain formulas similar to Eq. 6.124 but involving the 
Fourier coefficients of the variables, 




The Lagrangian state function generates Manley-Rowe formulas of 
the form of Eq. 6.2, with the powers P given by 




P a 



which, with the aid of Eq. 6.129, can be easily put into the form 

P a = - V '^a (6.131) 

where the vector "F has components 

< 6 - 132) 




and is therefore the contribution to time -aver age value of F caused 
by frequency oj^ . The Manley-Rowe formulas can now be integrated 
to the form of Eq. 6.4, and the powers can be evaluated on the bound- 
ary of the system. 

In general, the expressions for the "F contain the Fourier coef- 
ficients of the physically unmeasurable potentials and Lagrange mul- 
tipliers. Care must then be used in interpreting the Manley-Rowe 
formulas. The physical systems discussed in Sees. 6.1, 6.2, 6.3, 
and 6.4 obey Hamilton's principle, and it happens that the expres- 
sions for the 7^ do not contain the potentials that must be intro- 
duced [78]. One should not conclude, however, that this happens in 
all systems; the physical systems discussed in Sees. 6.6 and 6. 7 
have Manley-Rowe formulas including physically unmeasurable terms. 

6.6 Rotational Fluid Flow 

Here we discuss, in detail, a physical system with rotational flow 
of material, with or without electric charges. The analysis, which 
can be specialized to simpler systems merely by eliminating the ir- 
relevant terms, is for isentropic flow of a charged fluid, which is 
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compressible with a scalar pressure, under the action of arbitrary 
electric and magnetic fields and a gravity potential. 

This system is so complicated that the systematic procedure of 
Hamilton's principle actually simplifies the problem. In Sec. 6.6 
(Constitutive Relations) we discuss the properties of the system. The 
equations of motion are given in the second part of Sec. 6.6, and the 
power-flow theorem is derived. The Lagrangian is given in the third 
part of this section, and in Sec. 6.6 (Hamilton 1 s Principle) it is shown 
that the Euler-Lagrange equations coming from this Lagrangian are 
equivalent to the ordinary equations of motion., The Manley-Rowe 
formulas are given in the last portion of this section. 

Free charges and currents are not included in this system, but in 
Sec. 6.7 the modifications that are necessary are discussed. 

The Constitutive Relations. The electromagnetic variables are 
the electric field 35, the magnetic field IT, the electric displace- 
ment vector "D , and the magnetic induction vector B . At the out- 
set we recognize that "E and B" are related by two of Maxwell's 
equations 

V * B" = (6.133) 

and 

- (6.134) 



so we can define [30, Sec. 1.5] the usual vector and scalar potentials 
A and $, so that 

B = VXA (6.135) 

and 

"E = - V* -T (6.136) 

The fluid variables are chosen on a n per-particle" basis. It is as- 
sumed that each can be defined macroscopic ally as a function of space 
Y and time t. We let p be the number density of particles, IT the 
pressure (a scalar), 8 the temperature, s the entropy, T = p v the 
particle current density, "p the momentum per particle, and the grav- 
itational force per particle is -VQ, where & is the gravitational 
potential, assumed to be time -invariant. 

The momentum "p is, as in Sec. 6.3, the "canonical momentum, " 
which is not physically measurable; however, p - qA is physically 
measurable, where q is the charge per particle. We assume, as 
in Sec. 6.3, that "p - qA is a function of the velocity, so we can de- 
fine a new variable, the kinetic coenergy per particle, 

T 1 = C(p - qA) . dv" (6,137) 

We also define the "kinetic energy per particle" T by the Legendre 
transformation 
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T = (p - qA) v - T (6.138) 

In the theory of the motion of a single particle, T and T 1 play the 
role of energy state functions, but in our analysis they are considered 
as variables. We think of T 1 as a function jointly of the particle den- 
sity p and the particle current density T , rather than of the velocity 
v*. 

The functional dependence of p is, nonrelativistically, 

p = qA + mv (6.139) 

so 

1 2 

T = ^ (P - qA) (6.140) 

where m is the mass per particle. The relativistic dependence is 
more complicated: 



P = qA + = f r (6.141) 

^~^ 

thus 

I 
T 1 = - mc(c 2 - v 2 ) 2 (6.142) 

and 

1 

T = (m 2 c 4 -f (p - qA) 2 c 2 ) 2 (6.143) 

where m is the rest mass of the particle, and c is the velocity of 
light. 

The electrical constitutive relations most often used are those for 
free space, 

IT = e E (6.144) 

and 

B* = fxH" (6.145) 

where is the permittivity of free space, and u- Q is the permea- 
bility of free space. However, to allow nonlinear, dispersion-free, 
anisotropic, inhomogeneous, reciprocal, stationary electromagnetic 
media to occupy some parts of our system, we merely require (as 
in Sec. 6.1) that H be a reciprocal function of B" and .spacial co- 
ordinates r", and that 15 be a reciprocal function of "E and "r, so 
that we can define the energy state functions 

W (B , "r) = Magnetic energy density 



= f H(B , r) * dB (6.146) 
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W e ! (E, "f) = Electric coenergy density 

= VB(E, r) dE (6.147) 

and 

W e ("D, r) = Electric energy density 

= "E D - W e '(E, r) (6.148) 

The remaining constitutive relations are among the thermodynamic 
variables p, s, 6, and IT. We consider pressure and tempera- 
ture to be functions of p and s . On thermodynamic grounds we de- 
fine the "internal energy per particle" or "intrinsic energy per par- 
ticle" T by 



Equations of Motion. We choose to treat only isentropic flow, by 
which we mean that the entropy of a given particle remains constant. 
This restriction is expressed as the equation of isentropy, 

~~ = s + 7- Vs = (6.150) 

where D/Dt indicates the substantive derivative, taken in a refer- 
ence frame that at each instant moves with the particle. 
The equations of motion include the force equation, 

qT) = ? - qi:+ (^ V)p - q(v . ViX 

= Force per particle 

- - + q"E + qv X B" - Vfl (6.151) 

We note from Eq. 6.138 and the preceding discussion that (just as in 
Sec. 6.3) 

VT = (v * V)p - q(v * V)A + vX (VXp) - qvX (VX A) (6.152) 
From Eq. 6.149 it follows that 



) = + 0Vs (6.153) 

so that with the aid of Eq. 6.135 

V(T + T + - + fl) = -f + qS + qE + v X (V X p) + 0Vs (6.154) 

which is an alternate form of the force equation. This is one of the 
equations of motion. In addition, there is the equation of isentropy, 
here multiplied by p , 
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ps +1 * Vs = (6.155) 

We also have the equation of continuity, or particle conservation, 

v ' ^ = - P (6.156) 

and.Maxwell's equations, Eq. 6.133, 6.134, and 

v D = qp (6.157) 

VXH = qt +^ (6.158) 

We have assumed that the only charges and currents are those caused 
by the charged particles and their motion; that is, we allow no "free" 
charges or currents, such as appear in magnetohydrodynamic sys- 
tems. These are discussed in Sec. 6.7. 

The power -flow theorem follows from the equations of motion just 
given. Dot-multiply Eq. 6J.54 by I, multiply Eq. 6.156 by (T + T+ -+2 
dot-multiply Eq. 6.134 by H, and dot -multiply Eq. 6.158 by -. Then 
add the four resulting equations, combining the left-hand sides into the 
divergence of the power-flow vector IF 

The right-hand side is 

- Jf (W e + W m + pT + pi + pa) (6.160) 

The power-flow vector F" is a sum of terms representing flow of 
kinetic, intrinsic, and gravitational energy of the particles, and the 
term TTV, which represents a flow of power caused by motion of 
particles in overcoming a pressure force, and finally the electro- 
magnetic power-flow term. The expression for energy density, Eq. 
6.160, includes electric, magnetic, kinetic, intrinsic, and gravita- 
tional energy densities. 

We have now recognized the power-flow vector "F of Eq. 6.159; 
when the Manley-Rowe formulas are obtained, we will try to inter- 
pret the formulas for F in terms of Eq. 6.159. 

The Lagrangian. The Lagrangian density L is formed by taking 
the kinetic coenergy density, and subtracting from it the gravitational 
potential energy density and the intrinsic energy density. Added to 
this is the Lagrangian for the electromagnetic field, with coupling 
terms. The result is 

L = pT 1 - pT - p&+ W f - W m - qp$ + qT- "K (6.161) 

As it stands, the Lagrangian does not predict the correct equations 
of motion, because there are kinematic constraints. Two of the kin- 
ematic constraints, Eqs. 6.133 and 6.134, have already been used to 
define the potentials A and $ . Another kinematic constraint is the 
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continuity equation, Eq. 6.156, which we account for with a Lagrange 
multiplier, . Another kinematic constraint is the equation of isen- 
tropy, Eq. 6.155, for which we define the Lagrange multiplier (3. 
If these were the only kinematic constraints used, the resulting 
equations of motion would not be general, but instead would predict 
only a limited class of flows. One of the resulting Euler-Lagrange 
equations would be 

V = - p + |3Vs (6 .162) 

and so, for uniform entropy, only irrotational flow would be pre- 
dicted. 

Early attempts to predict more general motion are those of Ito 
[45], Ziman[l21], and Thellung [105] . These authors employed a 
Clebsch transformation [51, pp. 248-249] expressing the momentum 
p in the form 

p = - V<f - XV4/ (6.163) 

Neglecting for the moment the electrical and thermodynamic vari- 
ables, the Lagrangian is taken to be 

p( + X+ - T - T - n) (6.164) 

and is therefore a function of p, $ , X, and i|j, and their deriva- 
tives. The Euler-Lagrange equations are the force equation, the 
equation of continuity, and the two equations 

|f = (6.165) 

and 

f = (6.166) 

This derivation, using the Clebsch transformation, has been extended 
to include electromagnetic forces [67] and relativistic motion [112] . It 
does not require the introduction of the continuity equation as a kine- 
matic constraint, but the form of the Lagrangian and the potentials 
are not easily given a physical interpretation. We prefer to use an 
alternate method, which is similar analytically, but different in phil- 
osophy. 

This method, due to Lin [54, 92], shows clearly the particle-like 
basis of fluid mechanics. At some time t Q we give to each particle 
of our system a label which is its position at that time, "a . The par- 
ticle keeps this label throughout its travel, and consequently we can 
define a vector field "a as a function of space 7 and time t. Opera- 
tionally, we determine a" corresponding to a given point Y at a given 
time t by looking at that point, at that time, and noting what particle 
is ttuere. The value of a" associated with that particle then is defined 
as (r, t) . 
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Stated another way, we go to a point T and at time t we ask the 
particle passing by, "Where were you at time t ?" The particle 1 s 
answer defines the value of the vector field "a at "r and t. 

The equation expressing the fact that each particle keeps the same 
value of "a is 

^| = i + (v o V)a = (6.167) 

or, when multiplied by p , 

pt + (I . V)a" = (6.168) 

We use thisjts a kinematic constraint, to be handled with a Lagrange 
multiplier X 

The Lagrangian of Eq. 6.161 is then replaced by 

L = pT 1 - pT - ?Q + W e ' - W m - qp$ + qT - A 

- Up + V I) - p(ps +T - Vs) 

-T - [pt -f (t * V)a] (6.169) 

which is considered a function of the variables p , T, s , $, "a, t, , 
p,T, and their space and time derivatives, as well as space ex- 
plicitly. 

Hamilton 1 s Principle. Hamilton 1 s principle, Eq. 6.122, predicts 
Euler -Lagrange equations in the form of Eq. 6.124. We omit the de- 
tails and write down the Euler-Lagrange equations. There is one 
such equation associated with each variable. 

Corresponding to the variation in p we find 

I = T + T+- + a + qS+7C-i+ps (6.170) 

the equation obtained by varying T is 

VI = - p + S.X.Va. + (3Vs (6.171) 

where, as in the following, Z. indicates a sum over the three di- 
rections in space; the variation in s yields 

P 6 - (p?) - V (Ip) = (6.172) 

which, with the aid of Eq. 6.178, can be put into the form 

D6 

^-=(3 + v-Vp = 6 (6.173) 

the formula resulting from the variable A is 

VXH = qT +15 (6.174) 

associated with the potential < is the equation 
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V - TJ a qp (6.175) 

the formula from the variation in a" is 

JL (pX.) + V (TX.) = (6.176) 

ot j j 

or, in view of Eq. 6.178, 

Iji = I + (v * V)X = (6.177) 

corresponding to the variations in the Lagrange multipliers, , (3, 
and "x we obtain merely the three kinematic constraints, 

p + V - T =0 (6.178) 

ps +1 . Vs = (6.179) 
and 

pi + (T V)a = (6.180) 

It is not obvious that these Euler -Lagrange equations are identical 
to the ordinary equations of motion. To show that they are, we must 
derive the ordinary equations of motion from the Euler -Lagrange 
equations, and vice versa. 

It is convenient to introduce a lemma to help do this. Let us sup- 
pose that "a and X are any two vectors whose substantive deriva- 
tives vanish: 

t+ (v . V)a = (6.181) 

and 

T + (v V)T = (6.182) 

also let us suppose that s, (3 , and 6 are scalar functions of space 
and time such that 

s 4- v- Vs = (6.183) 

and 

P + V - V|3 = 6 (6.184) 

Then it is identically true that 

v X [v X (S.X.Va.) + (3Vs)l + 0Vs 

= - (S.X.Va. + pVs) - V(X t + (3i) (6.185) 

* J J J 

We now show that the ordinary equations of motion can be derived 
from the Euler -Lagrange equations. The only problem is in deriving 
the force equation, Eq. 6.154. We recognize the left-hand side of 
Eq. 6.185 as 
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vX (VXp) + evs (6.186) 

by use of the Euler-Lagrange Eq. 6.171. We now take the gradient 
of Eq. 6.170, and subtract from it the time -derivative of Eq. 6.171, 
to eliminate the Lagrange multiplier . The result is 






- V(X f + (3s) - ? + -- (S.X.Va. + j3Vs ) 
P dt J J J 

(6.187) 
which, using Eq. 6.136 and the lemma just proved, reduces to 

V(T + T + ~ + n ) as - p + qA + qE + v X (V X p) + 0Vs (6.188) 

which is identical to the force equation, Eq. 6.154. We have thus 
shown that Hamilton's principle predicts the correct equations of 
motion. 

The fact that the Euler-Lagrange equations predict the correct 
equations of motion does not imply that the converse is true that 
is, that the Euler-Lagrange equations are general enough to predict 
all possible types of motion. This is shown by deriving the form of 
the Euler-Lagrange equations from the known equations of motion. 
The earliest formulations of the fluid-flow problem in Eulerian co- 
ordinates failed to predict the most general type of motion. Lin's 
method, which we have used here, does predict the most general 
motion; this was first shown by Serrin, in a letter [91] to Lin. We 
use here a method somewhat different from Serrin 1 s. 

The only Euler-Lagrange equations that we need verify are Eqs. 
6.170, 6.171, 6.173, 6.177, and 6.180, since the others are identical 
to the equations of motion. Further, Eq. 6.180 is automatically sat- 
isfied if we define the vector field "a(r , t) f as before, as the initial 
position vector of the particle at "r at time t. 

Equation 6.173 can be satisfied by defining p properly. Let p be 
the time-integral of the temperature of a given particle, starting at 
time t Q . Thus Eq. 6.173 holds, and p(r,J; ) =0. 

Similarly, we can define a vector field X so that Eq. 6.177 is sat- 
isfied. At time t Q , let T be any vector with a curl identical to the 
curl of the momentum p. For example, X could be chosen equal 
to the momentum at time t Q . Now, at each point we assign the value 
of "X to the particle there, and let the particle carry this value as a 
label, just as it does the value of "a. From this definition it is clear 
that the substantive derivative of "X vanishes, and Eq. 6.177 holds. 

Now we need only prove Eqs. 6.170 and 6.171. We use the lemma 
just proved above. We have chosen "a, X, p, and 6 so that Eqs. 
6.181, 6.182, 6.183, and 6.184 hold, so the result of Eq. 6.185 is true. 
We define a vector field u(r", t) by 

u = V X (p - S.X.Va. - pVs) (6.189) 
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noting that at time t y u vanishes everywhere. From the force 
equation, Eq. 6.154, 

V x ^ = V X [v x (V X p)] + V x (9Vs) (6.190) 

from Eq. 6.185 we can evaluate V X (v Xu) as 

VX (vXu) = VX [vX(VXf)] + VX(0Vs) --- [VX (S.X.Va. + (3Vs)l 

01 L J J J J 

(6.191) 

thus 

"u" = V X (v X u) 

= (u - V) v - (? V) u - u(V - v) (6.192) 

Equation 6.192 then shows that 

~i = {u. V)?-u(V* v) (6.193) 

With specified flow, this is a set of three coupled linear time -vary ing 
total differential equations for the three components of "u; the solu- 
tion is thus unique. A possible solution is the trivial one u = , and 
because u = everywhere at time t Q , we see that this is therefore 
the solution. Thus 

u = V X (p - S.X.Va. - pVs) = (6.194) 

J J J 

so that for some scalar ~u(r , t) , 

SA.Va. + (3Vs - p = Vv (6.195) 

In view of Eqs. 6.194 and 6.195 the lemma of Eq. 6.185 reads 

v x (V x p) + 9Vs = f + Vi - V(5C t + ps) (6.196) 

so that the force equation, Eq. 6.154, can be altered to the form 

V(T + T + - + Q + q*) = Vv - V(T - t + PS) (6.197) 

P 

But two scalar s with equal gradients can differ only by a function of 
time, so for some v(t) , 

T+T+-+Q-fq$-!-T*a + ps = u + v (6.198) 

We now define the scalar function of space and time (r", t) as 

r * 

;(r, t) = \ v(t)dt + u(r, t) (6.199) 

*o 
so that 
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= T + T+~ + a-fq$ + X-i-f(3s (6.200) 

and 

Vt, = Vi; 

= - p + S.X.Va. + pVs (6.201) 

j j j 

Because these are exactly the Euler-Lagrange Eqs. 6.170 and 6.17! y 
we have shown that it is possible to derive the Euler-Lagrange equa- 
tions from the equations of motion. Thus Hamilton 1 s principle does 
in fact predict the most general kind of motion. 

^ We now have a physical interpretation for the Lagrange multiplier 
X . If the flow is irrotational at time t Q , we can take T = ; if not 
we can consider \ to be the original momentum. 

We have just shown that Hamilton's principle leads to Euler-La- 
grange equations that are identical in meaning to the ordinary equa- 
tions of motion. The ordinary equations of motion are nonlinear and 
contain only physically measurable variables; the Euler-Lagrange 
equations, on the other hand, are linear (if one considers the quan- 
tities that appear as variables) and contain terms with physically un- 
measurable potentials and Lagrange multipliers. The Euler-Lagrange 
equations lead to Manley-Rowe formulas expressed so that the powers 
can be calculated on the boundary of the system. 

Manley-Rowe Formulas* The Manley-Rowe formulas have not 
previously been given for this system, except for the case of uni- 
form entropy and irrotational flow. The reason the general case 
could not be done earlier is that the physically unmeasurable quan- 
tities X , "a , and (3 appear in the formulas; without the motivation 
of Hamilton's principle we would not be led to define these quantities. 

The systematic procedure outlined in Sec. 6. 5 yields Manley- 
Rowe formulas in the form of Eq. 6. 2 with the P a in the form of 
divergences: 



P a = -V T a (6.Z02) 



where 



F = iRe/E XH*+T fl * + X*+(l) *1 
a 2 \ a a a I a a v P 'a J 

+ S.[l (X.a.) * - (IX.) a *1 

J L a- J J a J a ja J 

+ I(p5) *- (IP) s *\ (6.203) 

a a a a J 

The Manley-Rowe formulas can now be integrated over space, and 
put into the form of Eq, 6.4. Although the powers can be evaluated 
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at the boundary of the system, the expressions for F are not what 
one would expect from a knowledge of the power -flow vector 7 of 
Eq. 6.159. The physically unmeasurable terms, which disappear for 
fluid with uniform entropy and irrotational flow, cannot be given a 
satisfying physical interpretation. 

It should be stressed that in deriving the frequency-power form- 
ulas, it was assumed that each of the variables could be written in 
a Fourier series of the form 



p(r, t) = p (r) +Re S^Cr) e (6.204) 

Whereas there may be physical reasoning to support this assumption 
for the physically measurable variables, no such reasoning applies 
to the physically unmeasurable variables, including a~, "X , and (es- 
pecially) j3. It may be true that in most practical cases the flow is 
such that the formulas hold, but one should make sure of that. 

This treatment can be specialized to systems of practical inter- 
est simply by omitting the irrelevant terms. For isentropic fluid 
flow, including acoustics, eliminate the electrical variables by set- 
ting "E, B*, U, H, q, ~K, and $ all to zero. If the entropy is 
uniform, drop all terms that involve s , 9 , and (3 . If the flow is 
irrotational, set V X "p to zero, and eliminate all terms involving 
"a or X . On the other hand, if the fluid is incompressible, make p 
a constant; if, in addition, the entropy is uniform, all terms includ- 
ing T can be neglected. 

For electron clouds, including electron beams, the effects of en- 
tropy, gravity, and pressure are usually ignored; therefore, elim- 
inate all terms involving IT, s, 6, T, (3, and ft . For irrotational 
flow, ignore the terms containing "a or "x, and set V X "p = . 

6. 7 Magnetohydrodynamic Fluid 

Here we discuss a magnetohydrodynamic fluid, a fluid with free 
charges and currents [17 j. We assume the conductivity is infinite, 
so that there are no losses. This assumption forces the electric 
field in the reference frame moving along with the fluid to vanish. 
Because this is a significant qualitative difference from the system 
discussed in Sec. 6.6, one cannot obtain the equations of motion of 
that system in the limit as the free charges and currents are van- 
i shingly small. 

For generality, we include in the analysis all the physical effects 
discussed in Sec. 6.6, and, in addition, free charges and currents. 
Thus we allow the fluid to be charged, so that the total electric charge 
and current is composed of parts caused by the fluid charge and parts 
caused by the free charge. If the fluid charge and current can be neg- 
lected compared with the free charge and current, simply put the 
charge per particle, q, zero. If the entropy is uniform, drop all 
terms in the analysis that involve the entropy s , the temperature 
9, and the Lagrange multiplier p. If the fluid is assumed to be in- 
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compressible, simply treat the particle density p as a constant; i 
in addition, the entropy is uniform, all terms involving the intrinsic 
energy per particle 1 may be neglected, since X is constant. 

Instead of repeating much of the description of Sec. 6.6, we merely 
give the modifications that are necessary for that discussion to ap- 
ply to systems with free charges and currents and infinite conduc- 
tivity. Generally speaking, this consists of adding terms to exist- 
ing equations and adding a few new equations and some new discus- 
sion. 

Constitutive Relations. In addition to the variables defined in Sec. 
6,6 (Constitutive Relations), we have 

p = Free charge density (6.205) 

"j = Free current density (6.206) 

Therefore, the total electric charge density is 

qp + P f (6.207) 

and the total electric current density is 

qT + J f (6.208) 

The constitutive relations discussed in Sec. 6.6 (Constitutive Re- 
lations) are unchanged; a new constitutive relation is Ohm's law for 
the fluid. We assume the fluid is infinitely conducting, so that the 
electric field within the fluid vanishes. Expressed in terms of the 
electric and magnetic fields measured in the stationary frame of ref- 
erence, 

E; + v X B" = (6.209) 

Equations of Motion. The equations of motion given in Sec. 6.6 
must be altered. The forces on the fluid caused by free charge and 
current modify the force equation, Eq. 6.154, to 

V(T + T -f ~ + fl) = - p + q 4- qE + v X (V X p) 

p E J X 1" 

+ 0Vs + -4- + (6.210) 

p p 

The equation of isentropy, Eq. 6.155, and the equation of particle 
conservation, Eq. 6.156, remain unchanged. The remainder of the 
equations of motion are Maxwell's equations; Eq. 6.133 and 6.134 
are unchanged, but the other two are modified to include the new 
charge density and current density: 

V - D = qp + p (6.211) 
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VXH=qi+7 f + TF (6.212) 

In addition to these equations of motion is the equation of conserva- 
tion of free charge, 

V* J .+ p f = (6.213) 

Also, Ohm f s law, Eq. 6.209, appears as one of the equations of mo- 
tion. 

The derivation of the power-flow theorem 

~(W e + W m + pT + pT + pQ) + V * F = (6.213 

of Sec. 6.6 remains valid even in the presence of free charges and 
currents, when Eq. 6.209 is taken into account. The physical power- 
flow vector "F is the same as the power-flow vector found in Sec, 
6,6, and is given by Eq. 6.159* 

The Lagrangian^ The Lagrangian of Eq. 6.161 must be modified 
to include the coupling terms of the free charges and currents. It 
becomes 

L = pT 1 - pT - p2 + W e f - W m - qp$ + qT - A - p $ + J^ A (6.215) 

Besides the kinematic constraints discussed in Sec. 6.6, we con- 
sider Eq. 6.213, the free charge continuity equation, as a kinematic 
constraint. We choose to take account of it by defining a suitable 
potential. One possible choice for a potential might be the "pseudo- 
polarization" vector R of Panofsky and Phillips [7 6, Sec. 13-4] , with 
the property that 

J f = F (6.216) 

"and 

p = - V - F (6.217) 

We call this a pseudo-polarization vector because of the similarity 
to the polarization vector 3? of a dielectric medium. The polari- 
zation charge density pp aSd the dielectric current density J are 
related to the polarization vector P" by [23, Sec. 5. 2] p 

p -V- P (6.218) 

P P P 

and 

J = F (6.219) 

P P 

if the dielectric medium does not move. 

However, the use of the vector R as a potential is not satisfac- 
tory when fluid motion is present. We can obtain a hint about the 
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proper potential by noting that in the case of a moving dielectric, the 
polarization current is no longer given by Eq w 6.219, but rather by 
[23, Sec. 9. 2] 

J = F + V X (F X v) (6.220) 

We define a pseudo-polarization vector ll by analogy with the dielec- 
tric case, by requiring that 

p = - V - 5 (6.221) 



and 



j = n + V X (IIXv) (6.222) 



This representation for the free current and charge was suggested 
by Haus[38]. 

It is not clear as yet that arbitrary currents can be represented 
in the form of Eq. 6.222. To realize this, consider the differential 
e quati on 

5 = (H V)v 1(V . v) + ^ - vp (6.223) 

For a given particle, for a given flow, this is a set of three coupled* 
time-varying, total, differential equations. A unique solution ex- 
ists with given initial conditions 5(t Q ) specified. Putting together 
the solution for each particle, we define n(r , t) as a function of 
space and time, uniquely, for a given initial condition II(r", t ) . We 
select the initial condition in which 

V. E(r, t ) = -p f (r, t ) (6.224) 

Now Eq. 6.223 can be written in the form 

5+ VX (5x v) = 7 - vp - v(V S) (6.225) 

of which the divergence is 

-=(P + v - n) = (p + v - n) (v * v) (6.226) 

where Eq. 6.213 has been used to eliminate ~3 . But Eq. 6.226 is, 
for a given particle and given flow, a time -vary ing, total, differen- 
tial equation, and has a unique solution. By inspection, a possible 
solution is 

p + V - 5 = (6.227) 

In view of the initial condition of Eq. 6.224, this is the unique solu- 
tion. Thus we have defined !T(T, t) in such a way that 

p f = - V - n ( 
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and, from Eq 6.225, 

J = IT + V X (5 X v) (6.229) 

Note that there is no restriction on T and p other than the contin- 
uity requirement of Eq. 6.213. 

"When the kinematic constraints are accounted for, the Lagrangian 
of Eq. 6.215 becomes, instead of Eq. 6.169, 

L = pT r - pT - pQ + W e ' - W m - qp$ + qt A 

+ $(V - n) + A * 5 + A vx(5x r/p) 

-(P + V* T) - p(pi+T. Vs) 

-I * [pi + (I - V)a] (6.230) 

which is considered a function of the variables p , I, s , A, $ , aT, 
EFy > P> X , and their space and time derivatives, as well as space 
explicitly. 

Hamilton 1 s Principle. Some of the Euler-Lagrange equations of 
Sec. 6. 6 (Hamilton's Principle) are changed, and others are not. Un- 
changed are Eqs. 6.173, 6.177, 6.178, 6.179, and 6.180. Equations 
6.170, 6.171, 6.174, and 6.175 are changed to the formulas: 

; =T + T+~ + a + q $ + X.i + |3s + 5 ' ^ XB (6.231) 

V = - p + S.X .Va. + (3Vs + -3-^-2 (6.232) 

J J J r 

VXH = qI+T f +"C (6.233) 

and 

V - "5 = qp + p (6.234) 

Also, there is^a new Euler-Lagrange equation, resulting from the 
variation of II . It is 

"E + vXB" = (6.235) 

We now have the problem of showing that the Euler-Lagrange equa- 
tions are equivalent to the ordinary equations of motion. It is easily 
verified that all the ordinary equations of motion except the force 
equation follow from the Euler-Lagrange equations. Furthermore, 
the only Euler-Lagrange equations that require verification are Eqs. 
6.231, 6.232, 6.173, 6.177, and 6.180. 

It is convenient to introduce a lemma similar to Eq. 6.185. If 

J f =S-hVX(lx7) (6.236) 

P f = - V - ff (6.237) 
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VX "E = - "5 (6.238) 

"E + VXB" = (6.239) 

p + V* (pv) = (6.240) 
and 

V - "B = (6.241) 

then it is identically true that 

7 X if pi 

~- +4" ( 6 -242) 



The force equation can be derived easily now. Take the gradient 
of Eq. 6.231, and subtract from it the time -derivative of Eq. 6.232, 
thus eliminating the Lagrange multiplier . Use the two lemmas, 
Eq. 6.185 and Eq. 6.242, and calculate the quantity 

vX(VXp) (6.243) 

from Eq. 6.232, to obtain the force equation in the form of Eq. 6.210. 
We thus conclude that Hamilton's principle predicts the correct equa- 
tions of motion. 

The generality of Hamilton's principle is shown by predicting the 
form of the Euler-Lagrange equations from the known equations of 
motion. We define "afr, t) and (3(r , t) as in Sec. 6.6. The defi- 
nition of X (r, t) is altered only in its initial value; we put*X(F, t ) 
equal to any vector so that at time t 



VXX = vxp-. (6.244) 

We have derived all the Euler-Lagrange equations except Eqs. 6.231 
and 6.232. Note that the variables "a, X, s, (3, and 8 have been 
chosen so that the lemma of Eq. 6.185 holds, and note that Eqs. 6.236, 
6.237, 6.238, 6.239, 6.240, and 6.241 are predicted from the equa- 
tions of motion and the definition of the pseudo -polarization U, so 
that Eq. 6.242 is known to be valid. 

Now let us define a vector field ufr , t) by an extension of the def- 
inition of Eq. 6.189: 



p - S.X.Va. - pVs - ~p. j (6.245) 

noting that at time t Q , "u vanishes everywhere. As redefined here, 
"u obeys Eq. 6.192. 

The reasoning given in Sec. 6.6 is only slightly altered to fit the 
new function ufr", t) . We find from the differential equation, Eq. 
6.193, that u must vanish, so that for some scalar v( , t), 

S.X.Va. + pVs +5211. . p = vv (6.246) 

J J J P 
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Using the two lemmas, Eqs. 6.185 and 6,242, and the force equation, 
Eq. 6.210, this result can be put into the form 



v (T 



+ T + + a + q ^ ) = vS - v(\ i + pi) - v * (6.247) 



From this point, the same reasoning used in Sec. 6.6, and the 
same definition for ,(r", t), lead to the Euler-Lagrange equations, 
Eqs. 6 231 and 6.232. Thus we have shown that the most general 
type of motion is predicted by Hamilton 1 s principle. 

Manley-Rowe Formulas. The simplest form of the Manley-Rowe 
formulas is obtained not when L is used as the state function, but 
rather when 

) (6.248) 



is used. This differs from L only by a perfect divergence; it is 
known that both state functions can be used in Hamilton's principle, 
and lead to the same Euler-Lagrange equations [70, Section 3. 4] . The 
Manley-Rowe formulas generated by this state function are of the 
form of Eq. 6.4, with the vector F given by 

F n = ReE XH*+T ** * 

2 



>a 



a. 1 
J 



. * 
a ja 



Observe that this differs from the corresponding result of Sec. 6.6, 
Eq. 6.203, only by the last two terms. The physical interpretation 
of these terms is not clear. 
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APPLICATIONS 



Chapter 7 
ROTATING MACHINE APPLICATIONS 



It is shown in Sec. 4.6 that some rotating machines (those with- 
out commutators) obey the Manley-Rowe formulas if losses within 
the machine are neglected. Here we demonstrate how these fre- 
quency-power formulas can, in simple cases, lead directly to re- 
sults of practical importance. In Sec. 7.1 we discuss the conven- 
tional two-phase induction machine with balanced excitation, and 
show that the Manley-Rowe formulas predict the existence of three 
regions of operation a motor region, a brake region, and a gen- 
erator region. 

In Sec. 7.2 the balanced treatment is extended to cover the unbal- 
anced machine, or the machine run with unbalanced excitation on 
the stator, but a balanced rotor. There are now four frequencies 
of interest, so we require an additional constraint; this is supplied 
by a simple, yet general, rotor model. As special cases of the gen- 
eral solution we discuss the nearly balanced machine and the com- 
pletely unbalanced machine. 

The d-c induction machine is discussed in Sec. 7.3, and the gen- 
eral expressions obtained from the Manley-Rowe formulas are also 
shown from the actual solution. 

7.1 The Balanced Induction Machine 

As shown in Sec. 4.6, the induction machine [116, Sec. 3.6.4] obeys 
the Manley-Rowe formulas, provided stator and rotor loss and 
mechanical friction are taken out. Of course, rotor loss is essen- 
tial if the machine is to operate, but for the purpose of applying the 
Manley-Rowe formulas this resistance must be considered as "out- 
side" the system. This requirement is purely conceptual, however, 
and the analysis is in no way restricted to machines with a wound 
rotor, in which this separation can (in part) be carried out physi- 
cally. 

We excite the stator with balanced currents at frequency co . Be- 
cause the stator and rotor are both assumed balanced, the only cur- 
rent in the rotor windings is at frequency 

_ = .- m (7. 1) 

where w is the (assumed constant) mechanical speed. 

If we call P_ the power input to the rotor (that is, the negative 
of the power dissipated by the rotor resistance), P g the power in- 
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pat to the stator, and P the mechanical power input, and if we 
consider GO and to to be the two independent frequencies, the 
Manley-Rowe formulas are 



= 



(7,2) 



and 



n 

+ =0 



(7.3) 



Without loss of generality, we assume that <o is positive; however, 
either co or 00 can "be negative. 

In the normal operation of the induction machine we do not put 
power into the rotor, so P_ < . For operation as a motor, we put 
power into the stator and get power out the shaft, so that P g > 
and P < . For operation as a generator, on the other hand, we 
put power in the shaft and get power out the stator, so P g < and 
P > . And finally, there is a third type of operation, known as 
"brake" operation, in which power flows in both the stator and the 



shaft, so P g > and 
in Table 7.1. 



m 



These conditions are summarized 





P_ 


P s 


P 

m 


MOTOR 


Neg. 


Pos. 


Neg. 


GENERATOR 


Neg. 


Neg. 


Pos. 


BRAKE 


Neg. 


Pos, 


Pos. 



Table 7.1. Balanced machine power flow. Each power 
shown is the power into the machine 

It is clear from Eqs* 7.2 and 7.3 that motor action results only 
when both GO and co_ are positive, or when 

< m < tt s (7-4) 

On the other hand, we have generator action only if a> m is positive 
and oj_ is negative, that is, 

"s < "m ( 7 * 5 ) 

When oa m is negative we must have a positive <o_ , and Eqs. 7, 2 
and 7.3 predict that both P and P are positive, so we get brake 
operation. 

A common dimensionless parameter that expresses the shaft speed 
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co is the slip s : 

m * 
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For negative slip we have a generator, for < s < 1 we have a 
motor, and for s greater than 1 a fcrake is obtained. These results 
are conveniently summarized in Fig. 7.1, which shows, as a func- 
tion of co m (or of s ), the three regions of operation. 
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Fig. 7.1, The three regions of operation of a balanced 
induction machine the brake region (B) , 
the motor region (M) , and the generator 
region (G) . The Greek letters show the cor- 
respondence to the more general case of Fig. 
7.2 

In the motor region, we define the motor efficiency rj m as the 
ratio of mechanical power out to stator power in, so that by Eqs, 
7.2 and 7.3 



-P 



(7.7) 



Furthermore, the starting torque T is expressed as 



"m * s->l 

where the notation s -* 



(7.8) 



s 



indicates the limit as 



ing torque figure of merit sometimes used is f Q 
the starting torque T times synchronous speed 



o ~+ , A start- 

m. 
the negative of 

co , divided by 

the stator power input P required to sustain this torque. In view 
of Eq. 7.8, this figure of merit for the balanced motor is unity, so 
for the unbalanced cases discussed in Sec. 7.2, f will express, in 
a sense, the "starting efficiency" of the motor, as compared to an 
"equivalent" balanced motor. 

The "free-running speed" co,. is that speed for which the torque 
vanishes, or the speed the motor approaches if unloaded. This 
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speed is easily found by supposing that the torque is a continuous 
function of shaft speed. We notice that the applied torque must be 
negative for motor action and positive for generator action, so the 
torque must vanish when the slip is zero, or when w m = co g . This 
point is identified on Fig. 7.1 as w^. 

In the generator region, we define the generator efficiency 77 g as 
the ratio of electrical power out the stator to mechanical power in, 



" P s w s 

j - ._ - 

P m "m 



Although defined for different regions of operation, the generator ef- 
ficiency 77 is the reciprocal of the motor efficiency r] m . In terms 
of r\ we may classify regions of operation as motor regions 
(0 < 7^ < 1) , generator regions (r] m > 1) , and brake regions (r7 m < 0) , 

These facts about balanced induction machines are well known, and 
are derived here in this way partly as an indication of the method we 
use in treating the unbalanced machine in Sec. 7.2. 

Use has already been made [86, 85, 21] of the formulas of Eqs. 7.2 
and 7.3 for balanced induction generators. References 86, 85, and 
21 deal with configurations of coupled balanced induction generators, 
and show for some particular cases that the efficiency constraints of 
Eq. 7.9 hold for the overall systems, as well as for each individual 
machine. 

There is some need for a variable- speed, constant-frequency gen- 
erator for aircraft use, and configurations of induction generators 
have been proposed in the past to do this at high efficiency, in viola- 
tion of the Manley-Rowe formulas. Patents have even been issued 
[86,85,21] for devices claiming to do this, although such machines 
have never been found to work. 

It is interesting to note that our general formulation indicates how 
to get around this efficiency limitation merely use devices that do 
not obey the Manley-Rowe formulas. The two simplest such devices 
are the switch and the rectifier. Efficiency can be improved by using 
commutator machines, or by rectifying and recovering the otherwise- 
lost rotor power. It cannot be improved, however, by using networks 
of balanced induction machines. 

7.2 The Unbalanced Induction Machine 

The general analysis of Sec. 7.1 is extended here to unbalanced in- 
duction machines [116, Sec. 4.6.3; 25, Art. 10.8] . The same physical 
device is used for both the balanced and unbalanced induction ma- 
chines; the difference is in the stator excitation. We no longer re- 
quire a strictly balanced stator input, although we will insist on a 
balanced rotor, to limit the number of frequencies present. As in 
Sec. 7.1, we neglect stator winding resistance, mechanical friction, 
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and other losses, and place the rotor dissipation "outside" the sys- 
tem. Again, however, this requirement is purely theoretical, and 
the analysis holds for squirrel-cage motors, in which this cannot be 
done physically, as well as for motors with wound rotors. 

In unbalanced operation the shaft turns at a (constant) speed o) , 
and the stator is excited with current of frequency co . The current 
in the rotor is a sum of currents at frequency 

o> + = w g + w m (7.10) 

and at frequency 



Let us call P. and P the power input to the rotor (that is, the 
negative of the power dissipated by the rotor resistance) caused by 
currents in the rotor at frequency co and frequency co_ respect- 
ively. The stator power input, which (because the rotor is assumed 
to be symmetrical) is entirely at frequency co g , is P g , and the 
mechanical power input, P m , is the torque r multiplied by the 
speed co m . 

In an actual machine, the stator power input differs from P only 
because of losses we have neglected, primarily stator resistance; 
and the mechanical power input differs from P m only because of 
mechanical losses. 

Because there is power flow only at the four frequencies co , co m , 
co, , and co , the frequency-power formulas become 

P P, p 

IT + ^r~if" = (7J2) 

and 

P P, P 

-1 + + = (7.13) 

co g co + <*_ 

Just as in the balanced induction machine, no power is put into the 
rotor. Thus P, < and P_ < . For operation as a motor, we 
put power in the stator and get power out the shaft; as a generator^ 
power goes in the shaft and out the stator terminals; and as a brake, 
power flows in both the stator and the shaft. 

We put Eqs. 7.12 and 7.13 in the form 

P aw f^.-t 

m m \ co co 
\ - + 

and 

P s = -s(^ + ^) (7 ' 15) 
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The motor efficiency r) m is defined just as in Sec. 7.1. By the 
preceeding equations it is 



-P 



(7.16) 



This efficiency is used to determine the regions of motor action 
( < 77 < 1 ) generator action (' rj m > 1 ), and brake action ( r? 
< ). Table 7,2 summarizes these conditions. 





P + 


P. 


P s 


P m 


^m 


MOTOR 


Neg. 


Neg. 


Pos. 


Neg. 


o<n m <i 


GENERATOR 


Neg. 


Neg. 


Neg. 


Pos. 


^m> 1 


BRAKE 


Neg. 


Neg. 


Pos. 


Pos. 


^m <0 



Table 7.2. Unbalanced induction machine power flow 

< g or < s < 2, then 



We notice immediately that if -co g < 
both w and w , are positive, so that the quantity 



_ , 



(7.17) 



is between -1 and +1, so that the 



is less than unity. Thus 77 m 

machine is either a motor or a brake. Thus when - co < co < o> , 

o II JL o 

the machine cannot act as a generator. 

Similarly, if w does not lie between -co and co g , either o) + 
or o> (but not both) will be negative, so that the quantity of Eq. 7.17 
is greater than unity. Thus r7 m is either larger than one or less 
than minus one, so that the machine cannot operate as a motor. It 
must be either a generator or a brake in these regions. 

In words, this means that for shaft speeds less than synchronous 
speed (in either direction) the device cannot operate as a generator, 
and for shaft speeds in either direction greater than synchronous 
speed it cannot act as a motor. 

The torque r is given by 



P co, - P, co 

"* T T "" 

P CO + P CO 



so that points of zero torque occur when P_& 



(7.18) 



The start- 



ing torque T O is merely the expression of Eq. 7.18 in the limit as 
the slip s-*l . The starting torque figure of merit defined in Sec. 
7.1 then becomes 



f o = 



(7.19) 
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The analysis thus far is general for unbalanced induction machines. 
To simplify the equations we examine the ratio of P^OJ^ to Pj.w_ . 
We adopt a rotor equivalent circuit, which is quite general. In ad- 
dition, we define an "unbalance parameter" , which describes how 
close to balanced or unbalanced operation the machine is. 

The Rotor Model. Let us break up the magnetic field in the air 
gap caused by stator currents into a component rotating in the neg- 
ative direction, and a component rotating in the positive direction* 
The former generates rotor current at frequency to, , and the latter 
generates rotor current at frequency co_ . Let t, be the square of 
the ratio of the magnitude of these rotating field components. Thus 
for balanced operation, there is only the field rotating in the posi- 
tive direction, so = . On the other hand, for completely unbal- 
anced operation, that is, either with only one stator winding excited, 
or with only single -phase current in the stators, we have = 1 . 
For balanced operation in the other direction, = + oo. 

As an example of the calculation of in a practical case, sup- 
pose that a two-phase servomotor has one stator excited with cur- 
rent I a cos o g t and the other excited with current 1^ sin w g t . It 
is easily seen that 

12 

(7.20) 

The induced voltage on the rotors is proportional to the time de- 
rivative of the field, so using this decomposition of the gap field 
into rotating components, the ratio (squared) of induced voltages 
at the two frequencies, E + at oo^, and E^ at c^ , becomes 

(7.21) 

As a rotor model we assume an inductance L in series with a 
resistance R . Since only their ratio is important, they may be de- 
fined on a per-phase basis, if desired. The impedances (squared) 
presented to the rotor current at the sum and difference frequencies 
are, respectively, 

Z + z = R 2 + c* + 2 L 2 (7-22) 

and 

Z_ 2 = R 2 + o>_ 2 L 2 (7.23) 

from which we may determine the ratio of rotor currents (squared), 
and thus the ratio of P + to P_ , as 

22 7 - Z 2 nr 2 

P+ E + Z_ R 
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so that 

_LJ1 , t 2_l f * a * (7.25) 

F_+ s X 2 + {2 - S ) 2 

where we define 

x - ~- (7.26) 

0) JLi 

independent of the speed (or slip) of the machine. 

Using this rotor model we express the efficiency, torque, etc., in 
terms of the three parameters x, s, and . These may or may 
not be independent for any given excitation. The rotor inductance 
and resistance and the excitation frequency determine x , and so 
in many cases of practical interest this is a constant. The remain- 
ing two parameters, s and * are not > * n general, independent. 
They are (i) for the balanced machine ( t, = ), (ii) for the com- 
pletely unbalanced machine ( = 1 ), and (iii) when the stators are 
driven from current sources. For other types of operation, as for 
example when one stator winding is excited through a phase -shifting 
network, or where voltage sources are used, will be found to de- 
pend on s . Thus any general examination should take into account 
the (possibly rather complicated) dependence of on s . 

The General Solution. In terms, then, of the parameters x, t> , 
and s, we express the motor efficiency as 

s [x 2 + (2 - s) 2 ] - t>(2 - s)(x 2 + s 2 ) 
^ = (1 " S) s[x 2 + (2- s) 2 ] +t(2- s)(x z + s 2 ) (7 " 2?) 

the torque as 



"s s[x 2 + (2 - s) 2 ] + ;(2 - s)(x 2 + s 2 ) 
the starting torque as 







(7.29) 



and the starting figure of merit as 

f o TTI 

We find the free-running speed(s) by solving the equation 
2 - s x 2 + s 2 



and selecting the stable solution(s). 
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We can interpret these results further by looking for regions of 
motor, generator, and brake operation, for fixed x and , , as a 
function of s, or shaft speed. We do this under the assumption 
that , is determined independent of s . We further assume, with 
no loss of generality, that < , < 1 . Cases with t> greater than 
unity can be obtained by replacing s with Z - s and with ljt> . 
Physically this means we consider only cases where the machine 
"tends" or "prefers" to go in the positive direction; by reversing 
the definition of mechanical speed we easily treat the other case by 
symmetry. 

For convenience we define the parameter 



X 2 +8 2 

x 2 + (2 - s) 2 



so that 



(7.32) 



(7.33) 



For MOTOR operation, < rj m < 1 . For < s < 1 this occurs 
when ~1<X<1, and for l<s<2 this occurs when JXJ > 1. 
For GENERATOR action, rj m > 1 . For s < this occurs when 
|x| < 1 , and for s > 2 this occurs when |\| > 1 . And finally, 
for BRAKE operation, 77 < 0. For s > 1 this occurs when 



|x| < 1 , and for s < 1 this occurs when 
tions are summarized in Table 7.3. 



x > 1 . These condi- 







M <i 


|\| > 1 


co m < -^s 


2 < s 


BRAKE 


GENERATOR 


-co < w < 
s m 


1 < s < 2 


BRAKE 


MOTOR 


< o> m < u g 


< s < 1 


MOTOR 


BRAKE 


w s <w m 


s <0 


GENERATOR 


BRAKE 



Table 7.3. Unbalanced induction machine operation, 
in terms of the parameters s and X 

We have reduced the problem to one of determining when X is 
less than one in magnitude. A few facts about X(s) bring this out; 
X(s) has a zero at s = 2, and a pole at s = . Furthermore, X(s) 
is the ratio of two cubic expressions in s, and so has only five 
points of zero slope. Since the points of zero slope are also points 
of zero slope for I/X(s) or 2 /X(s) = X(2 - s) , we see that these 
points are located symmetrically about the point s = 1 . 



114 



r e quency ~ po we r formulas 



For high values o s, X approaches - & from below, and for 
high negative values of s, X approaches - from above. Because 
X is negative except between and 2, there must be exactly one 
point of zero slope between -co and 0, and exactly one point of 
zero slope between 2 and +co , for otherwise there would have to 
be three such points on either side, or a total of six too many. And 
between and 2 there may be either two points of zero slope, or 
no such points. 

Because by assumption & < 1 , there must be exkctly one (no more 
no less) point where X = - 1 in the region - oo < s < . Since X(l) = 
< 1 , and because there is at most one point of zero slope between 
and 1, X = 1 exactly once in this region. For I < s < 2, X 
may equal unity either two times or no times, and for 2 < s < +00, 
X = -1 either two times or no times. 

Using this reasoning, we subdivide the range of speeds into re- 
gions of operation. The most general type of subdivision of this 
sort (assuming that is independent of s ) is shown in Fig. 7.2. 
Depending on the values of x and , regions A and 3 may or 
may not be present. No attempt has been made to have Fig. 7. 2, 
or any other such diagram, drawn to scale. 
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Fig. 7.2. Most general division into regions of opera- 
tion, under the assumption that x and are 
independent of the slip 

As explained in Sec. 7.2 (Almost-Balanced Operation), the dia- 
gram of Fig. 7.2 reduces to that of Fig. 7.1 as balanced operation 
is approached. 

The free-running speeds w- are indicated in Fig. 7.2; the point 
separating regions S and II is also a point of zero torque, but it 
is not stable. 

Some of the regions of operation shown in Fig. 7.2 are well known. 
It is not known to what extent, if any, the others, especially the gen- 
erating and braking regions at high mechanical speed, are important 

Almost-Balanced Operation. This analysis gives a glimpse into 
the operation of nearly -balanced machines. If a machine is intended 
to be operated balanced, but there is a defect, either in the device 
or in the excitation, the operation might be described by a value of 
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close to, but not quite, zero. 

The various regions of operation for the almost-balanced machine 
are shown in Fig. 7.3 e The only difference between Fig. 7. 3 and 
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Fig. 7.3. Most general break-up into regions of 

operation for an almost-balanced machine 

Fig. 7.1 (for the balanced machine) is the new braking region close 
to synchronous speed. The torque, starting torque, starting figure 
of merit, and efficiency are all continuous functions of , but not 
uniformly continuous. The almost-balanced machine approximates 
the balanced machine for t> small enough, except near regions where 
the functions have possible singularities. The only such region is 
close to synchronous speed, and as we have seen, there is an im- 
portant qualitative difference there. Aside from this, however, the 
balanced case is obtained in the limit as the unbalance disappears. 

In summary, a slight unbalance introduces a small braking re- 
gion near synchronous speed. The free-running speed is slightly 
less than synchronous speed. The same effect is obtained when 
there are small losses in the system, so in practical cases an ob- 
served small braking region like this might not be completely elim- 
inated by merely reducing the unbalance. 

The Completely-Unbalanced Machine. By "completely -unbalanced' 1 
we mean that the air-gap field due to stator currents is a purely 
pulsating field, so = 1 . This occurs when only one stator wind- 
ing is excited or when the current through all stator windings is of 
the same phase. 

With = 1, the motor efficiency reduces from Eq. 7.27 to 

(7.34) 



" v ' 8 (2-s)+x i 
and the torque becomes, from Eq. 7.28 

D s (1 - S) a(2 - s) - x 2 



(7.35) 



s(2 - s) + x 2 



The starting torque and starting torque figure o merit both vanish. 
The freerunning speed is found by setting the torque equal to zero, 
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for nonzero o> , that is, by solving the equation 



s(2 - s) x 2 



(7.36) 



so 



(7.37) 



In this special case of completely unbalanced operation, some of 
the regions of Fig. 7.2 disappear. The machine is symmetrical in 
the sense that if it operates in one mode at speed co , the opera- 
tion at speed -to m will be identical; there is no "preferred" di- 
rection of rotation. As a consequence, regions F and II mujst 
vanish; the most general chart of the form of Fig. 7.2 for the com- 
pletely unbalanced machine is Fig. 7.4. The regions 3 and S van- 
ish if x > 1 , orifR>coL. Note the free-running speeds co are 
indicated. 
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Fig. 7.4. Most general division into regions of opera- 
tion for a completely unbalanced machine. The 
free-running speeds to- are indicated. Note 
the symmetry about the point s = I 

7.3 The d-c Induction Machine 

The d-c induction machine [11 6, Sec. 4.4.2] uses a commutator in 
an essential way, but nevertheless we can derive an expression for 
efficiency of the device from the Manley-Rowe formulas. We must 
be careful to define our system in such a way that the commutator 
and brushes are outside. The commutator is used to convert fre- 
quency; hence the frequencies applied to the system we are allowed 
to consider, do not coincide with the frequencies applied to the ac- 
tual machine. Direct current is applied to the rotor brushes in the 
actual operation, but we have to consider the power input at these 
terminals to be a-c. 

In the operation of the device, the brush carriage, with position 
indicated by i|j(t) , is rotated at a constant rate, so that 

iKt) = to^t (7.38) 

and the rotor, whose position is <|>{t) , also rotates at a constant 
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speed w m . If a d-c current I is injected into the "Y rotor brushes 
and the 6 rotor brushes are left open-circuited, there will be a 
net torque developed and the device will operate as a motor, pro- 
vided the stator windings are terminated in resistors. 

We now inquire as to the power relations. The current and volt- 
ages in the stator will be of frequency the same as that supplied by 
the rotating field, o^ Furthermore, the motor driving the brush 
carriage does no work on the system other than to overcome fric- 
tional losses, which in this analysis we neglect anyway, so the power 
input from that source can be neglected. 

The current through a given rotor wire will not be direct current, 
in spite of the fact that the rotor is supplied with d-c. The relative 
speed of the commutator brushes over the rotor is merely o>, - o , 
and the current each rotor winding experiences is a square wave 
with this fundamental frequency. Thus for application of the Mau- 
ley -Rowe formulas to the d-c induction motor, we must consider 
the rotor power to be at frequency co - w and all its odd harmon- 
ics. 

However, in applying the formulas it is known that we may con- 
sider power at harmonics of a given frequency to be power at that 
given frequency (see Appendix C). Thus we consider all the rotor 
power input to be at the fundamental frequency 01, - o> , and not 
spread among the various harmonics, as it physically is. 

Our procedure is then clear. Power to the rotor we consider to 
be at frequency w r = o>^ - co m ; power to the stator is at frequency 
co, ; and power at the shaft is d-c power, which by Appendix D we 
may consider to be power at frequency co m . 

Let the rotor, stator, and shaft power inputs be, respectively, 

P_ , P , and P . Then the Manley-Rowe formulas are 

r s m 

"P "P 

_i+_i = o (7.39) 

b "m 

and 

P- P^ 

-* _ -H (7.40) 



For motor operation we want power to flow in only at the rotor, 
o that P r > , 
.40 shows that 1 
ficiency 77 is 



so that P > 0, P < 0, and P < . A glance at Eqs. 7.39 and 
7.40 shows that this is possible, for o> m < ; in fact, the motor ef- 



. 



which shows that the faster the rotor turns over, the higher the ef- 
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ficiency. Equation 7.41 further shows that for < o> m < co the de- 
vice runs as a brake, and for co, < u> m it acts as a generator. 

The predicted power relationships are verified by the exact sol- 
ution [116, Sec. 4.4.2] . The power output ( - P ) is given by T &> m , 
where T is the mechanical torque of electrical origin. The stktor 
power output (-P fl ) is given by 2(f)R s | i* | 2 , the factor of 2 arising 
from the two stator windings with identical power, and the factor of 
one-half coining from averaging the sinusoidal waveform of i and 
i . Using the formulas for stator current and torque given in ref- 
erence 116, Eqs. 4.170, 4.171, and 4.172, we can immediately verify 
Eq. 7.39. Equation 7.40 can be verified from the expression for 
rotor input power 

P = v r i r = v r l (7.42) 

r v v Y 

The brush voltage v r is not given explicitly in reference 116, but 
can be calculated easily from Eq. 4.136a of reference 116 and the 
solutions for the stator currents, Eqs. 4.170 and 4.171. It turns out 
to be (neglecting rotor resistance) 



v r = * >* ' _ ux u ~ (7.43) 



from which Eq. 7.40 is immediately verified. 



Chapter 8 
COMMUNICATIONS APPLICATIONS 



Frequency-power formulas are conservation theorems, and are 
of the most value when used as such: to find constants of the motion, 
to base new definitions on, to check analytic results, and, especially, 
to give an intuitive grasp of the frequency conversion operation. 
Only occasionally can the formulas be used to obtain directly results 
of practical interest, and these cases usually involve only a few fre- 
quencies. A few typical applications to communications systems are 
discussed here. We cannot, of course, pretend to discuss all ap- 
plications or even to discuss any thoroughly. 

In Sec. 8.1 the constraints imposed by the formulas on harmonic 
generators are given; it is found from the real-power constraints 
that nonlinear reactance generators are not limited in efficiency by 
the frequency-power formulas, whereas nonlinear resistor genera- 
tors are. These results are not new, but the reactive -power con- 
straints given have not, to our knowledge, been used. 

Many frequency converters are used in communications for pro- 
cessing small signals, and linearized time-varying equations of 
motion are almost universally used in analyzing such applications. 
In Sec. 8. 2 we show how these linear time-varying equations are 
derived for a parametrically-pumped nonlinear system, of the type 
that obeys any of the four types of frequency-power formulas; we 
also give the resulting frequency-power formulas. Many such de- 
vices are analyzed by matrix methods, and in Sec. 8.3 the conver- 
sion matrices are shown to be constrained by the frequency-power 
formulas. In Sec. 8.4. we demonstrate a typical use of frequency- 
power formulas to determine fundamental limits of performance 
of a special type of frequency converter. We do not treat even 
this small problem completely, but merely show one method by 
which the frequency-power formulas can be used. 

Several other interesting applications of the Manley-Rowe form- 
ulas to communications have been made already. We will not dis- 
cuss these in detail, but only briefly mention them here. 

If there are two independent frequencies and only three frequen- 
cies at which power is exchanged, then the frequency-power form- 
ulas are very simple, and in fact ratios of the powers can be cal- 
culated. Rowe[60,89], Uhlir[l07], Duinker [18] , and others have 
discussed, from this point of view, the use of nonlinear reactances 
as parametric amplifiers, upconverters, modulators, mixers, down- 
converters, etc. 

119 
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Uhlir [107] has made an interesting observation about the use of 
the Manley-Rowe formulas in three -frequency devices. He com- 
pares them to other "second laws, " such as the second law of ther- 
modynamics, or the law of conservation of momentum in classical 
particle dynamics. These "second laws, " together with the law of 
conservation of energy, lead directly to expressions of interest, 
such as efficiencies, in the same way the Manley-Rowe formulas 
do. For example, one can combine the first and second laws of 
thermodynamics into a form similar to Eq. 1.12. 

Heat Input _ Work Output _ Heat Rejected (8.1) 

= = 



for a reversible Carnot engine, where T^ and T C are the tem- 
peratures of the hot and cold resevoirs. 

The theory of noise in parametric amplifiers is both simplified 
and made more general by using the Manley-Rowe formulas [40, 93]. 
The theory is simplified because the formulas can be used in place 
of detailed solutions, and it is made more general because specific 
circuits need not be used, and because the results are shown to hold 
for any physical system that satisfies the Manley-Rowe formulas. 

Kurokawa and Hamasaki [50, 49] have made interesting use of the 
formulas by showing that they predict a type of orthogonality of sol- 
utions of a par ametrically- pumped distributed system, in a way 
similar to the way the ordinary Poynting's theorem predicts orth- 
ogonality of modes in lossless distributed systems [4], The nonlinear 
theory of a limiter made from a weakly-pumped reactance was dis- 
cussed by Olson, Wang, and Wade [72], who used the Manley-Rowe 
formulas to simplify their derivation. 

One of the most interesting applications is to magnetic amplifiers. 
A magnetic amplifier can be thought of as a combination of a non- 
linear reactance upconverter and a nonlinear resistance downcon- 
verter, and Mauley [58] has published a discussion along these lines, 
using his formulas. 

8.1 Harmonic Generators 



Nonlinear elements can be used as harmonic generators, and it 
is interesting to ask what constraints, if any, the frequency- power 
formulas impose. Power is put into the device only at frequency co, 
and we assume there is no other independent frequency present (note 
that this last statement is an assumption, not a statement of fact). 
Power then flows out only at frequencies of the form ko>, with k 
positive. If P. and Cl are the real and reactive powers into the 
device at frequency ku, the four frequency -power formulas reduce 
to 

Type I: S fc P k = (8.2) 

Type II: S k kQ k = (8 * 3) 
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Type in (Page form): S P (1 - coske) > (8.4) 

K. K. 

Type III (Pantell form) : S k 2 P >0 (8.5) 

Type IV (Inductive): S Q (1 - cos k6)/k > (8.6) 

S k kQ k~ {8 ' 7) 

Type IV (Capacitive): ^v^V 1 " coske ^ k - * 8 - 8 ) 

S kCI < (8.9) 

Jx Ix 

Equation 8.2 is merely the law of conservation of energy, and the 
Manley-Rowe formulas have nothing more than this to say. Because 
of this, it is often thought that efficiency arbitrarily close to 100% 
is possible for a nonlinear capacitor harmonic generator[l3, 52, 53] . 
Actually, such a conclusion must [46] come from a circuit analysis; 
the Manley-Rowe formulas do not say that this is possible they 
merely do not say that it is impossible. 

On the other hand, for nonlinear resistance harmonic generators, 
the frequency-power formulas of the third type, Eq. 8.5, give a def- 
inite restriction on harmonic generation efficiency. The power out 
-P at frequency nco is bounded by Eq. 8.5: 

S k k2p k p 

- ' ' f (8.10) 



where the sum is over all k except n; thus the conversion effi- 
ciency is limited to (1/n 2 ) , as reported specifically by Page[74,75] 
and noted by others. 

As far as we know, the second and fourth types of frequency-power 
formulas have never been applied to harmonic generators. 

Subharmonics can be generated with nonlinear reactance s[52, 53]; 
if power is put in at frequency co , power can be extracted at fre- 
quencies co/2, co/3, etc. Note that the Manley-Rowe formulas, Eq. 
8.2, do not prevent this. Much of the parametric excitation observed 
in the 19th century is in reality subharmonic generation. The para- 
metron computing element [109, 31,117] is one practical device that 
makes use of this phenomenon. 

Subharmonic generation is impossible, however, with nonlinear 
resistors [74] that obey frequency-power formulas of Type HE. Pan- 
tell 1 s form does not predict that it is impossible, but Page's form 
does. Setting 6 = 2ir in Eq. 8.4, we find a weighted sum of power 
inputs at all frequencies except co and its harmonics, is positive. 
Because by assumption power is put in only at frequency co, it can 
therefore come out only at harmonics of GJ . Moreover, as Page has 
shown [74] , if power is put in a device that satisfies frequency-power 
formulas of Type III only at the two frequencies tdj and o> , power 
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can come out only at frequencies of the form mtOj + nw Q , with m 
and n integers. 

8.2 Linear Parametrically-Pumped System 

A nonlinear system of the type discussed in this report, when strongly 
pumped, appears, to superimposed small signals, just like a linear 
time -varying system. In fact, one can even calculate the parameters 
of the time -varying system, from a knowledge of the actual nonlinear 
system and the pumping. As an example of this, we consider a case 
of some practical interest, a nonlinear capacitor pumped by a volt- 
age e (t) , given as a function of time. The treatment can then be 
generalized. We will see that if the overall nonlinear system can be 
shown to obey any of the four types of frequency -power formulas by 
the proofs of Chapter 3, then the equivalent time-varying system 
obeys formulas of the same type. 

The nonlinear capacitor is characterized by its charge q being a 
function of the voltage e ; let us define the "incremental capacitance" 

Cl (e) - fj (8.11) 

as a function of e . Now let the unperturbed variables (voltage, 
charge, and current), as pumped, have a subscript p (for pump), 
and the small perturbations have a subscript s (for small-signal); 
thus 

e(t) = e p (t) + e g (t) (8.12) 

q(t) = yt) + q g (t) (8.13) 

and 

i(t) = i (t) + i fl (t) (8.14) 

where 

dq (t) 



Our goal is to obtain equations of motion involving the small- signal 
variables alone. 

This is easily done; we expand the charge in a Taylor's series 
about the unperturbed solution: 

q(t) = 

where the subscript p on the derivative denotes evaluation at the 
unperturbed solution. We note from Eqs. 8.13 and 8.15, then, that 

(8.17) 
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if the small signals are small enough, which is in the desired form 

* = c(t) 



if we identify the equivalent time -varying capacitance with the in- 
cremental capacitance, as pumped, 

C(t) = C.(e (t)) (8.19) 

We have succeeded, then, in replacing a heavily -pumped nonlinear 
system with a linear time-varying system that is equivalent, as far 
as small perturbations are concerned. It is clear that we cannot, 
within the framework of the time -varying model alone, predict its 
range of validity, nor can we predict any of the very interesting 
limiting and saturating properties the finitely -pumped nonlinear 
system would exhibit. 

This technique can also be used with more complicated systems. 
It is our intent to show that it can be used for any system for which 
any of four proofs in Chapter 3 apply. Although we use notation ap- 
propriate to systems that obey the first type of frequency-power 
formula, the reader can easily check that the following paragraph 
also applies to systems that can be discussed in Sees. 3.3, 3.4, or 
3.5. 

Suppose we have a number of variables x- , on which an equal 
number of variables f . depend. Then, suppose the system de- 
scribed by these variables is pumped with the variables x i p W ai1 ^ 
f - (t) . We wish to replace the heavily-pumped nonlinear system 
with a linear time-varying system that is equivalent as far as the 
small- signal variables x ig (t) and fj s (t) are concerned. This is 
easily done by writing a Taylor's series expansion for f,(t) = fj p W +fj g (t): 

x. s (t) + ... (8.20) 

Taking only the first-order terms, we find linear time-varying equa- 
tions of motion, 



which accurately describe the small perturbations of the system if 
they are small enough, and if the higher derivatives of the f . with 
respect to the x i are small enough. Note in particular that the per- 
turbations must be small for all time, so that if the X|(t) are of 
the form of Eq e 3.17, then v. , the average value of dx i (t)/dt l must 
not change, for otherwise x is would not be bounded. 

We could have allowed significant space variation of the variables, 
so this derivation is not limited to lumped systems. In the descrip- 
tion of many systems in Chapter 6, however, physically-unmeasur- 
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able potentials and Lagrange multipliers are used, and it is not cer- 
tain that these variables can be decomposed into an unperturbed part 
and a bounded perturbation, even if the measurable variables can be. 
One must therefore use caution in dealing with pumped distributed 
systems. 

We now wish to show that the time-varying model thus produced 
obeys the same type of frequency-power formulas that the pumped 
nonlinear system does; or, more precisely, that if any of the four 
proofs of Chapter 3 holds for the nonlinear system, then the same 
proof holds for the equivalent time-varying system, with the ex- 
pressions for real or reactive power being of the same form. We 
do this for each of the four types in turn. 

First, suppose the system has an energy state function U(x^) , 
with 



(8 - 22) 



so that the proof of the Manley-Rowe formulas of Sec. 3.2 holds. In 
that case, 



and we may evaluate the line integral 

S.ff. dx. (8.24) 

1J 18 IS 

using Eq. 8.21, by any path and get the same result, 

U = fS.f. x. 
t i is is 

8f. \ 

i- x. x. (8.25) 

&K / 1S J s 

Thus U t is a state function of the x^ g and time explicitly (through 
the pumping), with 

au (x. , t) 
f js " ax." ' 8 - 26 > 

so the time-varying small- signal system must obey the first type of 
frequency -power formula, with the P given by the same form, Eq. 
3.25, but with subscripts s on the variables: 

P a = iRe 2. J a tf u )>!.>a (8.27) 
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It does not, of course, obey any such formula corresponding to an 
independent frequency on which the pump waveform depends. Note 
that there is no d-c contribution to the time-varying systemfs form- 
ulas. 

When these steps are applied to a distributed system that obeys 
Hamilton 1 s principle, it is found that the resulting P a can always 
be expressed as the divergence of a vector of the form of Eq. 6.132, 
but involving only the small- signal components of the variables: 

PQ = -V T a (8.28) 

where each component of J? is 




The same sort of analysis holds for systems known to obey the 
second type of frequency-power formula. Here, the dissipation 
function is G(v.) , with 

8G(v.) 
f. = - (8.30) 



and therefore, 

8f . \ / 8f . 

1 ^ - I * \ (8.31) 



Thus the dissipation function 

G = iS.f. v. 

t i is is 



8f . \ 

_L] v v (8.32) 

8v. / is js 
3 P 

is a state function of the v ig and time explicitly, with 

8G (v. , t) 

t 1S = f (8.33) 



8v. g js 



and so the linear time-varying system obeys the second type of fre- 
quency-power formula, with the Q a of the form of Eq. 3.26 with 
s subscripts: 



Q = 4lmS.(f. ) *(v. ) 
a 2 i v is'a is a 



Also, if the original nonlinear system obeys the third type of fre- 
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quency-power formula, the time-varying model for the pumped non- 
linear system also obeys this type, provided the pump does not force 
the nonlinear system into regions in which the matrix of Eq. 3.51 is 
not positive definite. The expression for the small- signal variables 

/8f. \ 
f. = S . (-J- v. (8.35) 

js ^Vp 1S 

gives each f. as a function of the v. and time explicitly (through 
the known functions (afr/av^ ). It is then known from the proof of 
Sec. 3 4 that frequency-power formulas of the third type hold, pro- 
vided the matrix with entries 

8. (v. , t) / 8f . \ 

-11-12 = Ui (8.36) 

8v is ^Vp 

is positive definite, which it is if the pump does not force the orig- 
inal nonlinear system into regions where the matrix of Eq. 3.51 is 
not positive definite. 

And finally, a system that is shown to obey the frequency -power 
formulas of the fourth type by the proof of Sec. 3.5, if pumped heav- 
ily, can be replaced for small signals by a linear time -varying 
model; and this model itself obeys the same type of formulas pro- 
vided the pump does not force the original nonlinear system into 
regions where the matrix of Eq. 3.66 is not positive definite. 

These results can be summarized as follows: 

Let a system that can be shown by the methods of Chapter 
3 to obey one type of frequency -power formulas be pumped 
strongly. Then, to small signals, the pumped nonlinear 
system behaves like a linear time -vary ing system, whose 
parameters can be calculated from a knowledge of the 
nonlinear system and the pumping. Also, the equivalent 
linear time-varying system obeys frequency -power form- 
ulas of the same type, where the expression for real or 
reactive power is of the same form, but is calculated 
from the first-order small-signal variables alone. 

With this theorem established, we can write explicitly the 
four types of frequency-power formulas for linear, parametric ally- 
pumped systems. Let us suppose that pumping is at frequency co 
and its harmonics, and that signal and generated sidebands appear 
only at frequencies of the form ko> + o> , where k is positive or 
negative. Denoting the real and reactive power at frequency kco -f w 
by P and Q, respectively, we have 

P 

Type I: 2 = (8.37) 

k kto + co 
p s 



matrix constraints J27 

Type II: S k Q k = ( 8 - 38 ) 

Type III: S k P k - (8a39) 

a 

Type IV (Inductive): S ; > (8 40) 

k kco + co ; 

P s 

S k 

Type IV (Capacitive): S 1 < (8.41) 

k kto + co ~~ ' 

P s 

Note that in Eq. 8.37 there is no contribution from d-c power, be- 
cause the small- signal components of the variables necessary to 
generate it must vanish. 

If k is far enough negative so that kco + o) g is negative, then we 
may wish to sum over positive frequencies. This may be done in 
Eqs, 8.37, 8.39, 8.40, and 8.41, because the real power associated 
with a negative frequency is the same as the real power associated 
with the same positive frequency. However, Eq. 8.38 must be re- 
placed by 

kco + co 
S_ i % L_ Q - o (8.42) 

k ]kw +w | k 
p si 

if we are to sum over the corresponding positive frequencies, be- 
cause the reactive power associated with a negative frequency is the 
negative of the reactive power associated with the corresponding 
positive frequency. In Eq. 8.42 the Q, are to be interpreted as 
reactive power inputs at frequency | kco + co g | . 

8.3 Matrix Constraints 

Much linear time -invariant single -frequency circuit analysis is 
done by matrix methods. For example, to describe an n-port linear 
time -invariant network such as the one in Fig. 8.1, we may use an 
n X n impedance matrix, or admittance matrix, or scattering ma- 
trix, or other type of matrix. The same techniques are useful in 
describing linear par ametrically -pumped networks. 

Just as conditions of losslessness, passivity, etc., place con- 
straints on the matrices that describe linear time -invariant net- 
works, so the appropriate frequency- power formulas place constraints 
on the matrices that describe time-varying networks of certain types. 
In Sec. 8.3 (Time -Invariant Matrix Constraints), we derive the well 
known matrix constraints for linear time -invariant networks of var- 
ious types, as an indication of the method used. Then in Sec. 8. 3 
(Linear Time -Varying Matrix Constraints), we derive the constraints 
on matrices for time -varying systems of the various types. Some 
of these constraints can be expressed in a variety of forms, which 
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we tabulate for convenience. A few of these time-varying system 
matrix constraints have been used in the past for special purposes 
[40, 93, 50, 49] . 



I, 



LINEAR 
TIME- 
INVARIANT 
NETWORK 



Fig. 8.1. Schematic representation of a linear time- 
invariant system, with n ports. Although 
electrical variables are shown, some of the 
ports may be mechanical 

Time -Invariant Matrix Constraints. Consider the network of Fig. 
8.1. It has n ports, each of which has a voltage of frequency co 
and (complex) amplitude E i (i = 1, 2, . . . n) , and current of am- 
plitude L. If we define the column matrix E with entries E. , in 
order, and the column matrix ^1 with entries I. , in order, we may 
express the linear equations of motion of the terminal variables in 
the form 



E = ZI 



(8.43) 



where the n X n square matrix Z is the impedance matrix of the 
network. The admittance matrix' is Y = Z"" 1 



I = YE 



(8.44) 



Alternatively, we may define the scattering variables for each 
port 



a. = 



Z.I. 
11 



(8.45) 



and 



V. - Z.*I. 
i 11 



Z.* 



(8.46) 



where Z i is an arbitrary complex number with positive real part 
(often, for convenience, Z i is taken to be the terminating imped- 
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ance at port i). The linear equations of motion can then be written 
in the form 

b = Sa (8.47) 

where b is the column matrix of b. , a is the column matrix of 
a. , and the n X n square matrix S is called the scattering matrix 

[in- 

The total power input P is the sum of the power inputs at each 
port, or 

P s f Re S. E.I.* 

* 111 

= JS (E I * + E *I) (8.48) 

4 i i i i i 

which can be written in matrix form as 
P = i(ETi+ Itg) 

= JET (Yt + Y)E = ilt (Zt + Z)I (8.49) 

where the dagger indicates the complex conjugate transpose of a 
matrix; or P can be written in terms of the scattering variables 
as 

P = S.(|a.| 2 - b.| 2 ) 



= at(l - STS)a (8.50) 

where 1 is the unit matrix. Similarly, the reactive power input Q 
can be written 



= 1 Et (Yt - Y)E = -~r It (Z - Zt)I (8.51) 

4j ~ ~ ^J ^ 

The properties of the network place constraints on the matrices 
Z, Y, and S. First, suppose the network is lossless. Then P 
= r for any "choice of E, , or a, so that from Eq. 8.49 

Z + Zt = (8-52) 

and 

Y + Yt =0 * 8 - 53) 

and from Eq. 8.50 

1 - S1S = < 8 - 54 ) 

or S is unitary [43, Sec. 1.16]. Note that Eqs. 8,52, 8.53, and 8.54 
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do not require that Z , Y , or S be symmetrical. 

Secondly, suppose that the network is pure lossy, in which case 
the total reactive power input Q vanishes. Then, from Eq. 8. 51, 

Z = Zt (8.55) 

and 

Y = Yt (8.56) 

or both Z and Y are Hermitian [43, Sec. 1.16] . 

As a third example, suppose that the network is passive. Then P 
is nonnegative, so that from Eq. 8.49 

Z + Zt > (8.57) 

and 

Y + Yt > (8.58) 

where we mean in setting a matrix greater than or equal to zero, 
that it be positive definite [43, Sec. 1. 17] . Similarly, from Eq. 8.50, 
the scattering matrix describing a passive system fulfills the con- 
dition 

1 - StS > (8.59) 

And finally, a network composed of only resistances and induc- 
tances has a positive reactive power input, so from Eq. 8.51 

j(Zt- Z) > (8.60) 

and 

j(Y - Yt) > (8.61) 

Similarly, a capacitive network has a negative Q, so 

j(Z - Zt) > (8.62) 

and 

j(Yt- D > (8.63) 

These constraints for special types of linear time -invariant cir- 
cuits are given in Table 8.1. The constraints are well known; they 
were derived here primarily to indicate the method used in the next 
section for networks obeying one of the four types of frequency-power 
formulas. 

Linear Time-Varying Matrix Constraints. Suppose a linear para- 
metrically-pumped system, of the type described in Sec. 8.2, has a 
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Lossless 


P - 


z = -zt 

Y = -Yt 
S" 1 = Sf 


Pure Lossy 


Q = 


Z = Zt 

y = Y| 


Passive 


P > 


z + zt > o 

Y + Yt > o 

1 - StS > 
1 - SS t > 


Inductive 


Q > 


j(Zt - Z) > 
j(Y - Yt) > 


Capacitive 


Q < 


j(Z - Zt) > 

j(Yt "Y) > o 



Table 8.1. Matrix constraints for linear time-invariant 
networks of various types 

finite number of ports, as shown in Fig. 8.2. In general, each port 
exchanges power at a number of frequencies. Let us represent this 
system by Fig. 8.3, where we show extra terminal pairs, so that 
each terminal pair now exchanges energy at only one frequency. 



LINEAR 

PARAMETRICALLY- 
PUMPED 
NETWORK 



Fig. 8.2. A linear parametrically- 
pumped system, with a small 
number of physical ports. Each 
port exchanges power at a num- 
ber of frequencies 



LINEAR 

PA R A METRIC ALLY - 
PUMPED 
NETWORK 



Fig. 8.3. The same system shown 
in Fig. 8.2, but with additional 
fictitious ports added, so that 
each new port exchanges power 
at only one frequency 
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This is a nonphysical representation for the network, in the sense 
that one physical element may appear more than once if currents of 
various frequencies flow through it. In addition, the loads placed on 
the ports are not all physically separate, so any given physical ele- 
ment may appear in several of them. 

For example, consider the network of Fig. 8.4, which is a rea- 
sonably good approximation of a par ametrically- pumped varactor 
[107], including the series resistance. This physical device has 
only one terminal pair, but when it is used as a frequency converter 
it is convenient to represent it as a poly-port network. The network 
of Fig. 8.5 would be suitable if power is exchanged at only two fre- 
quencies. Note that the series resistance appears at both terminal 
pairs, even though it is only one physical element. 

Assume the pumping is at fre- 
quency co , and the small signals 
at frequencies (kco + to g ) . Be- 
cause the system of Fig. 8. 3 is 
linear, the equations relating the 
terminal variables are linear and 
can be put in matrix form. Each 
terminal pair has voltage and cur- 
rent with frequency co- (of the 



C(t) 



Fig. 8.4. Model for a parametric - 
ally pumped varactor, includ- 
ing series resistance 



form kco + co ) and amplitudes E. 
and I. . Then we may define E 
and ,[ as column matrices of the 
Ej and I., and write the "conver- 
sion matrices" Z and Y 



that [10 6, Ch. 5] 



Fig. 8.5. A schematic representa- 
tion for the parametrically- 
pumped varactor. Each port 
shown is at a different fre- 
quency; the voltage and cur- 
rent at each port are merely 
the components of actual var- 
actor voltage and current at 
the corresponding frequency. 
Note that because the series 
resistance is linear, it may be 
put separately in each lead, so 
that the box shown consists 
only of the varying capacitance, 
and therefore satisfies fre- 
quency-power formulas of 
Types I and IV 



and 



E = ZI 



I = YE 



such 



(8.64) 



(8.65) 



Alternatively, we can define the 
scattering variables at each port 
by Eqs. 8.45 and 8.46, and write 
the equations of motion in the form 

b = Sa (8.66) 

using the scattering matrix S. Al- 
though Z, Y, and S describe a 
time-varying system, they do not 
vary in time themselves. 

If the network of Fig. 8.3 obeys 
frequency-power formulas of any 
type, these place constraints on the 
matrices Z, Y, and S, similar 
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to the constraints of Table 8.1. We now calculate these matrix con- 
straints, and tabulate them. 

Suppose the network of Fig. 8.3 consists of reactive-type elements, 
gyrators, transformers, and other devices that obey the Manley- 
Rowe formulas. Then we write, from Eq. 8.37, 

iRe(E.I.*) 
TJ i-i_ = (8.67) 



w i 



where the sum over frequencies is replaced by a sum over ports, be- 
cause each port passes power at only one frequency. In terms of the 
diagonal matrix K with entries 1/c^ along the main diagonal, Eq. 
8.67 may be written 

i (ItKE + ETKJ) = ( 8 - 68 ) 

or 

= rf(KZ + ZtK)I 

= Et(YtK +KY)E ( 8 69 > 

and because I and E can be arbitrary in Eq. 8.69, 

KZ + Z-fK = ( 8 ' 7 ) 

and 

KY + YtK = < 8 - 71) 

In addition, Eq. 8.67 can be written 



using the scattering variables, so 
= atKa - btKb 

= at(K - StKS)a 
whereupon we conclude that 



Equations 8.70, 8.71, and 8.74 are the constraints the Manley- 
Rowe formulas place on the system matrices Z, Y, and S Many 
other equivalent forms can be obtained from these constraints; for 
example, premultiply Eq. 8.74 by IT'Sf 1 and postmultxply it by 
K^S-", to obtain 



(8 - 75) 



a form that has been used before [40] . In Table 8.2 ^we list several 
variations of the constraints of Eqs. 8.70, 8.71, and 8.74. 
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KZ + Z-fK * 


KY + YtK = 




Z = -K" x ZtK 


Y = -K^^tK 




Zt = -KZK" 1 


Yt = -KYK*" 1 




(Z) n = (--ifK^CZtfK 


af-<- 1 fK-. ( vt,-K 


Type I 


K = -ZtKY 


K" 1 = -ZK^Yt 


(Manley-Rowe 
Formulas) 


= -YtKZ 


= -I5- 1 5T 




= (~lf(ZtfK(Y) n 


I "\\ t "V \ "u^ "* / *7 4* \ 




K = StKS 


K" 1 = SK -1 St 




S = K-'(S t )-'K 


St = gS- 1 ^ 1 




(S) n = K^Stf 1 ^ 


(Stf = KCS)"^- 1 


Type II 


Y = Yt 


z = zt 


Type III 


Z + Zt > 
Y + Yt > 


1 - sst > o 

I - gtS > 


Type IV Inductive 


j{ZtK - KZ) > 


j(KY - YtK) > 


Type IV Capacitive 


j(KZ - ZtK) > 


j(Ytg - KY) > 


Types I and IV 
Inductive 


jZtK = -JKZ > 


jKY = -jYtK > 


Types I and IV 
Capacitive 


JKZ = -jgtK > 


jYtK = -jKY > 


Types Hand III 


Z = Zt > 


X = Xt > o 



Table 8.2. Various forms of the matrix constraints for 
linear time-varying circuits that obey one or 
more frequency-power formulas; n is a pos- 
itive or negative integer 
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Now suppose the network of Fig. 8.3 obeys frequency-power form- 
ulas of the second type. Then by Eq. 8.38 the sum of the reactive 
powers in at each port vanishes; thus 

= S.Q. = ^T E| (Yt - Y)E 

= jr It (Z - Zt)I (8.76) 

^J 

and hence 

Z = Zt (8.77) 

and 

Y = Yt (8.78) 

so both Z and Y are Hermitian. Note that some of the frequen- 
cies co^ may be negative, and for these E^ and 1^ are the complex 
conjugates of the amplitudes of the corresponding positive frequency, 
and the matrices Z and Y must be written accordingly, for Eqs. 
8.77 and 8,78 to hold. 

Now suppose the network of Fig. 8.3 obeys the third type of fre- 
quency-power formulas. Then Eq. 8.39 predicts that the total real 
power input is positive. Just as for passive networks, this condi- 
tion predicts that 

Z + Zt > (8.79) 

Y + Yt > (8.80) 
and 

1 - t > (8.81) 

Now suppose the network of Fig. 8.3 obeys the fourth type of fre- 
quency-power formulas, and is inductive. Then, from Eq. 8.40, 

S. > (8.82) 

iw. - 

and the matrix constraints are 

j(Z|K - KZ) > (8.83) 

and 

j(KY - YtK) > (8.84) 

On the other hand, if the system obeys these formulas but is capaci- 
tive, the matrix constraints become, from Eq. 8.41, 
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and 



j(KZ - ZtK) > 



j(JtK -KY) > 



(8.85) 



(8.86) 



We have found matrix constraints imposed by each of the four types 
of frequency-power formulas. These are summarized in Table 8. 2. 

Many devices obey more than one type, however, and it is inter- 
esting to note the constraints that apply to them. If a device satis- 
fies formulas of Types I and IV simultaneously, then, if it is induc- 
tive, the matrices -JKZ and jKY are Hermitian and positive def- 
inite; if it is capacitive, the matrices JKZ and -JKY are Hermitian 
and positive definite. Similarly, if a system obeys frequency -power 
formulas of both Types n and III, then Z and Y are Hermitian and 
positive definite. 

8.4 Impedance Translation 

It is our intent here to mention briefly one typical application in- 
volving frequency-power formulas of more than one type. We do not 
propose to discuss it thoroughly, but only to show that the frequency- 
power formulas can be used to gain insight and even to obtain definite 
fundamental limits on performance. 

Suppose a one -port linear parametrically-pumped network is con- 
nected to a passive load Z^ at one frequency, nw + GO , where GO 
is the pump frequency and GO is another frequency, and suppose 
that the external network is adjusted so that power flows only at fre- 
quencies noo + oo g and Go g . Then, by the method of Sec. 8.3, we 
represent this system by Fig. 8.6, with a terminal pair for each 
frequency. Our desire is to determine the input impedance at the 
co terminal pair. 



Z fi 



LINEAR 

PARAMETRICALLY- 
PUMPED 
NETWORK 



Fig. 8.6. A linear parametrically-pumped one-port 
system exchanging power at only two fre- 
quencies. This device is suitable for im- 
pedance translation 

This device might be used to translate an impedance from one 
frequency to another, for example to obtain an impedance function 
at frequencies close to oo g , that cannot be realized directly with the 
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present state of the art, using convenient elements. Do the frequency- 
power formulas indicate any fundamental limits on the performance 
of such a device? 

At least two applications of this idea have been made, A band-pass 
filter using a similar technique was described by Franks and Witt 
[27]; whereas the simple circuit of Fig. 8.6 does not describe their 
apparatus, it can do so if more terminals are added. And secondly, 
a striking example is the parametric amplifier, in which a passive 
idler impedance Z is translated to an impedance at co with a 
negative real parti 

To see in a simple example how the frequency -power formulas 
might predict fundamental limits of operation for such devices, sup- 
pose the parametrically-pumped network consists of nonlinear re- 
sistors, so that it obeys frequency -power formulas of Types II and 
III. Then, if P and Q are the real and reactive power flowing 
into the network at frequency w , and if P and Q are the real 
and reactive power flowing into the impedance Z^, the frequency- 
power formulas predict 



P s - P n (8 * 87) 



and 



Q s = Q n (8.88) 

These conditions constrain the possible impedances Z we may 
observe at the left-hand terminals, for 

P s + JQ S = iE s l s * 



= H I s| 2 Z s (8.89) 



and 



Thus the real and imaginary parts of Z g are restricted, for 
Re Z ; 



Re Z n (8.91) 
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and 



ImZ s = 




2Q_ 



I 1 .! 2 



frequency-power formula s 



Im 



(8.92) 



The phase angle of Z must therefore lie between zero and the 
phase angle of Z . If we mark the point Z n in the complex plane, 
as in Fig. 8.7, the region of possible values of Z Q is that region to 
the right of the ray passing through Z n , or the shaded area of Fig. 
8.7. 



IMAGINARY 
AXIS 




REAL AXIS 



Fig. 8.7. The complex plane, showing the impedance 
Z n . The frequency-power formulas show 

: that Z is restricted to lie in the shaded 
s 

region 

An interesting result of this is that if the impedance Z^ is that 
of a tuned circuit with a certain quality factor, the quality factor of 
the translated impedance function near frequency oo g is less than 
or at best equal to the original quality factor multiplied by the ratio 



These results are probably not important in themselves, but they 
do serve as an indication of the ability of the frequency-power form- 
ulas to predict fundamental limits. 
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We show here, with a counterexample, that not all lossless devices obey the 
Manley-Rowe formulas; that is, the frequency-power formulas of the first type. 
The device is a simple switch, thrown periodically with a fundamental fre- 
quency w , and connected to a voltage source E cos co t (with oo < w ) and 
a load resistor R. The circuit is shown in Fig. A.I. According to Sec. 5.2, 

the switch does obey the second 
and third types of frequency -power 
formula, but we show by direct 
calculation that it does not satisfy 
the first type. 

The current through the switch 
is (E/R)cos w t when the switch 
is closed and when the switch 
is open. Thus if 8(t) is a func- 
tion that is 1 when the switch is 
closed and when the switch is 
open, the current i(t) is 




Fig. A.I. Circuit used to show that 
a switch does not obey the Man- 
ley-Rowe formulas, even though 
it is lossless 



i(t) = 6(t) - cos co t 
R s 

and the switch voltage e(t) is 

e(t) = [1 - 6(t>] E cos tot 



(A.I) 



(A. 2) 



If the switch is open half the time and closed half the time, and is closed 
during a region symmetrical about the origin, a simple Fourier analysis shows 



9(t) = i + S 8 n cos aco t 
n=l 



(A. 3) 



where 



_2_ ( 
mr 



n odd 



6 = 



Thus 



and 



E6 



i(t) = -^~ cos w t + S -=- cos(nw 4- 5 )t + S -^- cos(nw - 



2R 



n=:1 



(A.4) 



(A. 5} 



e(t) = cos w_t - S 8 co 
s nsl ^ " 



8 n cos(no) - w g )t 

L p 



(A. 6) 
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Thus power flows into the switch at frequency o> s of magnitude 



P s = 8R 

and out at all frequencies of the form (nw w g ) of magnitude 



(A. 7) 



(n odd) 



(n even) 



(A.8) 



V. 

The reactive power input at each frequency is . 

Choosing w g and co as independent frequencies, the first type of frequency- 
power formula would predict 






n=l 




= 



and possibly 



00 

2 

n=l 



nP 




s 



(A.9) 



(A. 10) 



neither of which is true. 

Note, however, that the device does obey the second and third types of fre- 
quency-power formulas, as we would expect because its operation is defined 
in terms of its voltage and current. It is the limit of a time-varying resistor. 
Because the time -variation proceeds at frequency w (and its harmonics), 
there is only one frequency-power formula of Type if: 



00 

S 
n=l 



ap+s 



00 

- S < 
n-1 



and only one of Type III: 

00 00 



n=1 V- 



(A.11) 



(A. 12) 



Each of these is evidently consistent with the actual results (see Eqs, A,7 and 
A.8, and the remark following Eq. A.8). 

Other lossless devices that are known not to obey the Manley-Rowe formu- 
las include the ideal diode, variable transformer, and their mechanical coun- 
terparts. 
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Here we tabulate for convenience the four types of frequency- power formu- 
las. In Table B.I the general formulas are given, without any particular fre- 
quency constraints. In Table B.2 these formulas are specialized to the case 
of most practical interest, for frequencies of the form mWj + nco .If only fre- 
quencies of the form no> are present, the formulas assume the expressions 
of Table B.3. Linear time-varying systems with frequencies of the form nca 
+ o> , where n is positive or negative, might obey the formulas in Table B.4, 
where we assume that to is the pump frequency, and is the only frequency 
that enters into the explicit time dependence of the system; there is only one 
formula of each type. 

If only a few frequencies are present, the formulas are very simple. Table 
B.5 lists some of the formulas for a three -frequency system, and Table B.6 
for a four -frequency system. 

In using these tables, care must be taken in interpreting the Q^ as reac- 
tive powers when the associated co are negative; see the remarks following 
Eq. 3.26. 

Table B.7 contains a list of some of the devices that are known to obey the 
formulas of each type; the list is, of course, not complete. A few of the more 
important devices that are known not to obey formulas of each type are listed 
in Table B.8. 
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Table B.I. Frequency-power formulas of the four types, under 
general frequency constraints. The formulas hold 
for any values of the parameters 6 a and >^ . 



Type I 

Manley-Rowe Formulas 

(with d-c power) 



Type I 

Manley-Rowe Formulas 

(without d-c power) 



v io 9u a a w a 



8u 



Type II 



Type III (Page form) 



S P ( 1 - cos - 1 > 
a a \ 8w - 



1 - cos S 



> o 



(Pantell form) 



/8w n 

2 P -J > 
a a\8 w 



S P 

d CL 



P g (6 ) > 
6 a v y ~ 
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Table B. 1. (cont. ) 



Type IV Inductive 



a 
Q a 

w 



e, o 



- g a O a ) > 



Type IV Capacitive 



Q / 8w \ 

-5 i . cos a ^ 
<^ a V 3 a a / 



< 



< o 



< o 

- 



where 



rith |(x) ^r nonmcreasing positive fvmction defined from to +. 
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Table B. 2. Frequency- power formulas of the four types, with 

frequencies of the form mo^ + nuj . The formulas of 
types HI and IV hold for all values of B lf 9 , X i and Xo - 



Type I; 

Manley-Rowe Formulas 

(without d-c power) 



mP 



nP 



m n 



- n 

- U 



Type H 



2 2 mQ =0 

m n mn 

S 2 nQ = 

m n mn 



Type HI 



2 S P [1 - cos(m6 i +n9 ] 5:0 
m n mn u 

S S P (mXi+nXo) 2 ^ 
m n mn 

S S P g (61 , )> 
m n mn & mn A 



Type IV Inductive 



Q [1 - cosfmSj +n8 )] 



m n 



S 2 



m n muj } no) 



Type IV Capacitive 






m n 



+nX ) 2 



m n 



nco 



< 



m n 



+ nw 



where 



g mn ( 9 ! , 9 






+ n6 ) U (x) sinCmOj x + n9 x) dx 



where (x) is any nonincreasing positive function defined from 
to +00 . 
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Table B.3. Frequency- power formulas of the four types, with 
frequencies of the form nw . The formulas of 
Types III and IV hold for all values of 9 . 



Type I 

Manley-Rowe formulas 



P 

n n 



Type II 



nQ =0 

IL n 



Type III 



P (1 - cos n6 ) 

% ' 



Z n 2 P 2: 
n n 

Z P g (8) 2:0 
n n & n v ' 



Type IV 
Inductive 



T. (1 _ C os n 8 ) 
n n 

nQ 2=0 
n n 



Type IV 
Capacitive 



(i . C os n 8 ) 



nQ 



where 



oa 

g (9) = n9 \ g(x) sin 
n ^/ 



(n6x)dx 



with (x) any nonincreasing positive function defined from to oo 
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Table B.4. Frequency- power formulas of the four types, with 
frequencies of the form neo + co s ; it is assumed 
that w specifies the explicit time dependence of 
the system 


Type I 
Manley-Rowe Formulas 


- Pa o 


n nto + co g 


Type II 


s a Q a = 


Type in 


s n p a > o 


Type IV 
Inductive 


tt -, n 


n n " p + s ~ 


Type IV 
Capacitive 


n ^ n 


n xuo -f co g 



Table B.5. 



Frequency-power formulas of the four types, with 
only three frequencies: cj g , o> , and 0^=0) - co 
Only a few of the possible formulas of Types fc and 
and IV are given; they hold for all values of the par- 



ameters 9 



8' 






and 



appendix B 

Type I 

Manley-Rowe 

Formulas 

Type II 
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Q s = Q. = - Q p 



Type III 



pJi - cosej -f p [i - cose ] + P [i - cos(e - e s )] > o 



P -f P^ > 



Type IV 
Inductive 



Q p [l - cos9 p ] ^ Qjl cos(e p 6 8 )] 



*P 



. 

- + > o 

0) to 

p 



Type IV 
Capacitive 



Q s [l-cos6 s ] 



Qjl - cos(8 p - 9,)] 



w_ 



w s ~P 

2 r /Y 



0) 



Q s Q - 

_ + < o 

0) W_ 
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Table B. 6. Frequency- power formulas of the four types, with only 
four frequencies: w g , Wp, o)_ = w- - a> g and u> + = w 4- w g . 
Only a few of the possible formulas of types III an IV are 
given; they hold for all values of the parameters 8 S , 9p , 
X s and X p- 



Type I 

Manley-Rowe 

Formulas 



p p 

s 



p 

w , 



0) 03 

P 






Type II 



Q s - Q_ + 

Q + Q + 
p - 



= 
= 



Type HI 



p g [i - cos e g ] +p [i - cos e ] + p_[i-cos(e p -e s )] 

+ P [1 - cos(6 +0 )] > 



P X 2 

s*s 



P 

P + P 



P > 

+ - 



> 

+ - 



P + P + 4P > 
s p - - 

P + P + 4P > 
s p + - 
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Type IV 
Inductive 



Q [1 - cos 9 



] Q [1 - cos 6 ] Qjl - cos(6 p -6 s )] 



Q + [i - cos(6 p + e g )] 

+ _ - - - > 



> 



Q Q Q, 

_ + _I +_I > o 

CO W W 

s - + 
Q Q Q, 



Q Q Q 

_! + E + 4 > 

0) 0) 00 

s p 
Q Q Q, 



CO CO CO 

s p 4 



Type IV ' Q J i . cos 6 J Qj 1 - cos p ] Q_[ 1 - cos(0 p - 8 g )] 
Capacitive 



s . P 

Q [1 - cos{6 +6 )] 



< o 



< 
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Table B.7. Some devices that are known to obey frequency- 
power formulas of each type. The reference given 
is to an earlier section of this book. The exact 
form of the expressions for power at each frequency 
should be consulted before the formulas are used, 
especially for distributed devices 



Type I (Manley-Rowe Formulas) 

6.2 Acoustic system 

4.2 Capacitor, linear 

4.2 Capacitor, nonlinear 
4.4 Capacitor, time-varying 

Circulator 

4.3 Coils, coupled 

6.1 Dielectric, nonlinear 

4.9 Differential gear 

6.1 Electromagnetic medium 

6.3 Electron beam, irrotational 

6.6 Electron beam, rotational 

4.7 Energy conversion device, 

conservative 

6.4 Ferrite 

6.2 Fluid flow, irrotational 
6.6 Fluid flow, rotational 
4.9 Gear 

4.9 Gyrator 

6.4 Gyromagnetic medium 

6.2 Hydrodynamics, irrotational 
6.6 Hydrodynamics, rotational 
4.9 Ideal transformer 

4.6 Induction motor 

4.3 Inductor, linear 

4.3 Inductor, nonlinear 

4.4 Inductor, time -varying 

4.7 Loudspeaker 

6.7 Magnetohydrodynamic system 

4.9 Maser model 

4.5 Mass 

4.5 Mechanical spring 
4.7 Microphone 

Moment of Inertia 

4.10 Network 

4.7 Phonograph cartridge 

6.6 Plasma model 

4.6 Rotating machine model 
Rotational spring 

4.7 Solenoid 

4.5 Spring, mechanical 

4.8 Traditor 

4.9 Transformer, ideal 
4.9 Transformer, physical 

Transmission line 
Waveguide 



Type II 

5.1 Avalanche diode 

5.2 Carbon microphone 

5.2 Cryotron 

5.1 Crystal diode 

5.1 Dashpot 

5.1 Diode 

5.1 Diode, ideal 

5.1 Friction 

5.1 Ideal diode 

5.3 Ideal transformer 

5.4 Network 

5.2 Piezoresistive device 
5.2 Photoconductor 

5.1 p-n-p~n diode 

5.2 Potentiometer 
5.1 Ratchet 

5.1 Rectifier 

5.1 Resistor, linear 

5.1 Resistor, nonlinear 

5.2 Resistor, time-varying 
5.2 Rheostat 

5.1 Semiconductor diode 

5.2 Strain gage 

5.2 Superconductor 

5.2 Switch 

5.2 Thermistor 

5.3 Transformer, ideal 

5.3 Transformer, time -varying 

5.1 Tunnel diode 

5.1 Varistor 
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Table B.7 (continued) 


Type III 




Type IV, Inductive 


5.1 Avalanche diode 


4.3 


Coils, coupled 


5.5 Capacitor, linear 


4.9 


Ideal transformer 


5.2 Carbon microphone 


4.3 


Inductor, linear 


Circulator 


4.3 


Inductor, nonlinear 


5.2 Cryotron 


4.4 


Inductor, time -varying 


5.1 Crystal diode 


4.10 


Network 


5.1 Dashpot 


4.11 


Resistor, linear 


5.1 .Diode 


4.7 


Solenoid 


6.1 Electromagnetic medium, lossy 


4.9 


Transformer, physical 


5.1 Friction 






5.3 Gyrator 






5.1 Hall effect device 






5.1 Ideal diode 






5.3 Ideal transformer 






5.5 Inductor, linear 






5.4 Network 






5.2 Piezoresistive device 






5.2 Photoconductor 






5.2 Potentiometer 






5.1 Ratchet 




Type IV, Capacitive 


5.1 Rectifier 






5.1 Resistor, linear 


4.2 


Capacitor, linear 


5.1 Resistor, nonlinear 


4.2 


Capacitor, nonlinear 


5.2 Resistor, time-varying 


4.4 


Capacitor, time -varying 


5.2 Rheostat 


4.10 


Network 


5.1 Semiconductor diode 


4.11 


Resistor, linear 


5.2 Strain gage 


4.9 


Transformer, ideal 


5.2 Superconductor 






5 . 2 Switch 






5.2 Thermistor 






5.3 Transformer, ideal 






5.3 Transformer, time -varying 






Transmission line 






5.1 Varistor 






Waveguide 
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Table B.8. Some of the raore important systems that are 
known not to obey frequency-power formulas 
of each type 



Type I (Manley-Rowe Formulas) 

Diode 

Friction 

Rectifier 

Resistor 

Switch 

Transformer, time -varying 



Type II 

Capacitor 

Circulator 

Distributed systems 

Gyrator 

Inductor 

Traditor 

Transformer, physical 

Transmission line 

Waveguide 



Type III 

Capacitor, nonlinear 
Capacitor, time -varying 
Distributed systems, nonlinear 

(with exceptions) 
Inductor, nonlinear 
Inductor, time -varying 
Traditor 
Transformer, physical 



Type IV 

Circulator 

Diode 

Distributed systems 

Gyrator 

Resistor, nonlinear 

Resistor, time -varying 

Switch 

Traditor 

Transformer, time-varying 

Transmission line 

Waveguide 



APPENDIX C 



We wish, to show here that for the purposes of using the frequency-power 
formulas of the first type, the Manley-Rowe formulas, it is permissible to 
lump together power at some frequency CD and all its harmonics, and to 
consider it all as power at the fundamental frequency co . 

If we denote by P the power at frequency nco , the contribution to Eq. 
3.36 from frequency co and all its harmonics is 

p pi S P 

v 9 (nco) _ n n9co _ n n 8o> ,^. .. 

11 nto 8co n nto ^a w ^^a 

which is precisely the contribution we would get by considering the power 
S P to be put into the system at the fundamental frequency co . Thus one 
neecfnot, for the purpose of applying these formulas, distinguish between 
power at a frequency and power at its harmonics. 

Note that this is not true for frequency-power formulas of the other three 
types. 
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APPENDIX D 



. Here we demonstrate a fact of some importance in applying frequency-power 
formulas of the first type to rotating machines. Generally there is d-c power 
flow into such machines via the shaft, and a rigorous application of the Manley- 
Rowe formulas requires the use of the d-c terms. However, a rotating shaft 
with an average speed co ought, we feel, to be putting power in at frequency 
co, since that is the rate at which the physical elements pass each other. We 
show here that it is permissible to consider the d-c power input at a shaft as 
power at frequency co , provided the shaft angle cf> * s OELe of the x. variables 
that describe the energy state function (as it generally is). 

Expressing the shaft angle cb in the form of Eq. 3.17, we must equate the 
v. quantity with the average shaft speed co , whether or not the shaft runs at 
constant speed. Thus co is some function of the independent frequencies co , 
and in the rigorous application of the Manley-Rowe formulas the d-c power in- 
put term in Eq. 3.36 would be 

^ 

where P is the average -torque multiplied by co . However, if we were to con- 
sider this power input as occuring at frequency co , the contribution to Eq. 
3.36 from P would be 

fit 

which, being the same as Eq. D.I, is correct. We therefore conclude that no 
error is made, in applying the formulas, if the actual d-c power input at the 
shaft is considered as input at frequency co . 

Note that this problem does not arise for the other three types of frequency- 
power formulas. 



154 



APPENDIX E 



There are many possible sets of independent frequencies one can use in 
writing frequency-power formulas. As indicated in Sec. 3.1 (The Frequen- 
cies), it is desirable to choose the independent frequencies so that as few as 
possible are necessary to specify the explicit time variation of the system. 
Even so, there are many possible choices of the remaining w a independent 
frequencies; frequency-power formulas hold for any such set of independent 
frequencies. In this appendix we determine when new information is obtained 
by selecting a different set. 

Suppose we have two (equally large) sets of independent frequencies on 
which the explicit time dependence of the system does not depend; let one set 
be indexed by a, and the other by b. Then each oj a can be expressed in 
terms of either set, and 

8cx> 8u) 8w a 

__ = 2 i (E.I) 

8u b a 8 a 8u b 

Also, for devices that obey the Manley-Rowe formulas, 

8v. Q 8v. 9w 

10 - S (E.2) 

9w a 8u 9w 

Frequency-power formulas of Type I, using the o> a independent frequencies, 
are 

s ^^o +S ^ (E.3) 

i v iQ 8o a a CO Q 8co a 

If each is multiplied by 9w a /8w , the sum of the formulas becomes 

s 5i !!H + s ^ ^ = (E.4) 

i v iQ 8co b a Wa 9Wb 

or simply the Manley-Rowe formula using c^ as an independent frequency. 
Thus we can get Manley-Rowe formulas using one set of frequencies from those 
using any other set, so no new information is obtained by using new independent 
frequencies. 

The same is true of frequency-power formulas of Type II the same in- 
formation is contained in the formulas using each set of independent frequen- 



cies 



a . 

However, in using frequency-power formulas of Types III and IV, one can 
sometimes get new information from a new choice of independent frequencies. 
This occurs when the formulas used are Eqs. 3.52 or 3.53, the most common 
forms for Type HI, and Eqs. 3.67 or 3.68 for Type IV. However, when one 
uses the more general forms, Eqs. 3.55, 3.56, or 3,58, or Eqs. 3.69, 3.70, 
or 3.71, then no new information is obtained. To conclude this, for example 
for Eq. 3.55, merely put the constants a equal to 
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e a = \ aZT 

D 



for any set of constants 6^ . Then Eq. 3.55 becomes 



a 

S P I - cos S, ~ 6,_ > (E.6) 

a a \ b 9co, b/ - ; 

which is merely the corresponding formula using the o>, set of independent 
frequencies. 

In summary, new information is obtained by a new choice of independent 
frequencies 03 only when using formulas of Types III or IV, and even then 
only when using the simple forms of Eqs. 3.52, 3.53, 3.67, or 3.68. 



APPENDIX F 

We examine here some conditions under which 

V X p = (F.I) 

for systems with fluid flow, where "p is the momentum per particle. We show 
from the equations of motion that if Eq. F.I is ever satisfied, then under cer- 
tain conditions it holds for all time. Flow with this property is called "irro- 
tational; " the equations of motion are simplified and, as shown in Sees. 6. 2 
and 6.3, the Manley-Rowe formulas are more easily obtained. 

For simplicity we discuss at the same time both the acoustic system and the 
electron beam system. To specialize to either case, merely omit the irrele- 
vant terms. To be quite general, we add viscous force to the general force 
equation of Sec. 6.7, Eq. 6.210: 



_*. _ 

V(T + T+^+Q) = -p-FqA + qE + v X u + 9Vs 4- - + -= - + F V 



X B _ 
+ F V 

(F.2) 



where for simplicity we define 

u = V Xp" (F.3) 

The force equation for fluid flow, Eq. 6.42, is obtained by dropping the terms 
involving & , ~A. , IS, B", 9, s, p , and Tr from Eq. F.2; the electron beam 
force equation, Eq. 6.74, is obtained by omitting the terms involving T, TT , 
Q , 9 , s , p , "3fr , and F . 

We take the curl of each side of Eq. F.2, obtaining 

or, using a vector identity for V X (v X H) , 
= "u + (v V)u 

___ p f "E -f J f X B~ 
= (u V)v - u(V * v) + V X F V + V9 x Vs -f V X (F.5) 

If there is no viscous force, and if either the temperature or the entropy is 
uniform, and if there are no free charges or currents, then Eq. F.5 reduces to 

= (u * VJv 7 - ~u(V * v) (F.o) 

For any given particle, this is a set of three coupled linear total differential 
equations, with time-varying coefficients. For given initial conditions, there 
is a unique solution. But by inspection a possible solution is the trivial one 
u = ; therefore if a particle has zero curl of its momentum at any time, it 
always does. It makes sense, then, to talk about M ir rotational flow, because 
if the flow is ever ir rotational, it remains so. 
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On the other hand, if there are viscous forces, or If the entropy and the 
temperature are both nonuniform, or if there are free charges and currents, 
then Eq. F.6 does not hold. These effects spontaneously introduce a curl of 
the momentum; if any of them is present, it does not make sense to speak of 
irrotational flow. 

Eq. F.6 is the differential form of Lagrange's Theorem [28; 101, Sec. 2.3] , 
which is usually written in terms of a surface integral of u being an invari- 
ant of the motion. 

Nonviscous fluid flow is irrotational if it starts from irrotational motion, 
for example if it starts from rest. An electron beam is irrotational if the 
electrons come from a source that is constructed properly. Gabor [28] has 
shown that electrons coming from a cathode have only a normal component of 
H, that is, a component perpendicular to the surface of the cathode, when they 
leave. Furthermore, even this component vanishes if there is no magnetic 
field normal to the cathode. Thus, at any point in space occupied by electrons 
that have come from a cathode with no normal magnetic field, u. = . Any 
electron beam device with a magnetically- shielded cathode therefore has ir- 
rotational flow [28] . 
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